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SUPREM  II  — A PROGRAM  FOR  IC  PROCESS  MODELING  AND  SIMULATION 


Dimitri  A.  Antoniadis,  Stephen  E.  Hansen,  and 
Robert  W.  Dutton 


PREFACE 

SUPREM  II  is  the  second  version  of  the  Stanford  University  IC 
process  simulation  program,  SUPREM.  Several  modifications  and  refine- 
ments relative  to  SUPREM  I [l]  have  been  implemented  but  the  basic 
program  architecture  and  I/O  philosophy  has  not  been  altered.  The 
input  language  syntax  remains  the  same,  but  various  semantic  changes 
have  been  made  mainly  to  accommodate  the  new  or  improved  physical 
models  used  in  SUPREM  II.  Also,  CPU  time  saving  features  have  been 
included,  the  most  important  one  being  the  introduction  of  latent 
regions  in  the  simulated  space  where  the  impurity  distributions  are 
not  processed  until  they  exceed  a physically  significant  level  of 
concentration.  In  this  new  version  of  the  program  the  probability  of 
numerical  instability  such  as  the  occurence  of  negative  impurity  con- 
centrations, has  been  greatly  reduced  by  the  implementation  of  a new 
segregation  model.  Among  new  physical  processes  that  are  now  simulated 
by  SUPREM  II  are,  phosphorus  diffusion  under  high  concentration  and 
the  associated  base  push,  arsenic  clustering  which  affects  both  arsenic 
migration  and  final  electrical  activity,  and  enhanced  oxidation  rates 
due  to  high  surface  concentration  of  n-  and  p-type  impurities. 

This  report  consists  of  three  parts.  In  the  main  part  the  physi- 
cal process  models  and  their  numerical  implementation  in  SUPREM  II 
are  described.  Also,  some  details  of  the  program  and  data  structure 
are  discussed.  The  second  part  (Appendix  1)  contains  the  SUPREM  II 
operating  manual  valid  at  the  time  of  the  program  release.  The  last 
part  of  the  report  (Appendix  2)  consists  of  a set  of  examples  of 
program  usage  in  multi-step  process  simulation  together  with  some  of 
the  resulting  outputs. 


1 . INTRODUCTION 


The  Stanford  University  PRocess  Engineering  Models  (SUPREM)  pro- 
gram is  a computer  simulator  capable  of  simulating  most  typical  IC 
fabrication  steps.  The  program  is  designed  so  that  these  steps  can 
be  simulated  either  individually  or  sequentially,  just  as  they  would 
occur  during  the  actual  fabrication  process.  The  output  of  the  program, 
available  at  the  end  of  each  step,  consists  of  the  one-dimensional 
profiles  of  all  the  dopants  present  in  the  silicon  and  silicon-dioxide 
materials.  These  profiles  may  be  displayed  in  various  formats  including 
line-printer  output,  line-printer  plots,  and  high-resolution  (Calcomp 
type)  plot.  It  is  understood  that,  in  sequential  step  simulation,  the 
output  of  a processing  step  constitutes  the  initial  conditions  for 
the  subsequent  step.  The  junction  depths  and  sheet  resistances  of  all 
n or  p layers  formed  during  the  process  are  also  calculated. 

The  fabrication  step  simulation  is  based  on  several  process  models. 
The  models  implemented  in  SUPREM  are  the  following: 

(a)  ion  implantation 

(b)  chemical  predeposition  through  the  surface  (gaseous  or  solid) 

(c)  oxidation/ drive-in 

(d)  epitaxial  growth 

(e)  etching 

(f)  oxide  deposition 

Migration  of  impurities  is  fully  accounted  for  in  any  of  the  above  models 
involving  high  temperature. 

Communication  between  the  user  and  SUPREM  II  has  been  kept  as  simple 
as  possible.  The  input  file  consists  of  free  format  statements  involving 
key  words  and  numbers.  Typically,  a single  processing  step  can  be  simu- 
lated with  an  input  specification  involving  less  than  60  alphanumeric 
characters.  The  details  of  input  specification  are  given  in  the  SUPREM 
II  user  manual  (Appendix  1) . 

Several  physical  parameters  have  been  stored  in  the  program  and 
constitute  the  default  values  used  by  the  models.  These  values  are 
given  in  Appendix  1 with  the  appropriate  references.  They  may  be 
overriden  by  user-specified  parameters. 
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The  various  process  models  have  been  implemented  in  SUPREM  as  sub- 
programs, each  consisting  of  a number  of  subroutines  and  special 
functions.  Figure  1 is  a schematic  of  this  arrangement.  Input  para- 
meter specifications  have  been  designed  to  resemble  actual  process 
runsheet  data  and  documentation.  A run  may  consist  of  a series  of 
process  steps.  A typical  input  file  is  listed  in  Appendix  2.  The 
sequence  of  steps  and  the  specification  of  the  correct  model  parameters 
are  controlled  by  a supervisor  program  that  evokes  the  appropriate 
step  subprogram.  All  communication  between  the  various  subprograms 
is  directed  through  the  common  variable  area  of  the  computer  memory. 


One  essential  part  of  the  common  area  contains  the  impurity  con- 
centration arrays.  In  the  present  version  of  SUPREM,  there  is  a 
capability  for  handling  up  to  three  different  impurity  species.  The 
impurity  concentration  is  stored  in  terms  of  a discrete  profile  with 
a maximum  number  of  400  points.  Each  concentration  value  corresponds 
to  a point  in  the  discrete  space  (spatial  grid)  defined  along  a 
vertical  axis,  with  its  origin  at  the  surface  of  the  material  — silicon 
(Si)  or  silicon  dioxide  (SiC^)  • Typically,  the  spatial  grid  is  divided 
into  the  SiC^  grid  points,  high-resolution-depth  grid  points,  and  low- 
resolution  grid  points,  as  illustrated  in  Appendix  1.  Within  each  of 
the  above  three  areas,  the  grid-point  separation  is  uniform.  The 
spatial  steps  for  each  area  and  the  three  interface  point  indexes  are 
stored  at  the  end  of  the  concentration  array,  together  with  additional 
information  pertaining  to  the  various  physical  characteristics  (such 
as  dif fusivities , segregation  coefficients,  etc.)  of  the  impurity 
species  represented  in  this  array.  Storage  (and  recovery)  of  interme- 
diate or  final  results  in  nonvolatile  records  is  therefore  greatly 
facilitated. 

Very  often,  during  a processing  step,  the  physical  dimensions  of 
the  simulated  discrete  space  may  change,  as  happens  for  example  during 
oxidation,  etch,  deposition,  and  epitaxy.  The  distance  between 
spatial  grid  points  may  not  be  uniform  during  any  of  the  above  proces- 
ses. For  this  reason,  a cubic  spline  interpolation  routine  has  been 
incorporated  into  the  program  and  is  used  at  the  end  of  each  of  the 
processing  steps  to  restore  the  uniformity  of  the  spatial  grid. 


’ 


2 . PROCESS  MODELS 

2. 1 Ion  Implantation  Model 

The  simplest  description  of  an  implanted  impurity  profile  in 
silicon  or  silicon  dioxide  is  a symmetrical  gaussian  curve  with  first 
two  moments  the  projected  range,  R^,  and  the  standard  deviation,  o^, 
calculated  from  the  LSS  theory  [2].  However,  experimental  distribu- 
tions of  many  ions,  such  as  boron  or  arsenic  are  found  to  be  assymmetrical. 
The  simple  gaussian  approximation  of  those  implanted  profiles  is  often 
inadequate  so  that  higher-order  moments  must  be  used  to  construct  range 
distributions.  Gibbons  and  Mylroie  [ 3 ] have  shown  that  the  third 
central  moment  is  enough  to  provide  sufficient  information  to  construct 
accurate  distributions  when  the  asymmetry  is  not  excessive  (less  than 
the  standard  deviation).  In  these  cases,  the  distribution  can  be  repre- 
sented by  two  half-gaussian  profiles,  each  with  a different  standard 

deviation,  a.  and  CJ„,  joined  together  at  a modal  range  R as  shown 
1 2.  m 

in  Figure  2.  This  method  is  used  to  specify  the  implanted  profile 
of  all  impurities  in  SUPREM  II,  with  the  exception  of  boron  where  a 
modified  Pearson  IV  distribution  is  used  as  described  later. 
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Figure  2 The  joint  half  gaussian  representation  of  as-implanted 
impurity  profiles 

For  the  joint  half-gaussian  distribution,  the  two  sides  are  given  by 


L(y)  = px  exP  [-  (y  - Rm)2/2°12]  0-y-Rn 

] V 


C2(y)  = PX  exp  ■ “ Rm)2/2°22 
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In  SUPREM  II,  the  values  for  R , a,  and  0_  are  obtained  from  the 

mi  l 

Gibbons-Mylroie  algorithm  which  makes  use  of  R^,  o^,  and  of  the  third 
moment  ratio,  y^.  These  last  three  quantities  are  obtained  by  inter- 
polation from  look-up  tables  for  each  element  in  silicon  or  silicon 
dioxide  as  a function  of  implant  energy  [ 4 ].  In  addition  to  the 
acceleration  energy,  the  normal  way  of  specifying  an  implant  is  by 
selecting  the  ion  dose.  As  a result,  a relationship  between  the  dose 
and  peak  concentration  must  be  used.  Defining  the  normalized  ion 
dose  by 


Q(a,R,a,b,)  = / exp  [-  (y  - R)2/2o  ] dy 


The  peak  concentration,  P , and  the  total  dose,Q,  are  related  by 


Py  = Q/[Q(olt  Rm,  0,  Rm)  + Q(o2, 
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In  implantation  through  silicon  dioxide,  two  sets  of  a, , a. , and  R 

i z m 

must  be  used  in  the  two  regions.  Because  the  stopping  power  of  silicon 
dioxide  is  somewhat  different  from  that  of  silicon,  the  modal  range 
in  silicon  is  modified  to  account  for  this  difference  in  the  following 
approximate  way, 

RmSi  = RmSi  + (1  ^Si^ox^  Zox 

I 

where  Z is  the  silicon  dioxide  thickness,  R and  R _.  are  the 
ox  mSi  mSi 

corrected  and  the  uncorrected  values  of  the  silicon  modal  range,  and 
Rpg^  and  RPox  are  the  silicon-dioxide  projected  ranges.  A possible 
small  discontinuity  of  the  profile  at  the  silicon-oxide  interface 


is  an  artifact  produced  by  this  approximation 


Boron  implanted  profiles:  Hofker  et  al  [5]  have  shown  that  the  implan- 
ted boron  profiles,  before  annealing,  may  be  described  by  a Pearson 
type  IV  distribution.  They  have  experimentally  determined  the  fourth 
moment  ratio,  62*  Although  this  distribution  describes  well  the  pro- 
file of  boron  from  the  surface  to  a distance  somewhat  beyond  the  peak 
of  distribution,  Figure  3 shows  that  the  approximation  fails  to  model  the 
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Figure  3 Boron  as-implanted  profiles  in  <11 1>  and  <100>  silicon 
in  a random  direction.  Pearson  IV  and  modified  Pearson 
IV  distributions  for  representation  of  these  profiles. 


observed  exponential  tail  that  is  due  to  random  scattering  of  boron 
ions  along  channeling  directions.  Based  on  experimental  results  [6,7,8] 
we  have  empirically  modified  the  Pearson  IV  distribution  by  adding  an 
exponential  tail  with  a fixed  characteristic  length  (0.045  pm), 
independent  of  dose  energy  and  crystalline  surface  orientation.  The 


tail  is  attached  to  the  shoulder  of  the  standard  Pearson  IV  distribution 


where  the  concentration  drops  to  50%  of  the  peak  value.  Of  course  after 
the  addition  of  the  tail,  renormalization  of  the  distribution  to  the 
implanted  ion  dose  is  carried  out.  Typical  resulting  profiles  from 
this  modification  are  shown  in  Figure  3. 


The  standard  Pearson  IV  distribution  is  calculated  as  follows: 


C(y)  = CQe 
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where 


f(y)  = 4r—  S.n  b + b.  x + b x2 
2b2  o In  2 n 
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a = - Y1(B2  + 3)/A 


bQ  = - (462  - 3Yl)/A 


bl  = 3 


b2  = -(2B2  - 3Yl  - 6) /A 


and  A = 1082  “ 12Y^  - 18 


The  constant  C is  found  from  a normalization  to  the  implanted  dose: 
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2.2  Impurity  migration  during  thermal  processing 


2.2.1  Introduction.  The  redistribution  of  impurities  in  the 
space  of  the  silicon-silicon  dioxide  system,  during  thermal  processing, 
is  governed  by  the  general  continuity  equation  which  can  be  written  as 


dt  f CdV  = y " 4)  dV  ' “ 

V(t)  V(t)  S(t) 


dS 


(1) 


where 

C = impurity  concentration 
S(t)=  closed  surface  (function  of  time,  t) 
V(t)=  volume  enclosed  by  S(t) 
impurity  flux  vector 
outward  unit  normal  to  S(t) 
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g * impurity  generation  rate  per  unit  volume 


l = impurity  loss  rate  per  unit  volume 

The  left  hand  side  of  equation  (1)  represents  the  time  rate  of 
change  of  the  impurity  content  in  any  volume,  V(t),  and  it  is  equated 
to  the  net  impurity  generation  rate  in  the  same  volume  minus  the  net 
impurity  outflow  through  the  surface  S(t)  enclosing  V(t).  The  reason 
for  using  the  continuity  equation  in  integral  form  is  that  it  simpli- 
fies treatment  of  volume  changes  which  occur  during  silicon  oxidation 
and  epitaxy  processes.  The  generation  and  loss  mechanisms  have  been 
included  to  account  for  the  exchange  of  impurity  atoms  between  different 
states  in  the  silicon  lattice  as  in  the  case  of  arsenic  where  atoms  may 
coexist  in  substitutional  and  clustered  states. 

For  one  dimensional  flow,  equation  (1)  can  be  written 


^7  Q(y1»y2)=u(yi>y2)”tF(y2)  ~ F(yi^ 


where 
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C(y)dy 
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^i»y2)  m J (s  - £)dy 


where  y is  the  direction  perpendicular  to  the  silicon  surface  and  is 
positive  inwards,  and  the  flux  F(y)  is  positive  in  the  y direction. 
Physically,  the  impurity  flux  may  arise  from  thermal  diffusion,  from 
interface  phenomena  such  as  evaporation  or  segregation  and  from  the 
motion  of  interfaces  as  in  the  case  of  silicon  oxidation  and  epitaxy. 

In  the  sections  that  follow  we  discuss  the  models  used  to  describe 
each  of  the  physical  processes  present  in  equation  (2)  as  well  as  their 
numerical  implementation  in  the  context  of  SUPREM  II. 
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2.2.2  Diffusive  flux.  Solid  state  diffusion  is  the  physical 
mechanism  responsible  for  impurity  migration  within  the  silicon  body 
during  high-temperature  processing  steps.  At  any  point,  y,  the  diffu- 
sive flux,  FD(y) , of  impurities  is  related  to  their  concentration  and 
diffusivity  gradient  by  the  modified  Fick's  first  law  [9].  For  one - 
dimensional  flow  the  relation  is 

D(y)  * " d7  (D(y)  ' C(y))  (5) 

where  D(y)  is  the  diffusion  coefficient  of  the  impurity.  Under  the 
assumption  of  single-state  migration  there  is  no  generation  or  loss 
within  the  material.  The  conservation  equation,  therefore,  becomes 

dQ(y,  ,y9) 

tz — = - <FD(y2)-  FD(yi))  (6) 


Under  the  further  assumption  of  uniform  diffusivity,  the  well-known 
Fick's  second  law  is  derived  by  using  equations  (3)  and  (5)  in  (6) 
and  differentiating  with  respect  to  y.  The  equation  becomes 


3_C 

at 


(7) 


Fick's  second  law  is  generally  adequate  for  the  calculation  of 
impurity  migration  under  low  impurity  concentration  conditions.  However, 
this  approximation  fails  as  the  impurity  concentration  increases  to  or 
above  the  intrinsic  carrier  concentration,  n^(T) , in  the  semiconductor  at 
the  process  temperature.  Figure  4 is  a plot  of  n^  vs  temperature  as 
defaulted  in  SUPREM  II  [10 ].  In  addition,  equation  (7)  may  fail  even 
at  low  concentration  conditions  if  another  impurity  species  is  present 
in  the  silicon  at  high  concentration.  We  refer  to  silicon  in  which 
any  impurities  exist  at  concentrations  lower  than  n^(T)  at  the  process 
temperature,  as  intrinsic  silicon.  If  the  opposite  is  true  then  we 
refer  to  it  as  extrinsic. 


One  of  the  first  attempts  to  explain  diffusive  flux  under  extrinsic 
conditions  was  to  include  the  "electric  field  effect"  of  the  introduced 
free  carriers  on  the  impurity  ion  migration  [ll],  in  a way  similar  to 
ambipolar  diffusion  in  plasmas.  The  effective  diffusion  coefficient 
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Figure  4 Intrinsic  carrier  concentration  in  silicon  vs  temperature. 
thus  becomes  a function  of  impurity  concentration  given  by 

D - Dife  -D1{l+[l+4  (£)’]  '"j  (8) 

where  is  the  intrinsic  diffusion  coefficient.  As  can  be  seen, 
the  maximum  value  of  f is  2,  for  C » n^  and  this  is  clearly  inade- 
quate to  explain  diffusivity  enhancements  of  the  order  of  10  to  20 
often  observed  with  most  of  the  common  impurities. 

At  present,  it  is  generally  accepted  that  extrinsic  diffusion 
phenomena  are  the  result  of  impurity  migration  by  interaction  with 
charged  point  defects  in  silicon  [12].  All  common  impurities  diffuse 
in  silicon  by  means  of  interaction  with  lattice  point  defects  such 
as  silicon  atom  vacancies  and  possibly  interstitials.  Thus,  the  diffu- 
sion coefficient  is  proportional  to  the  concentration  of  such  point 
defects.  Although  the  concentration  of  neutral  defects  at  any  given 
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temperature  is  dependent  on  the  impurity  concentration  (so  long  as 
it  does  not  approach  that  of  silicon  atoms) , the  concentration  of 
defects  at  various  charge  states  (which  have  been  identified  within 
the  silicon  bandgap) , depends  on  the  Fermi  level  position  in  the  band- 
gap  and  thus  is  a function  of  impurity  concentration.  Recently  Fair 
[13]  summarized  some  of  the  latest  work  on  diffusion  in  silicon  and 
suggested  that  the  effective  diffusion  coefficient  should  be  the  sum 
of  several  dif fusivities,  each  accounting  for  impurity  interactions 
with  different  charge  states  of  lattice  vacancies.  The  effective 
diffusion  coefficient  can  thus  be  expressed  as 

D = fe(DX  + D~  [V~]  + D=  [V=]  + D+  [V+]}  (9) 

v 

where  D is  the  diffusivity  due  to  each  identified  charge  state,  v, 
of  vacancies,  (v:  =,  -,  x,  + i.e  doubly  negative,  singly  negative, 
neutral  and  positive),  and  [Vv]  is  the  concentration  of  vacancies  in 
each  charge  state,  normalized  to  the  intrinsic  concentration  of  that 
state.  Using  Boltzmann  approximation  it  can  be  easily  shown  that  these 
normalized  concentrations  may  be  given  by 

[V_]  = — , [V=]  = (J-)2  and  [V+]  = -i  (10) 

n^  n^  n 

where  n is  the  free  electron  concentration.  Thus,  under  intrinsic 
conditions  (n=n^) , equation  (9)  becomes 

Di  = DV  + D_  + D=  + D+  (11) 

i.e  the  intrinsic  diffusivity  is  the  sum  of  the  diffusivities  resulting 
from  the  various  vacancy  charge  states . 

The  above  model  is  the  basis  of  what  has  been  implemented  in  SUPREM  II. 
However,  with  the  exception  of  phosphorus,  it  has  been  assumed  that  only 
the  neutral  and  one  charged  defect  state  is  responsible  for  the  diffusi- 
vity of  impurity  atoms.  The  specific  form  now  used  in  the  program  is 

D - D±  (1  +6  fy)/(l  + 0)  (12) 

where  is  the  measured  intrinsic  diffusivity  and  f^  = n/n^  for  donors 
and  n^/n  for  acceptors.  Thus,  under  intrinsic  conditions  f *•  1 and 
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D = D^.  On  the  other  hand  for  extrinsic  conditions  the  physical  meaning 
of  the  parameter  8 can  be  derived  by  combining  equations  (9) , (10)  and 
(12) , to  get 


i 1 + 6 


and  D 


i 1 + 8 


Thus,  8 = D /D  , is  an  index  of  the  effectiveness  of  charged  vacancies 
relative  to  neutral  ones  in  impurity  diffusion.  Figure  5 is  a plot  of 
normalized  diffusivity  vs.  f . Although  it  might  be  expected  that  8 
for  any  impurity  element  is  a function  of  temperature,  no  definite  charac- 
terizations exist  at  present.  Recommended  values  (also  included  in  the 
program  as  defaults)  are  8=3  for  boron  and  8 = 100  for  arsenic  while 
for  phosphorus  a completely  different  model  is  used  as  described  later. 

6 is  not  known  for  other  impurities  and  therefore  is  defaulted  to  1. 


10  * 
10 


Figure  5 Normalized  diffusivity  vs  normalized  carrier  concentration 
for  different  values  of  8. 
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It  is  well  known  that  when  different  impurity  atoms  are  present  in 
silicon  there  is  direct  interaction  among  them  [14].  Fig.  6 illustrates 
one  such  (simulated)  case  where  high  concentration  arsenic  effects  the 
distribution  of  boron*.  In  SUPREM  II  this  interaction  is  directly 
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Figure  6 Illustration  of  effects  of  heavy  arsenic  doping  on  the  redistri- 
bution of  boron. 


*The  case  of  phosphorus  is  quite  different.  As  explained  later,  phosphorus 
at  high  surface  concentration  indirectly  affects  the  diffusivity  of  other 
impurities  by  increasing  the  concentration  of  point  defects  in  silicon. 


modeled  through  the  use  of  equation  (5)  for  the  calculation  of  impurity 
flux.  For  the  example  shown  in  Figure  6 the  diffusive  flux  of  boron 
in  the  region  which  is  heavily  n-doped  is  reduced  because  fv (boron)  = 
n^n  ■>  0,  while  at  the  edge  of  the  arsenic  profile  the  dip  is  produced 
by  the  rapid  spacial  change  of  f and  thus  of  diffusivity  for  boron. 

Phosphorus : The  model  implemented  in  SUPREM  for  the  diffusive  migration 
of  phosphorus  has  been  presented  by  Fairand  Tsai  [15].  The  model  predicts 
with  reasonable  accuracy  the  phosphorus  kink  formation  as  well  as  the 
base  push  effect,  commonly  observed  during  heavy  emitter  diffusions  in 
bipolar  technology.  According  to  this  model  the  physical  explanation 
of  these  "anomalous"  effects  lies  in  the  enhancement  of  vacancy  concen- 
tration in  the  silicon  caused  by  the  dissociation  of  the  phosphorus- 
doubly  ionized  vacancy  pairs  that  flow  from  the  surface  into  the  silicon 
bulk.  A typical  high  concentration  phosphorus  profile  is  composed  of 
three  regions  as  shown  schematically  in  Figure  7a. 

(a)  The  surface  region  in  which  phosphorus  diffusion  is  assumed 
to  take  place  by  interaction  with  neutral  and  with  doubly 
charged  vacancies  forming  P+V  pairs.  In  this  region  the 
diffusivity  is  given  by 

D = f !oX  + D =(~r^2l 
e I W-  ) 

where  n,  the  free  electron  concentration  is  calculated  from 
the  total  phosphorus  concentration  using  equation  (48) . Thus 
in  this  region  the  diffusivity  is  monotonically  decreasing, 
as  shown  in  Figure  7a.  The  end  of  the  surface  region  is 
assumed  to  lie  at  the  point  where  the  Fermi  level  drops 
below  .11  eV  from  the  conduction  band  which  is  identified 
as  the  state  level  for  doubly  ionized  vacancies.  This  occurs 
at  a characteristic  electron  concentration,  n , given  by 

ne(T)  = 4.65  x 1021  exp  ( - -eV) 

(b)  The  tail  region  where  the  diffusivity  is  enhanced  relative 
to  the  intrinsic  value  due  to  supersaturation  of  the  silicon 
lattice  by  vacancies  resulting  from  the  P V dissociation. 
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Figure  7a  A typical  phosphorus  doping  profile  with  a 
demarkation  of  the  three  regions  considered 
by  the  Fair  and  Tsai  model.  Also  shown  is 
the  local  diffusivity  profile  derived  by  the 
model . 


In  this  region  the  diffusivity  is  given  by 

. . / 0. 3 eV\  I 

1 + eXP  ' — kT — /J  j (13) 

where  ng  is  the  surface  electron  concentration.  This  equation 
should  only  be  considered  as  an  empirical  fit  to  observations. 

As  can  be  seen  the  tail  diffusivity  depends  very  strongly  on 
ng  and  thus  the  surface  concentration  of  phosphorus.  Since 
the  model  predicts  that  there  exists  a vacancy  defect  super- 
saturation in  the  bulk,  a point  that  has  been  recently  verified 
by  Matsumoto  et  al  [16],  any  other  impurity  elements  present  in 
the  silicon  lattice  beneath  a region  of  heavy  phosphorus  surface 
concentration  should  also  experience  a diffusivity  enhancement 
approximately  equal  to  that  of  phosphorus.  Thus,  in  SUPREM  II 
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an  enhancement  factor,  f 


given  by 


I 


is  calculated  and  multiplies  the  intrinsic  diffusivity,  D , 
of  any  other  impurity  in  the  presence  of  phosphorus.  This 
treatment  permits  the  modeling  of  the  base  push  effect,  but 
also  predicts  an  enhancement  of  diffusivity  of  other  elements 
such  as  arsenic  or  antimony  that  may  constitute  the  buried 
collector  in  bipolar  processes.  This  latter  effect  is  proposed 
by  the  present  authors,  and  has  yet  to  be  experimentally 
confirmed . 

The  beginning  of  the  tail  region  is  determined  as  follows: 

It  has  been  observed  by  means  of  Boltzmann-Matano  analysis  of 
phosphorus  profiles  under  heavy  doping  conditions  [15,17] 
that  there  exists  an  intermediate  region  between  the  surface 
and  the  tail,  where  the  diffusivity  increases  in  propor- 
tion to  ( ) . This  empirical  observation  is  used  in  SUPREM  II 

to  determine  the  beginning  of  the  tail  region.  Namely,  the 
diffusivity  calculated  at  the  end  of  the  surface  region,  i.e. 
at  the  point  where  n = ng,  is  multiplied  by  ^ni^.  as  the 
depth  from  the  surface  increases,  until  the  point  where  the 
value  D ^ is  reached.  At  that  point  the  tail  region  is 
assumed  to  begin. 

It  can  be  easily  seen  from  equation  (13)  that  the  model  just  described 
predicts  that  the  diffusivity  of  phosphorus,  at  least  in  the  tail  region, 
should  increase  monotonically  with  increasing  surface  concentration.  It 
is  however  well  known  that  as  the  surface  concentration  increases  above 

20  -3 

about  3-4  x 10  cm  the  tail  diffusivity  decreases.  This  effect  may 
be  attributed  (Fair  [ 1 8 ] ) , to  silicon  bandgap  narrowing  due  to  the  stress 
Induced  on  the  lattice  by  the  high  phosphorus  concentration.  In  SUPREM  II 
this  stress  effect  is  modeled  by  increasing  the  intrinsic  carrier  con- 
centration, n^,  according  to  the  formula 


n.  = n (unstressed)  exp 
l 1 


(-%) 


where  AE  is  the  stress-induced  bandgap  narrowing  given  by  the 

O 

equation  [18] 


AE  = - 1.5  x 10 
g 


22  r i „ in20  “31  „ 
CTS  “ 3 x 10  cm  J eV 


where  CTg  is  the  total  surface  concentration  of  phosphorus.  Of  course, 
at  the  same  time  ng  must  also  be  increased  by  the  same  amount.  The 
importance  of  including  the  stress  effect  is  illustrated  by  Figure  7b 
where  the  calculated  diffused  phosphorus  profiles  with  and  without 


AE  are  compared  with  experimental  data. 
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Figure  7b  Calculated  phosphorus  profiles  for  two  different  surface 
concentrations.  Slower  diffusion  of  the  higher  surface 
concentration  profile  is  due  to  the  silicon  stress  model. 
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One  additional  stress  effect  also  included  in  SUPREM  II  affects 
the  diffusivity  of  ion  implanted  phosphorus  profiles.  In  this  case 
it  is  observed  that  the  diffusivity  is  further  decreased  under  high 
dose  conditions,  probably  due  to  permanent  lattice  disorder  [18] 
produced  by  the  implantation  process.  Thus,  in  the  case  of  ion 
implanted  phosphorus  an  additional  bandgap  narrowing  is  calculated 
from  the  equation  [18] 

AEg  = - 2.3  x 10_6(Qp)  eV 


where  Q 

P 


is  the  phosphorus  atom  dose. 


2.2.3  Diffusion  under  non-equilibrium  conditions  - Oxidation 
Enhanced  Diffusion.  It  has  been  observed  by  several  authors  [e.g.  19, 

20]  that  the  diffusivity  of  boron  and  phosphorus  is  enhanced  when  the 
silicon  surface  is  oxidized.  This  phenomenon  of  oxidation  enhanced 
diffusion  (OED)  is  generally  attributed  to  enhancement  of  silicon  point 
defects  due  to  the  oxidation.  Perhaps  the  most  plausible  model  has 
been  proposed  by  Hu  [21]  and  subsequently  quantified  by  Antoniadis  et 
al  [20].  This  model  relates  OED  with  oxidation  stacking  fault  (OSF) 
growth,  by  invoking  a dual  diffusion  mechanism  for  impurities  in  sili- 
con whereby  both  vacancy  and  interstitial  defects  are  responsible  for 
diffusion,  and  postulating  an  enhancement,  due  to  oxidation,  of  the 
silicon  self  interstitials.  Thus,  according  to  the  model,  during  oxida- 
tion the  interstitialcy  component  of  impurity  diffusivity  is  enhanced 
leading  to  OED  while  increased  interstitial  precipitation  leads  to  OSF. 

At  present  the  relationship  between  oxidation  rate  and  OED  is  not  known 
with  certainty.  Thus,  for  boron  in  SUPREM  II  temperature  dependent  but 
time  independent  OED  is  assumed  for  an  ambient  causing  silicon  oxidation. 
For  phosphorus  a fixed  enhancement  of  a factor  1.8  is  assumed  for  dry 
O2  and  3.3  for  wet  O2  [18], 


2.2.4  Interfacial  fluxes.  An  interface  is  the  surface  that  sepa- 
rates two  different  homogeneous  regions  such  as  the  silicon  bulk  and 
the  silicon  dioxide  or  the  gas  ambient.  During  thermal  processing  there 
exists  at  least  one  such  interface  that  plays  important  role  in  the 
redistribution  oi  dopant  impurities.  This  may  be  a bare  silicon  surface 


in  an  inert  ambient  where  evaporation  of  impurities  is  taking  place, 
or  the  same  surface  in  a dopant  atom  rich  ambient  where  dopant  incorpo- 
ration takes  place  or  perhaps  the  silicon-silicon  dioxide  interface 
where  dopant  atom  exchange  between  the  two  materials  is  taking  place. 

Although  the  actual  details  of  the  chemistry  taking  place  at  such 
interfaces  may  not  be  known,  generally  the  dopant  atom  flux  across 
interfaces  may  be  phenomenologically  described  by  means  of  a first 
order  kinetic  model  as 

Fs  -MCj  -C2/»eq  l.2)  04) 

where  F is  the  dopant  flux  defined  positive  from  region  1 to  2,  C.  is 
s 1 

the  dopant  concentration  of  the  interface  in  region  1 and  C2  the  same 
in  region  2.  Figure  8 illustrates  the  setup  in  the  case  of  a SiC^/Si 
interface,  m , „ is  the  equilibrium  segregation  coefficient  for  the 

0CJ ) 1 2. 

specific  impurity  species  in  the  system  of  regions  i-2  and  is  defined  as 
c? 

m , y = 
eq,l-2 

Finally,  h which  has  units  of  velocity,  is  the  surface  mass-transfer 
coefficient. 

In  SUPREM  II,  equation  (14),  in  one  form  or  another,  is  used  to 
model  all  interfacial  dopant  flux  phenomena.  Following  is  a brief 
discussion  of  each  of  the  major  cases. 


REGION  1 


REGION  2 


Figure  8 Inter facial  flux,  F , in  the  silicon-silicon  dioxide  system. 

Also  shown  is  the  f?ux  induced  by  the  motion  of  the  interface 
during  silicon  oxidation* 


(a)  Evaporation.  The  mass-transport  coefficient  becomes  the  impurity 

evaporation  coefficient  which  is  a function  of  temperature.  The  same 

evaporation  coefficient  is  used  for  Si02  and  for  silicon.  Also,  is 

assumed  zero  and  m =1. 

eq 


(b)  Chemical  deposition.  Chemical  deposition  is  simply  modeled  by 

assuming  (arbitrarily)  that  h (actually  h is  set  equal  to  1 um/sec), 

m^^  = 1 and  C^  equal  to  either  the  dopant  solid  solubility  or  to  any 
other  (user)  specified  concentration.  Thus,  the  simulated  surface  con- 
centration of  silicon  becomes  very  rapidly  equal  to  C^.  It  is  recognized 
that  the  model  may  be  overly  simplistic  but  there  exist  too  many  different 
processes  by  which  dopants  may  be  deposited  in  silicon,  making  it  impracti- 
cal to  attempt  to  accurately  model  each  of  them  specifically  in  a general 
program  like  SUPREM. 

(c)  Si02  - Si  interfacial  flux.  Under  non-oxidation  conditions  the  Si02~ 

Si  interface  is  stationary  and  equation  (14)  is  sufficient  to  model  the 

impurity  flux  exchanged  between  the  two  regions.  Unfortunately,  there 

exists  practically  no  characterization  of  the  flux  in  this  stationary 

system  with  the  exception  of  the  phosphorus  doped  Si02  - Si  system  (used 

for  silicon  doping) , which  has  been  carefully  explored  by  Ghoshtagore  [22], 

In  this  case  h was  derived  as  a singly  activated  function  of  temperature, 

while  m was  assumed  infinite.  Actually  the  exact  value  of  m would 
eq  J eq 

not  alter  the  observed  result  as  long  as  it  is  kept  large  (say  > 50). 


In  the  case  of  a moving  interface  as  in  silicon  oxidation  there  also 
exists  a mot  ion- induced  interfacial  flux  resulting  from  the  different 
dopant  concentration,  across  the  interface.  This  flux  denoted  by  F^,  is 
given  by 


F 


b 


v 

ox 


(Cx  - aC2) 


(16) 


where  v = dZ  / dt  is  the  oxide  growth  rate  and  a is  the  ratio  of  oxi- 
ox  ox 

dized  silicon  to  resulting  oxide  thickness  (equal  to  0.44).  F^  is  also 

depicted  in  Figure  8.  Generally  this  flux  competes  with  the  flux  Fs • 

If  h » v then  C~  -*■  ni  C,  while  if  h <<  v then  C,  -*■  aC0.  In  all 
ox  2 eq  I ox  I 2 

characterizations  of  the  moving  interface  system  to  date,  the  first 
condition  has  been  implicitly  assumed,  i.e.  that  the  equilibrium 
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segregation  condition  prevails.  Thus,  in  the  absence  of  any  meaningful 

values,  h has  been  arbitrarily  set  equal  to  0.1  vim/min  in  SUPREM  II, 

and  thus  the  condition  h » v is  always  satisfied. 

ox 

Although  the  exact  values  of  the  coefficients  are  not  known  at  present, 
the  above  described  modeling  of  interfacial  fluxes  is  physically  more 
meaningful  than  the  classical  model  that  has  been  used  in  SUPREM  I.  In 
addition  the  new  model  lends  itself  more  to  numerical  evaluation  resulting 
in  substantial  improvement  of  the  solution  stability. 

2.2.5  Generation  and  loss  mechanisms.  Almost  all  dopants  may  exist 
in  silicon  in  more  than  one  state  particularly  when  the  concentrations  are 
high.  Typically  one  of  these  states  is  substitutional  and  therefore  mobile 
while  the  other,  if  present,  may  be  some  form  of  precipitate  or  cluster  which 
is  immobile.  Exchange  between  these  two  states  gives  rise  to  the  generation 
and  loss  terms  in  the  continuity  equation  (1)  of  the  mobile  species.  In 
SUPREM  II  a model  describing  arsenic  clustering  phenomena  [23,24]  is  inclu- 
ded. However,  due  to  lack  of  similar  models,  effects  of  clustering  and 
precipitation  on  the  migration  of  other  impurity  atoms  has  been  neglected. 

The  assumed  chemical  reaction  for  arsenic  is 


where  m is  the  number  of  atoms  in  a cluster  and  k and  k,  the  clustering 

c a 

and  declustering  rate  coefficients.  Defining  the  concentration  of 

clustered  atoms,  C , as 
’ c 

C = C_  - C (18) 

C 1 

where  C^,  is  the  total  concentration  and  C the  substitutional  concentration, 

♦■U  />  AWirn  *-  -1  nn  n/lairtt'  -»  f P TTIOIT  K n T-»  V ■{  ff-  on  O O 


We  can  rewrite  equation  (19)  as 


3 C 

^ = kd  [m(keC)  - Cc]  - 4 - g (21) 

The  above  equation  together  with  equation  (1)  are  used  to  describe  the 
thermal  migration  of  arsenic  atoms  in  silicon.  The  cluster  size  that 
gives  best  fit  to  recent  observed  data  [24  ] is  m = 4 (tetramers) . How- 
ever, there  is  evidence  that  clusters  of  larger  size  may  be  involved. 

The  declustering  rate  coefficient  k,  has  been  adopted  from  the  work  of 

a 

Schwenker  et.  al.  [23].  Under  equilibrium  conditions  (i.e.  3 t = 0) 
equation  (21)  yields 


I=1+4k4c3 


(22) 


This  equation  is  plotted  in  Figure  9.  It  can  be  seen  that  clustering 
imposes  a limit  to  the  maximum  concentration  of  substitutional  and 
therefore  electrically  active,  arsenic  atoms.  This  limit  is  a function 
of  process  temperature. 

Although  the  clustering  model  predicts  the  formation  and  dilution 
of  clusters  during  predeposition  and  subsequent  drive-in,  an  assumption 
must  be  made  about  the  initial  cluster  concentration  when  arsenic  is 
implanted.  In  this  case,  it  is  arbitrarily  assumed  in  SUPREM  II, 
that  the  initial  clustered  arsenic  profile  is  simply  that  determined  by 
thermal  equilibrium  at  the  annealing  temperature,  given  the  known  implanted 
profile  of  total  arsenic  atom  concentration. 


2. 3 Thermal  Oxidation 

The  rate  of  SiC^  growth  on  silicon  is  described  in  SUPREM  II  by  the 
well  known  formula  [25] 

Z2  + AZ  = B(t  + T) 

ox  ox  (2o) 

where  Z is  the  oxide  thickness,  t is  time  and  A and  B are  related  to  the 
ox 

linear  and  parabolic  growth  coefficients  and  Kp  and  normalized  partial 
pressure,  Pq  , of  Oj  by 
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Figure  9 Normalized  total  arsenic  concentration  vs  the  substitutional 
arsenic  concentration,  resulting  from  clustering  model,  as 
a function  of  process  temperature,  under  thermal  equilibrium 


A = po2Vkl 
B = p02kp 

The  parameter  T is  related  to  the  initial  oxide  thickness  by 
zi  (t=0)  + AZ  (t  = 0) 


Under  relatively  low  dopant  concentration  conditions,  Kp  and  depend  only 
on  silicon  crystal  orientation  and  on  the  oxidizing  ambient,  and  they  are 
singly  activated  functions  of  temperature.  All  relevant  data  are  given  in 
Appendix  1 together  with  appropriate  references. 

It  is  however  well  known  that  under  high  surface  concentration  condi- 
tions, such  as  in  MOST  source  and  drain  or  bipolar  emitter  region,  the  oxidation 
rate  of  silicon  is  enhanced.  A detailed  description  of  the  phenomenon  has 
been  given  by  Ho  et  al.  [26].  Briefly,  it  was  found  that  the  linear  rate 
coefficient  is  substantially  enhanced,  indicating  an  enhancement  of  the 
surface  reaction  rate  while  the  parabolic  rate,  which  depends  on  diffusion  of 
02  through  Si02  is  only  modesty  enhanced.  These  enhancements  are  decreased 
with  increasing  temperature.  In  SUPREM  II  we  have  incorporated  a model 
proposed  by  Ho  and  Plummer  [27]  which  relates  quantitatively  the  enhancement 
factors  for  and  Kp  to  the  surface  dopant  concentration.  The  underlying 
idea  is  that  the  rate  of  silicon  oxidation  depends  on  the  lattice  vacancy 
concentration  which  is  affected  by  the  presence  of  dopant  atoms  as  already 
discussed  in  Section  2.2.2. 


According  to  the  model,  the  linear  rate  coefficient  can  be  written  as 
\ = ^ [1  + Y(CT  - D]  (25) 

where  is  the  intrinsic  (i.e.  low  concentration)  coefficient,  y is  an 
experimentally  determined  parameter  given  by 


y = 2.62  x 10  exp 


I"  1.10  evl 

L"  kT  J 


(26) 


CA  the  normalized  total  vacancy  concentration  given  by 

.2 


1 


cT  = 


* •* 


1 + c + c + c 


(27) 


with 


C+  = exp  [ (E+  - E.)/kT] 


; E = 0.35  eV 
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.57  eV 


C”  = exp  C (E±  - E")/kT]  ; E_  = Eg  - 

C=  = exp  [(2E±  - E-  - E=)/kT];  E=  = Eg  - .11  eV 

The  above  expressions  should  be  recognized  as  the  normalized  intrinsic 
concentrations  of  vacancies  in  the  three  charge  states  with  their  cor- 
responding state  energies  in  the  silicon  bandgap.  Finally  the  silicon 

energy  bandgap,  E , and  intrinsic  level,  E . , are  given  as  functions  of 
8 i 

temperature  by 

Eg(T)  = 1.17  - 4.73  x 10~4  [t2/(T  + 636)]  eV 

E.(T)  = E /2  - kT/4 
i g 

For  n-type  dopants  the  enhancement  of  the  parabolic  oxidation  rate 
has  also  been  obtained  by  Ho  and  Plummer  [27].  Kp  is  given  as 

Kp  = Kj  (1  + 6 CT0,22)  (28) 

where 

6 = 9.63  x 10-16  exp  eV  (29) 

is  the  intrinsic  parabolic  rate  and  is  the  total  n-type  dopant 
atom  concentration. 


Since  during  oxidation  the  impurity  surface  concentration  changes 
due  to  diffusion  and  segregation  distribution,  the  calculated  enhanced 
values  of  and  may  generally  be  time- dependent . Thus,  rather  than 
using  the  classical  oxide  growth  equation  (22),  in  SUPREM  II  an  incremen- 
tal form  of  the  same  equation  is  used,  namely 

AZo*  -[!  -<2Zcx  +A)  +V(2Zox  + A)2  + 4B4tJ 

Thus,  as  the  simulation  time  proceeds  in  small  time  increments,  At,  and 

the  impurities  redistribute  in  the  silicon,  the  coefficients  A and  B are 

obtained  from  the  surface  impurity  concentrations  in  each  time  increment, 

and  from  those  values  the  corresponding  increment  of  oxide  thickness, 

Az  is  calculated, 
ox 

2.4  Silicon  Epitaxy 

2.4.1  Introduction.  Epitaxy  is  a processing  step  commonly  used  in 
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modern  semiconductor  fabrication  technologies.  As  with  oxide  growth, 
the  epitaxial  silicon  growth  step  defines  a numerical  moving  boundary 
problem.  The  objective  of  SUPREM  II  is  the  simulation  of  impurity 
redistribution  occuring  during  this  step. 

The  model  for  impurity  redistribution  is  vaguely  similar  to  that 
reported  by  Langer  and  Goldstein  [28].  Although  redistribution  caused 
by  diffusion  and  evaporation  is  well  simulated,  front-side  autodoping 
is  omitted.  Further  investigative  work  is  necessary  before  a reliable 
front  autodoping  model  can  be  incorporated  in  the  process  simulator. 

Backside  autodoping  has  also  been  neglected  to  conform  with  the 
one-dimensionality  imposed  on  SUPREM  II.  This,  however,  should  not  be 
a severe  limitation  because  backside  autodoping  is  often  avoided  either 
by  backside  sealing  or  by  the  use  of  ion-implantation  on  lightly  doped 
substrates.  On  the  other  hand,  frontside  autodoping  may  be  a nonnegli- 
gible  effect,  particularly  when  lightly  doped  layers  are  grown  on 
heavily  doped  ones. 

2.4.2  The  model  Figure  10  illustrates  the  epitaxial  growth  model 
used  in  SUPREM.  It  is  assumed  that  the  bulk  gas  phase  has  a uniform 
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Fb:  MOVING-BOUNDARY  FLUX 

Figure  10  Model  for  the  redistribution  of  impurities  during  epitaxial 
silicon  growth. 


MOVING 

INTERFACE 


Fs:  EVAPORATION  FLUX 


-27- 


i 


equivalent  concentration  of  single  dopant  species,  C*^»  This  is  a 
fictitious  concentration  equal  to  the  desired  uniform  epi-doping  con- 
centration. No  real  details  of  the  gas  phase  are  considered  in  SUPREM  II. 

At  the  solid  gas  interface,  there  exist  two  flux  components,  Fg, 
and  F^  as  indicated  in  Figure  10.  Fg  is  an  impurity  "evaporation"  flux, 

Fs  = h(kC*I  - C].)  (30) 

where  h is  the  impurity  evaporation  coefficient,  as  already  discussed  in 
Section  2.2.4.  As  evaporation  flux,  we  define  any  impurity  flux  exchange 
between  the  gas  and  the  solid,  other  than  direct  moving  boundary  incorpo- 
ration; 0^  is  the  impurity  concentration  at  the  solid  surface,  and  k is 
the  surface  equilibrium  segregation  coefficient  defined  as 


The  term  is  a flux  induced  by  the  interface  motion  and  can  be  calcu- 
lated as  in  the  oxidation  case  (Section  2.2.4).  It  is  given  by 

F,  = v(C*  - C_)  (32) 

b gl  I 

where  v is  the  interface  velocity  (growth  rate) 

Within  the  solid  (silicon)  body,  diffusive  flow  of  impurities  is 
accounted  for  as  has  been  discussed  in  previous  sections. 
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NUMERICAL  IMPLEMENTATION  OF  PROCESS  MODELS 


3. 1 Introduction 

The  previous  section  has  described  the  models  of  the  physical 
processes  that  are  simulated  by  SUPREM.  In  the  present  section,  the 
numerical  implementation  of  those  process  models  will  be  discussed 
with  particular  emphasis  on  the  discrete  formulation  of  the  impurity 
continuity  equation  under  the  moving  boundary  conditions  encountered 
during  silicon  oxidation  and  epitaxy.  However,  in  order  to  establish 
a basis  for  the  discussion  of  this  issue,  the  numerical  solution  of 
the  continuity  equation  under  stationary  conditions  will  be  outlined 
first. 

3.2  Stationary  Boundaries 

The  space  over  which  the  impurity  continuity  equation  (2)  is  to  be 
solved,  is  partitioned  into  discrete  cells.  The  impurity  concentration, 
C(y),  is  evaluated  at  points  (or  nodes)  lying  in  the  middle  of  each  of 
the  discrete  cells.  Figure  11  illustrates  this  space  discretization. 


C(y) 


Figure  11  The  partitioning  of  the  simulation  space  into  discrete 
cells.  Impurity  fluxes  across  cell  boundaries. 
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For  instance,  taking  cell  i = 6 as  an  example,  the  discrete  continuity 
equation  (neglecting  generation-loss  terms)  is  written  as 

dt  Q6  " FD(y5.5^  (33) 

where  FD(y6  5>  and  FD(y5  5>  are  the  impurity  diffusive  fluxes  at  the 
right  and  left-hand  boundaries  of  cell  6,  and  is  the  impurity  con- 
tent of  cell  6.  Thus,  for  the  general  ith  cell  not  lying  at  any  of  the 
space  boundaries,  the  continuity  equation  discretized  in  space  becomes 

ft  «1  - - [v w - V w]  <34> 

where  is  evaluated  at  the  cell  boundary  from  equation  (5)  and  is 
given  by 


For  simplicity  we  have  ignored  generation-loss  terms  in  (34).  Both  the 
spatial  derivative  in  equation  (5)  and  the  integration  in  (35)  are 
carried  out  using  numerical  approximations.  Specifically,  the  flux  F^ 
is  evaluated  by  replacing  (5)  by  a difference  equation  while  Qi  is 
evaluated  using  midpoint  integration. 

The  cells  at  the  two  space  extreme  boundaries  as  well  as  the  SiC^/ 
Si  boundary  (when  it  exists),  deserve  attention  at  this  point. 

(a)  Top  Boundary.  The  first  discrete  cell  is  actually  a half-cell  with 

its  node  at  the  physical  boundary.  Generally,  at  this  boundary  an 

evaporation  or  incorporation  flux,  F (0),  governed  by  equation  (14)  as 

s 

discussed  in  Section  2.2.4,  is  assumed.  At  the  inner  boundary  of  this 
cell  a diffusive  flux  is  evaluated  as  with  all  other  cells,  from  equa- 
tion (5).  The  continuity  equation  for  this  cell  becomes 

■ -[V’W  - Fs<°>]  <36> 


(b)  Deep  Boundary.  This  boundary  usually  lies  inside  the  silicon 
substrate  at  the  point  where  the  simulated  space  terminates.  The  last 


cell  in  this  end  is  also  a half-cell  similar  to  the  first  one.  However, 
a reflecting  boundary  (implied  by  setting  h = 0 in  equation  (14),  is 
assumed  here.  Because  the  depth  of  simulation  is  user  determined,  care 
must  be  taken  to  ensure  that  the  presence  of  a reflecting  boundary  at 
that  point  does  not  affect  the  simulation  results.  On  the  other  hand, 
this  reflecting  boundary  facilitates  significantly  the  use  of  SUPREM  II 
in  unconventional  cases  such  as  Si  on  saphire  or  poly-Si  on  SiC^,  where 
it  may  be  desirable  to  study  the  effect  of  the  reflecting  deep  boundary 
on  impurity  redistribution.  The  continuity  equation  for  this  last  cell 
becomes 


it  V VW  (37) 

(c)  Si02/Si  Interface.  Since  this  interface  is  an  extremum  point  where 

impurity  concentration  must  be  evaluated  it  must  always  lie  on  a node. 

This  node  is  shared  by  two  half-cells,  the  one  in  Si02  and  the  other  in 

Si.  Since  the  impurity  concentration  is  generally  discontinuous  across 

the  interface  (because  of  thermodynamic  segregation) , the  interface  node 

contains  two  different  concentrations,  one  for  the  oxide  half-cell  and 

one  for  the  silicon  half-cell.  An  interfacial  flux,  F , described  by 

s 

equation  (14),  is  assumed  to  flow  between  these  two  cells,  while  diffu- 
sive flux  calculated  from  equation  (5)  with  the  appropriate  dif fusivities 
for  the  two  materials,  is  flowing  across  the  other  two  boundaries.  The 
continuity  equation  for  the  two  interfacial  half-cells  becomes 


— Q 
dt  I, ox 


— Q 
dt  VI, Si 


[Fs  - fd 


[Fo(yi^)  - Fs] 


The  established  boundary  cell  equations  (36),  (37),  (38)  and  (39) 
together  with  equation  (34)  for  all  other  cells  constitute  a system  of 
equations  that  describes  the  temporal  evolution  of  the  impurity  distri- 
bution. Numerically,  this  evolution  is  obtained  by  replacing  the  time 
derivative  in  these  equations  by  a discrete  approximation.  All  the 
above  equations  are  of  the  type 

Hi(t)  ■ h Qi(t)  ' (40) 
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where  H^(t)  denotes  in  abbreviated  form  the  cell  boundary  flux  differ- 
ence. Assuming  that  at  time  t the  concentration  distribution  is  known, 

o 

the  distribution  at  a future  time  t^  may  be  derived  by  solving  the 
equation 


/ 


1 

H±(t)  dt 


o 


W 


(41) 


Various  methods  exist  for  performing  this  integration  numerically. 

The  method  used  in  SUPREM  II  is  first  order  implicit,  i.e.  it  assumes 
that  during  the  time  interval  (t^-tC)),  H (t)  = H^(tp.  Then  the  above 
equation  becomes 

Hi(ti)  = tQi(ti)  - Qi(t0)]/(trto)  (42) 


Now,  there  are  as  many  equations  of  the  form  given  above,  as  the  number 
of  discrete  space  cells.  Also,  since  for  each  cell  the  function  H^, 
is  evaluated  at  t^  (i.e.  at  the  future  time),  it  involves  not  only  the 
unknown  cell  concentration  (^(t.)  but  also  the  two  neighboring  cell 
concentrations  (t^).  Thus,  the  resulting  equations  are  mutually 

coupled  and  form  a system  with  a tridiagonal  determinant.  This  system 
of  equations  is  solved  successively  in  small  (simulated)  time  increments 
by  means  of  gaussian  elimination.  Also,  since  the  system  may  be  non- 
linear due  to  the  dependence  of  diffusivity  on  impurity  concentrations, 
Newton-Raphson  iterations  of  the  solution  are  performed  until  conver- 
gence of  the  results  (concentrations)  is  achieved. 


3.3  Moving  Boundary  - Si02/Si 


When  silicon  is  oxidized  the  interface  between  SiC^  and  Si  moves 
into  the  silicon  and  the  oxide  layer  expands  with  a velocity  v (t) 

OX 

given  by 


v (t) 
ox 


dZ 


_Ci 

dt 


(t) 


(43) 


where  is  determined  as  discussed  in  Section  2.3.  Under  these  condi- 
tions a moving  boundary  induced  impurity  flux,  Fb,  given  by  equation  (16) 


also  flows  across  the  interface. 

The  presence  of  the  moving  boundary  complicates  the  numerical 
formulation  of  the  continuity  equation  in  two  ways.  Firstly,  because 
of  the  non-unity  volumetric  ratio  of  Si  to  Si02  (a  = 0.44),  there  is 
expansion  of  the  discrete  volume  cells  as  they  become  part  of  the 
Si02.  Figure  12  illustrates  this  effect.  Thus,  for  the  cells  around 
the  interface,  volume  is  a function  of  time  and  the  integration  bounda- 
ries of  equation  (35)  change  between  times  tQ  and  t^.  Secondly,  the 
existence  of  the  two  interfacial  fluxes  Fg  and  F^  creates  a jump  dis- 
continuity in  the  impurity  flux,  that  propagates  through  space  during 
oxidation.  Referring  to  Figure  13  the  flux  across  the  interface  from 
time  tQ  to  t^  can  be  considered  as  a superposition  of  two  separate 
processes.  The  first  process  assumes  that  the  interface  moves  instan- 
taneously (at  tQ+) » to  its  new  node  and  that  Fg  flows  for  the  time 
interval  (tQ,t^).  The  second  process  consists  of  the  redistribution 
of  the  impurity  contents  of  the  two  SiC^  cells  nearest  the  interface 
due  to  the  moving  boundary  flux  F^  that  flows  across  the  boundary  of 
these  two  cells. 
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Figure  12  Relationship  of  Si02  and  Si  volumes  when  silicon  is 

oxidized  at  four  different  times.  Here  the  oxidation 
rate  is  assumed  purely  parabolic.  The  tick  marks 
identify  cell  boundaries. 


! 


Numerically  assuming  that  the  discontinuity  occurs  in  the  immediate  vici- 
nity of  the  ith  cell  at  a time  t in  the  interval  (tQ,tj),  the  first  order 
implicit  approximation  to  the  integral  equation  (41)  yields 

1 1~  ^ 

Hi(tl}  1 = [<*!<*!>  - Qi<t  )]/(t  -t  ) (44) 

1 o 

Referring  to  Figure  13,  i may  be  equal  to  I or  to  1-1,  where  I identi- 
fies the  interface  node.  The  + sign  is  used  when  the  equation  represents 
the  continuity  in  the  interfacial  cell  I in  SiC>2  and  the  - sign  when  the 
continuity  in  cell  1-1,  also  in  SiC^,  is  represented. 

To  complete  this  discussion,  it  remains  now  to  outline  the  inter- 
relation of  time  and  space  discretization  that  arises  from  the  considera- 
tion of  the  moving  boundary.  Given  the  constraints  that  the  interface 
must  always  lie  on  a node  and  that  the  number  of  oxide  nodes  may  increase 
by  at  most  one  node  in  any  time  step,  two  basic  possibilities  for  advance- 
ment of  the  interface  exist  in  SUPREM  II,  depending  on  the  oxidation  rate 


Figure  13  Interfacial  cells  and  moving  boundary  flux  at  two  incre- 
ments of  time.  For  simplicity  the  volumetric  ratio,  a, 
is  assumed  unity. 
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which  is  typically  a quadratic  function  of  time.  These  two  possibi- 
lities are  illustrated  in  Figure  14.  Note  that  the  figure  has  been 
drawn  assuming  a = 1 for  clarity.  In  Figure  14(a)  the  number  of 
oxide  nodes  increases  by  one,  while  in  14(b)  it  remains  fixed.  In 
the  first  case,  (a),  the  interface  crosses  the  left  hand  boundary 
of  the  cell  that  contains  the  interface  in  time  t^,  while  in  the 
second  case,  (b),  it  does  not  cross  any  boundary.  When  case  (a) 
arises,  equation  (44)  must  be  used  in  the  two  cells  that  share  the 
crossed  boundary  while  in  case  (b)  the  cells  in  the  neighborhood  of 
the  interface  see  no  special  boundary  flux  i.e.  equation  (42)  is 
valid  just  as  in  all  other  cells.  Of  course,  in  both  cases  the 
fact  that  the  volumes  of  the  two  cells  change  in  the  interval  (tQ,  t^) 
must  be  accounted  for,  in  establishing  the  cell  impurity  content. 

Thus,  as  already  mentioned  the  limits  of  integration  in  equation  (35) 
are  time  dependent  for  the  two  cells  in  the  interface  neighborhood. 


S1O2  Si 


(b) 

Figure  14  Examples  of  interface  motion  (a)  when  the  number  of  oxide 
nodes  increases  by  one  and  (b)  when  the  number  does  not 
increase.  In  case  (a)  the  interface  crosses  the  cell 
boundary  while  in  (b)  it  does  not.  For  simplicity  a = 1 
has  been  assumed. 
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3.4  Moving  Boundary  - Epitaxy 


As  with  oxidation-diffusion,  discussed  in  Section  3.3,  the  epi- 
taxial growth-diffusion  algorithm  has  been  derived  from  the  impurity 
conservation  equation  [Eq.  (1  )].  Special  treatment  of  fluxes  is 
necessary  only  at  the  moving/gas  solid  interface.  Again  the  algorith- 
mic constrain  is  imposed  — the  interface  must  always  lie  on  a spatial 
grid  node  and,  during  any  simulated  time  increment  t,  no  more  than  one 
new  node  of  solid  can  be  generated.  This  is  illustrated  in  Figure  15- 
The  solid  vertical  lines  identify  the  gas/solid  interface,  and  the 
broken  lines  delineate  the  discrete  subcell  boundaries  across  which  the 
fluxes  flow.  The  concentration  in  each  subcell  is  considered  uniform 
with  value  equal  to  the  concentration  at  the  grid  node. 

The  numerical  method  used  for  the  solution  of  the  conservation 
equation  is  implicit,  first  order  in  time  and  uses  midpoint  integration 
in  space.  The  procedure  for  the  solution  of  the  resulting  system  of 
equations  from  one  time  increment  to  the  next  can  be  best  understood  by 
referring  to  Figure  15. 

(a)  At  time  tQ,  the  concentrations  in  the  solid  and 
the  gas  (near  the  interface)  are  known.  These 
concentrations  are  either  the  initial  conditions 
or  the  results  of  the  solution  up  to  the  simulated 
time  tQ(Figure  15a). 

(b)  Still  at  time  tQ,  an  intermediate  step  (Figure 
15b)  is  taken  at  which  the  contents  of  the  two 
cells  near  the  interface  are  rearranged.  Speci- 
fically, the  concentration  at  the  new  node,  i-1, 

to  be  added  during  the  next  time  increment  At, 

* 

is  fixed  at  the  gas  concentration  C^,  and  the 
concentration  at  the  interface  node,  i,  is 
modified  to  account  for  particle  conservation  in 
the  ith  subcell.  The  second  operation  is 
necessary  because  one  nodal  concentration  value 
must  be  used  to  account  for  the  current  in  each 


subcell.  The  results  of  this  intermediate  step  serve 
now  as  initial  conditions  for  the  next  time  increment. 

(c)  Time  is  incremented  to  tQ  + At,  and  the  system  of  coupled 
impurity  conservation  equations  for  each  subcell  is 
solved.  The  interface  has  now  advanced  by  one  grid  node. 

This  is  shown  in  Figure  15c.  From  this  point  on,  the 
cycle  of  operations  is  repeated  as  just  described. 

In  setting  up  the  system  of  equations,  in  step  (c)  above,  the  two 
interface  fluxes,  described  in  Section  3.3.2,  are  also  included; 
however,  although  the  evaporation  flux,  Fs»  flows  out  of  or  into  the 
composite  system  shown  in  Figure  15c,  the  moving  boundary  flux  F 

b 

does  not;  F^  only  results  in  rearrangement  of  the  contents  of  the 
(i-l)th  and  and  ith  cells  in  going  from  the  situation  shown  in 
Figure  15b  to  that  in  15c.  This  is  a key  point  if  particles  are  to 
be  fully  conserved  by  the  algorithm. 

The  equations  for  the  impurity  conservation  in  the  first  two 

cells  from  time  t to  t,  becomes 
o 1 

- + Fs  - i Fb  ’ k “5> 

- + VjW  + 1 Fb  ■ k rQi(ti>  - W3  <46> 

where  Fp  identifies  the  diffusive  flux  and  the  remaining  symbols  and 
indexing  have  been  defined  earlier  in  the  previous  sections. 

3.5  Spatial  Discretization,  Truncation  Errors  and  Time  Step  Generation 

The  impurity  concentration  resulting  from  a diffusion  process  is 
a continuous  function  of  space  and  time.  The  concentration  profile 
resulting  from  numerical  simulation  is  a discrete  approximation  both 
in  space  and  time,  to  this  continuous  function.  This  space  and  time 
discretization  naturally  leads  to  so-called  truncation  errors  of  the 
numerical  solution  of  the  continuity  equation.  These  errors  are  small 
provided  the  space  and  time  intervals  are  small. 
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In  SUPREM  II  the  maximum  number  of  discrete  space  intervals  is 
fixed  to  400  and  thus  there  is  a limitation  imposed  to  the  smallest 
space  increment.  Ay,  that  can  be  used.  However  it  has  been  found 
through  practice  that  spacial  steps  of  order  100X  provide  accurate 
results  under  most  circumstances.  Clearly  there  exists  a trade-off 
between  computation  (CPU)  time  necessary  for  a particular  simulation 
and  spacial  discretization  errors.  It  is  very  hard  to  derive  general 
forms  that  would  rigorously  satisfy  specific  trade-off  goals,  parti- 
cularly since  profiles  of  different  shapes  may  coexist  in  any  simulation. 
A good  rule  of  thumb  for  selecting  Ay  maybe  that  it  be  fine  enough  to 
adequately  represent  the  sharpest  final  profile  of  the  simulation. 

Generally,  as  time  progresses  during  a process  involving  diffusion, 
the  profiles  tend  to  become  less  steep,  and  thus  the  rate  of  change  of 
the  content  of  each  discrete  subcell  diminishes.  It  is  reasonable  to 
assume  that  temporal  discretization  errors  for  a particular  discrete 
space  would  be  proportional  to  the  rate  of  change  of  the  profile,  i.e. 
to  the  rate  of  change  of  impurity  content  of  each  subcell.  Thus,  if  in 
the  interest  of  CPU  time  economy,  temporal  truncation  errors  are  to  be 
kept  rather  constant  during  a simulation,  the  tine  steps  should  become 
longer  as  simulation  time  elapses.  This  scheme  has  been  adopted  in 
SUPREM  II  where  the  time  steps  (which  are  built-in  to  the  program), 
increase  parabolically  with  time  starting  from  a small  value  of  about 
half  a minute.  Again  through  practice  it  has  been  found  that  this  scheme 
provides  accurate  results  under  most  process  simulation  conditions. 

Finally,  it  should  be  mentioned  here  that  a latency  feature  exists 

in  SUPREM  II,  that  can  significantly  reduce  the  CPU  time  of  simulations. 

13  -3 

The  basic  idea  is  that  impurity  concentrations  below  10  cm  (this  is 
arbitrary),  are  not  physically  significant  and  thus  need  not  be  processed. 
An  algorithm  is  built-in  to  the  program  that  keeps  track  of  the  length  of 
simulation  space  for  each  impurity  species  and  increases  it  as  the 
profiles  spread  out  by  diffusion.  Thus,  the  number  of  spacial  points 
over  which  the  continuity  equation  is  numerically  solved,  is  always 
optimum  and  is  controlled  by  direct  feedback  from  the  actual  profile 

- 

being  processed. 


t 
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4. 


CALCULATION  OF  ELECTRICAL  PROPERTIES 


4 . 1 Introduction 

The  previous  sections  have  described  the  models  and  their  numeri- 
cal implementation  used  in  SUPREM  II  to  obtain  the  impurity  profiles 
resulting  from  simulated  processes.  Ideally,  these  profiles  together 
with  some  lateral  geometry  information  would  serve  as  an  input  to  a 
device  analysis  program  that  may  calculate  electrical  parameters  and 
predict  the  behavior  of  devices  under  development , thus  completing  the 
cycle  of  device  computer  design.  However,  a lot  of  insight  into  the 
performance  of  devices  and  circuits  is  often  gained  from  simple  elec- 
trical parameter  calculations.  For  this  reason  SUPREM  II  incorporates 
the  calculation  of  layer  or  sheet  resistance  and  of  MOS  threshold  voltage. 
Following  is  a discussion  of  the  implementation  of  these  two  calculations. 
In  this  section  the  symbol  N(y)  is  used  to  denote  the  net  carrier  con- 
centration under  quasi-neutrality  conditions.  It  is  defined  as  N(y)  = 
C.(y)  - C (y)  where  C stands  of  atomic  impurity  concentration  and  the 
subsripts  A and  D for  acceptor  and  donor  impurity  type. 

4 . 2 Sheet  Resistance  , 

Layer  resistance  is  calculated  as  the  reciprocal  of  the  average  con- 
ductivity of  the  layer  g iven  by  the  equation 

i y? 

— = a = —V  C u(C)  N(y)dy  (47) 

Pq  y2"yl  / 

where  y^,  y ^ define  the  boundaries  of  the  layer  (measured  from  the  top 
silicon  surface),  q is  the  elementary  charge,  N(y)  the  net  carrier 
concentration  defined  above  and  u(C)  the  appropriate  carrier  mobility 
which  is  a function  of  total  impurity  concentration  [20].  The  net 
carrier  concentration  is  calculated  by  assuming  total  impurity  electri- 
cal activity  except  in  the  case  of  arsenic  where  only  the  unclustered 
atoms  are  assumed  to  contribute  electrons  and  in  the  case  of  phosphorus 
where  the  electron  concentration,  n,  is  related  to  the  atomic  impurity 
concentration,  C,  by  the  expression  [l 5 ] 
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C = n + 2.04  x 10  n 


Both  of  the  above  exceptional  cases  are  important  for  heavy  doping 
conditions . 

Sheet  resistance  is  calculated  in  SUPREM  II  for  all  layers  of 
individual  type  resulting  from  a process  simulation.  For  example, 
in  the  case  of  simulation  resulting  into  an  npn  layer  structure  three 
sheet  resistances  are  calculated,  namely,  (a)  of  the  n-layer  extending 
from  the  silicon  surface  to  the  np  junction  (b)  of  the  p-layer  extending 
from  the  first  (np)  to  the  second  (pn)  junction  and  (c)  of  the  n-layer 
extending  from  the  second  (pn)  junction  to  the  end  of  the  simulated 
space.  The  sheet  resistance  of  the  layer  bounded  by  the  end  of  the 
space  may  or  may  not  be  of  physical  significance,  depending  on  the 
particular  process  that  is  simulated.  Also,  the  sheet  resistance  of 
the  middle  layer  in  the  above  example  which  could  be  interpreted  as 
a pinched  base  resistance  may  well  be  an  underestimate  since  depletion 
regions  are  not  considered  in  the  calculations. 

4.3  MOS  Threshold  Voltage 

The  calculation  of  M.0S  threshold  voltage,  V_,  is  implemented  in 

SUPREM  II  based  on  the  full  depletion  approximation  and  the  assumption 

of  quasi-neutral  impurity  profiles  [28].  Figure  16  shows  a typical 

p-type  implanted  impurity  profile.  It  is  assumed  that  all  the  carriers 

are  depleted  for  0<x<x,.  The  threshold  voltage  is  then  calculated 

d 


V = +0)  + V„ 

T FB  vs  B 


where  Vvn  is  the  flat  band  voltage  given  by 


V = <t>  „ - q N /C  (50) 

FB  yMS  4 ss  ox 

with  <J>MS  the  gate-semiconductor  contact  potential  dependent  on  the 

gate  material  and  on  the  channel  surface  doping,  N the  oxide  inter- 
° ss 

face  charge  density  and  C the  oxide  capacitance  per  unit  area. 

ox 

is  the  surface  potential  under  strong  inversion  conditions 


given  by  [28  ] 


W 


Figure  16  A typical  p-type  profile  resulting  from  boron  ion-implan- 
tation. N(x)  is  the  hole  distribution  under  quasi- 
neutrality conditions. 


= + 


kT 


In 


[Wj 


(51) 


where  is  the  carrier  concentration  at  the  bulk  contact  and  N(y^)  the 
same  at  the  edge  of  the  depletion  region.  may  be  the  surface  con- 

centration, the  substrate  concentration  or  an  arbitrary  user  determined 
concentration.  Obviously  and  N(y^)  must  be  of  the  same  type.  Details 
on  the  specifications  of  is  given  in  Appendix  1.  The  positive  sign 
refers  to  p-type  and  the  opposite  to  n-type  impurities. 


Finally,  V is  the  silicon  bulk  voltage  given  by 


V = —3— 
B C 

ox 


y 

f 


N(y)dy 


(52) 


The  one  unknown  in  the  above  equations  is  the  depletion  region  width 
This  is  determined  by  solving  the  Poisson  equation  using  the  full 


depletion  approximation. 

= - f-  (CD(y)  - CA(y))  = ^ N(y) 
dx  s s 


(53) 
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Double  integration  of  this  equation  and  use  of  the  appropriate  boundary 
conditions  yields 


N(y)dy 


-If 


N(y)dy' 


Thus,  y^  is  determined  by  numerical  solution  of  the  non-linear  system 
of  equations  (51)  and  (54). 

Within  the  accuracy  limitations  imposed  by  the  assumptions  used 
here,  can  be  calculated  for  any  arbitrary  impurity  profile.  The 
vertical  profile  under  the  gate  does  not  have  to  be  of  a single  doping 
type.  Thus,  a junction  may  exist  within  the  depletion  region  as  is 
often  the  case  with  implanted  channel  and  depletion  mode  devices.  The 
one  important  consideration  is  that  the  doping  type  at  the  bulk  contact 
be  the  same  as  that  at  the  edge  of  the  depletion  region.  If  the 
profiles  are  such  that  this  is  not  satisfied,  no  threshold  voltage  can 


be  calculated. 


5.  SUPREM  II  DATA  STRUCTURE 

{ 

5 . 1 General  Information 

The  SUPREM  II  data  structure  is  divided  into  five  main  blocks. 

The  first  block  is  made  up  of  six,  450  word,  real  arrays.  The  first 
three  of  these  arrays  contain  information  about  the  impurities  being 
processed  in  the  process  under  simulation.  The  discretized  impurity 
profiles  are  stored  in  the  first  401  words  of  each  array,  followed  by 
one  word  specifying  the  impurity  element,  nine  words  containing  spatial 
grid  information  about  the  profile,  (see  Figure  on  page  55)  and 
the  rest  containing  the  physical  constants  pertaining  to  redistribution 
of  the  impurity.  The  fourth  of  the  six  arrays  is  reserved  for  storing 
the  total  arsenic  distribution  (including  clusters),  while  the  fifth 
is  a general  working  array.  The  sixth  of  the  arrays  in  this  block  is 
used  to  hold  the  400  point  spatial  matrix  indicating  the  distance  bet- 
ween each  of  the  points  in  the  profiles  of  the  previous  arrays.  The 
last  fifty  points  of  the  sixth  array  contain  current  step  parameters 
such  as  temperature,  time,  ambient,  etc. 

The  second  block  is  made  up  of  three  arrays,  one  of  type  logical, 
one  of  type  integer  and  one  of  type  rcai.  These  arrays  contain  infor- 
mation regulating  the  type  and  quantity  of  input  and  output,  along  with 
some  grid  and  substrate  information  that  is  not  normally  variable  from 
step  to  step.  In  addition  some  temporary  grid  pointers  used  when  oxi- 
dizing or  growing  epitaxial  silicon  are  also  stored  in  the  second  block. 

The  third  block  consists  of  arrays  containing  model  card  parameters 
as  described  in  the  SUPREM  II  user’s  manual  (Appendix  1).  The  arrays 
are  initialized  to  the  internal  default  values  when  the  substrate  card 
is  read-in  and  are  modified  by  the  occurence  of  the  model  cards  in  the 
input  deck.  When  a model  card  is  referenced  in  a step  card,  the  model 
parameters  are  read  from  the  appropriate  array  in  this  block. 

The  fourth  block  is  made  up  of  four  410  word  arrays  which  are  used 
as  working  arrays  by  various  routines  throughout  the  program.  The  fifth 
block  is  reserved  for  future  implementation  of  an  optimization  routine. 
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Figure  17  Flowchart  of  SUPREM  II  program. 
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When  SUPREM  is  initiated,  it  first  lists  the  entire  input  deck, 
numbering  the  lines.  Secondly,  it  goes  through  the  deck,  line  by  line, 
checking  the  validity  of  the  syntax  and  card  order.  If,  after  checking 
the  entire  deck,  any  errors  were  found,  SUPREM  lists  the  errors  and  the 
line  numbers  in  which  they  occured,  and  then  terminates.  If  no  errors 
are  found  it  then  begins  processing  the  first  process  in  the  deck. 

This  process  is  illustrated  in  the  flow  chart  of  Figure  17. 

As  described  in  the  user's  manual  the  first  card  in  a processing 
sequence  must  be  a title  card,  this  is  because  when  a title  card  is 
processed  by  the  input  routine,  all  of  the  data  arrays  described  above, 
except  those  of  the  second  block,  are  zeroed.  The  second  block  is 
initialized  to  its  default  values.  The  next  two  cards  in  the  sequence 
must  be  a^gr id  and  a substrate  card.  When  the  first  grid  card  is 
processed,  the  spatial  matrix  is  set  up  in  the  sixth  array  of  the  first 
block,  and  the  spatial  grid  information  in  the  first  of  the  impurity 
arrays  are  set.  The  substrate  card,  by  specifying  the  crystalline 
orientation  of  the  silicon  surface,  causes  the  model  cards  to  be  ini- 
tialized. If  the  substrate  card  also  specifies  a substrate  impurity 
element,  then  the  first  impurity  array  is  opened  and  its  physical  para- 
meters initialized. 

As  impurities  are  introduced  in  subsequent  steps,  the  second  and 
third  of  the  impurity  arrays  in  the  first  block  are  opened  and  ini-' 
tialized.  When  an  element  model  for  one  of  the  impurities  present  is 
referenced  in  a step  card  of  the  process,  the  parameter  values  asso- 
ciated with  that  model  in  the  third  block  are  mapped  into  the  appropriate 
impurity  array.  The  parameters  that  are  associated  with  the  element, 
and  special  purpose  models,  once  set,  are  not  changed  except  by  a later 
reference  to  a model  card  of  the  same  type.  However,  epitaxial  and 
oxidation  model  cards  must  be  specified  at  each  step  in  which  their 
associated  parameters  are  to  be  used,  or  else  the  program  will  use 
its  default  values. 

5.2  Program  Segmentation 

Following  is  a list  of  the  SUPREM  II  main  program  and  segments 
together  with  the  routines  that  they  use.  A brief  description  of  each 
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of  the  routines  is  included  in  the  source  code  and  thus  will  not  be 
given  here.  The  program  structure  reflected  in  this  section  is  the 
original  one,  developed  on  a Hewlett-Packard  2100  minicomputer  under 
DOS-Ill  operating  system.  Although  this  segmentation  is  not  a unique 
solution  it  can  serve  as  guide  to  those  users  of  SUPREM  II  that  may 
wish  to  use  the  program  on  small  machines. 

SUPREM  MAIN  AND  SEGMENTS 

SUPREM  Main 

SUPRM  STPRC  OPTIM  FUNC 
Segment  SUSYN 


SUSYN 

GETLN 

ERSET 

GDSYN 

SBSYN 

PLSYN 

PRSYN 

SLSYN 

MDSYN 

STSYN 

OPSYN 

RSLSY 

STOSY 

ISCOM 

PARSE 

LGTOP 

GETRL 

IGTEL 

MODNN 

STPTP 

PAR  TP 

IFTYP 

ITYPE 

Segment  SUINP 

SUINP 

COMST 

ERSET 

PLTST 

PRNST 

MODST 

STPSV 

OPTST 

RSLTS 

STOPS 

GETLN 

PARSE 

INITE 

ISCOM 

LGTOP 

GETRL 

IGTEL 

MODNN 

STPTP 

PARTP 

IFTYP 

ITYPE 

Segment  TGSIN 

TGSIN 

GETLN 

FGRID 

SUBS 

ISCOM 

ERSET 

PARSE 

BNDRS 

INITE 

OXINT 

LGTOP 

GETRL 

IGTEL 

MODNN 

STPTP 

PARTP 

IFTYP 

TYPE 

Segment  SVLDF 

SVLDF 

PARSE 

SAVFL 

LODFL 

ERSET 

LGTOP 

GETLN 

ISCOM 

GETRL 

IGTEL 

MODNN 

STPTP 

PARTP 

IFTYP 

ITYPE 

Segment  GRDNW 

GRDNW 

PARSE 

ERSET 

CGRID 

LGTOP 

GETLN 

ISCOM 

GETRL 

IGTEL 

MODNN 

STPTP 

PARTP 

IFTYP 

BNDRS 

SPLIN 

SEVAL 

ITYPE 

Segment  STPST 

STPST 

PARSE 

ERSET 

INITE 

MODIN 

OXIDI 

EPTAX 

LGTOP 

GETLN 

ISCOM 

GETRL 

IGTEL 

MODNN 

STPTP 

PARTP 

IFTYP 

OXINT 

SOLUB 

ITYPE 

Segment  IMPLN 


IMPLN 

ERSET 

QDIMP 

IPLNT 

GAUSN 

IMPRS 

PRSON 

TRAPZ 

SPLIN 

SEVAL 

Segment  DIFOX 

DIFOX  PREQ 

OXTHI 

TIMOX 

SITHI 

DIFPR 

OXDEP 

INTRP 

COEF1 

SOLVE 

PRINT 

UPREC 

ACTVE 

DFPB 

DFFY 

PREC 

Segment  DIFEP 

DIFEP 

GRIEP 

PREQ 

CFEPI 

SLVEP 

PREPI 

UPREC 

SPLIN 

SEVAL 

DFFY 

Segment  OPTIO 
OPTIO 

Segment  OPTPR 

OPTPR  DJAC  LNEQS  LUDEC  LUELM 

Segment  OUTP1 

OUTP1  ACTVE  SHEET  VTHRS  QDINT  EOUTP  WRTQD  TRAPZ  SPLIN  INSPL 

Segment  OUTP2 

OUTP2  PROUT  PLTLP  EOUTP  PRVAL  SPLIN  SEVAL 

Segment  SUERR 

SUERR 
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APPENDIX  I 

8UPREM-H  USERS  MANUAL 


STEPHEN  E,  HANSEN  JUIV  04,  1978 

0.  A,  ANT0NIA0IS 
R,  W.  DUTTON 


FOR  THE  0-03  VERSION  OPERATING  ON 
AN  IBM  370/168  USING  THE  FORTRAN-4  LEVEL  H 
COMPILER 


i 


SUPREM-II  INPUT  FORMAT 


THE  INPUT  PARAMETERS  FOR  SUPREM  ARE  SPECIFIED  BY  A SEQUENCE  OF  CAROS  OR 
LINES  IN  A DATA  FILE.  EACH  OF  WHICH  BELONGS  TO  ONE  OF  THE  FOLLOWING  CARO  OR 
LINE  TYPES. 

INITIALIZATION  INPUT/OUTPUT  PROCESS/MODEL 


TITLE 

COMMENT 

GRID 

SUBSTRATE 

END 


PRINT  STEP 

plot  model 

SAVE 

LOAD 


EACH  CARD  CONSISTS  OF  ONE  OF  THE  ABOVE  CARO  TYPE  IDENTIFIERS  FOLLOWED 
BY  EITHER  A PARAMETER  LIST  OR  A CHARACTER  STRING  (OF  WHICH  ONLY  THE  FIRST 
FOUR  CHARACTERS  ARE  SIGNIFICANT).  THE  PARAMETER  LIST  OR  CHARACTER  STRING 
IS  SEPARATED  FROM  THE  CARO  TYPE  IDENTIFIER  BY  EITHER  A COMMA  OR  ONE  OR  MORE 
BLANKS.  ONLY  THE  FIRST  72  COLUMNS  OF  A LINE  ARE  READ  BY  SUPREMIS  INPUT 
PROCESSOR.  IF  NOT  ALL  OF  THE  PARAMETERS  OF  A CARO'S  PARAMETER  LIST  WILL 
FIT  ONE  ONE  LINE,  THEY  MAY  BE  PLACED  ON  THE  FOLLOWING  LINE  IF  A PLUS  (♦)  IS 
USED  AS  THE  first  NON-BLANK  CHARACTER  OF  THAT  LINE.  CARD  TYPES  THAT  USE  A 
CHARACTER  STRING  INSTEAD  OF  A PARAMETER  LIST  MAY  NOT  HAVE  A CONTINUATION 
LINE.  ALL  BLANK  LINES  ARE  IGNORED. 


THE  CARD  TYPES  THAT  HAVE  A CHARACTER  STRING  INSTEAD  OF  A PARAMETER  LIST 
ARF  TITLE,  COMMENT,  STOP  AND  END.  THE  STRINGS  IN  ANY  TITLE  OR  COMMENT  CAROS 
ARF  USED  BY  THE  OUTPUT  ROUTINES,  WHILE  TEXT  IN  STOP  OR  END  CAROS  IS  IGNORED. 
the  PARAMETERS  in  THOSE  CAROS  which  USE  a PARAMETER  list  are  specified  by  a 
PARAMETER  name  TO  WHICH  A VALUE  IS  EQUATED  (I.E.  kNAME»«<V ALUEw) , THE  PARA- 
METER/VALUE PAIRS  MAY  occur  in  ANY  ORDER  witmJN  THE  LIST,  SEPARATED  FROM 
EACH  OTHER  by  commas,  any  blanks  or  spaces  within  THE  LIST  are  IGNORED, 

THE  VALUES  assigned  TO  the  VARIOUS  parameters  may  be  of  several  differ- 
ent TYPES  (I.E.  NUMERICAL,  LOGICAL,  ETC.).  IN  THE  CARO  TYPE  DESCRIPTIONS  ON 
The  following  PAGES,  the  TYPE  OF  VALUE  that  is  ASSIGNED  to  a parameter  is 
indicated  by  one  of  the  following  symbols. 


SYMBOL 

TYPE 

EXAMPLE 

OR  (DESCRIPTION) 

«N» 

NUMERICAL. 

1 0 

. , 2 # 

1.26E-14  , 330.236 

«L> 

logical 

T 

, TRUE  , 

Y , YES  , F , FALSE  , 

<E» 

ELEMENT 

n 

(BORON)  , 

P (PHOSPHORUS),  SB  (A 

<3> 

STEP  TYPE 

(SEE  STEP 

CARO  DESCRIPTION) 

<M» 

MODEL  TYPE 

(SEE  MODEL 

CAPO  DESCRIPTION) 

<P» 

PARAMETER 

NAME 

(THE  NAME  OF  A OPIMJZING  OR 

«PT» 

file  type 

A 

. ASCII 

, B , BINARY 

n , NO 

NTIMONY),  AS  (ARSENIC) 


RESULT  PARAMETER) 
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ALL  LENGTHS,  AREAS,  AND  VOLUMES  ARE  GIVEN  IN  MICRON  BASED  UNITS,  EXCEPT 
FOR  IMPURITY  CONCENTRATIONS,  which  ARE  GIVEN  in  atoms/cubic  CENTIMETER,  ano 
implant  DOSES  GIVEN  IN  ATOMS/SQUARE  CENTIMETER.  TIMES  ARE  ALWAYS  GIVEN  IN 
MINUTES  AND  SO  VELOCITIES  ARE  IN  MICRONS/MINUTE  AND  DIPFUSIVITIES  ARE  IN 
MICRONS  SQUAREO/MINUTE. 

FOLLOWING  IS  A DESCRIPTION  OF  EACH  CARD  TYPE  AND  ITS  ASSOCIATED  PARA- 
METERS. IN  THIS  DESCRIPTION  OPTIONAL  PARAMETERS  OR  PARAMETER  GROUPS  ARE 
ENCLOSED  BV  SQUARE  BRACKETS,  ( ).  WHERE  ONE  OF  A LIST  OF  PARAMETERS  OR 
PARAMETER  GROUPS  MUST  BE  CHOSEN,  THE  CHOICES  ARE  ENCLOSEO  BY  ROUND  BRACKETS, 
( 1,  SEPARATED  BY  THE  WORD  "OR". 


A.  INITIALIZATION  CAROS 


TITLE  CARO 


THE  TITLE  CARO  SPECIFIES  THE  CHARACTER  STRING  USED  ST  THE  PRINT  ANO  PLOT 
ROUTINES  TO  LABEL  THE  OUTPUT.  THE  TITLE  CARO  ALSO  CAUSES  THE  INITIALIZATION 
or  many  or  THE  PROGRAM  > S DATA  BUPPERS  AND  SO  IT  MUST  BE  THE  P IR8T  CARD  IN 
ANT  PROCESSING  SEQUENCE.  IP  SEVERAL  PROCESSES  ARE  TO  BE  RUN,  THEY  MAY  BE 
PLACEO  IN  The  SAME  INPUT  PILE,  SEPARATED  BY  THEIR  TITLE  CAROS, 

TITLE  t«CHABACTER  STRING*) 


COMMENT  CARO 


COMMENT  CAROS  CAN  Bfc  PLACEO  AT  ANY  POINT  IN  A PROCESSING  SEQUENCE  APTER 
the  INITIAL  TITLE,  GRIO  AND  SUBSTRATE  CARDS,  THE  CHARACTER  STRING  SPECIPIED 
IN  THE  LAST  COMMENT  CARD  BEPORE  A STEP  CARO  IS  USED  BY  THE  PRINT  AND  PLOT 
ROUTINES  TO  label  THE  OUTPUT. 

COMMENT  (^CHARACTER  STRING*} 


END  CARO 


THE  END  CARD  IS  USED  TO  TERMINATE  THE  PROGRAM  AND  SO  MUST  BE  THE  LAST 
CARO  IN  THE  INPUT  PILE  OR  CARD  OEC«, 

END  (^CHARACTER  STRING*) 


■ -V 


GRID  C*RO 


THE  SILICON  WAFER  IN  9UPREM  19  REPRESENTED  BY  UP  TO  400  POINTS  SPACED 
ALONG  AN  AXIS  PERPENOICULAR  TO  ITS  THE  SURPACE,  THESE  POINTS  ARE  DIVIOED 
AMONG  THREE  REGIONS,  THE  GRID  SPACING,  (DISTANCE  BETWEEN  ADJACENT  POINTS) 

IS  UNIFORM  WITHIN  EACH  REGION  BUT  USUALLY  DIFFERS  FROM  REGION  TO  REGION. 

THE  FIRST  REGION  IS  THE  ONE  REPRESENTING  THE  OXIDE.  THIS  REGION  IS 
ALLOCATED  FROM  TEN  TO  FIFTY  POINTS  IF  AN  OXIOE  LAVER  EXISTS  AND  ZERO  POINTS 
IF  NO  OXIOE  EXISTS.  THE  GRID  SPACING  AND  THF  ACTUAL  NUMBER  OF  POINTS  USED 
THIS  REGION  ARE  DETERMINED  BY  THE  PROGRAM. 

THE  SILICON  DISTRIBUTION  IS  STORED  IN  the  next  TWO  REGIONS.  THESE 
REGIONS  DIFFER  ONLY  IN  THEIR  GRID  SPACINGS,  THE  SECOND  REGION  HAVING  A GRID 
SPACING  AT  LEAST  TWICE  THAT  OF  THE  FIRST.  THE  GRID  SPACING  IN  THE  FIRST 
SILICON  REGION  IS  DETERMINED  BY  THE  GRID  CARO'S  DYSI  PARAMETER,  ITS  BOUN- 
DRIES  EXTENO  FROM  THE  SILICON  SURFACE  TO  A DEPTH  SPECIFIED  TH|  GRID  CARO'S 
DPTH  PARAMETER.  (THE  SILICON  SURFACE  IS  EITHER  TmE  OX  I OE /S I L ICON  INTERFACE 
OR  THE  SURFACE  OF  THE  WAFER  IE  NO  OXIDE  IS  PRESENT),  THE  GRID  CARD'S  YMAX 
PARAMETER  SPECIFIES  The  maximum  SILICON  DEPTH  REPRESENTED  AND  SO  THE  BOUN. 
DBIES  OF  THE  SECOND  SILICON  REGION  ARE  AT  DPTH,  TmE  END  OF  THE  FIRST  SILICON 
REGION  , AND  AT  V H A X , THE  END  OF  THg  SILICON,  THE  PROGRAM  WILL  TRY  TO  USE  A 
GRID  SPACING  OF  OF  TWO  TIMES  OYSI  IN  THE  SECOND  SILICON  REGION,  BUT  IF  THIS 
WOULD  CAUSE  THE  NUMBER  OF  GRID  POINTS  TO  EXCEED  THE  400  POINT  MAXIMUM,  THE 
PROGRAM  WILL  INCREASE  THE  GRID  SPACING  TO  FIT  THE  NUMBER  OF  POINTS  AVAILABLE. 

THE  REASON  FOR  having  Two  regions  IN  The  SILICON  is  THAT  INCLUOING  a oeep 
STRUCTURE,  SUCH  AS  a BURIED  COLLECTOR,  MAY  REQUIRE  A GRIO  SPACING  THAT  IS  TOO 
COARSE  TO  ACCURATELY  REPRESENT  A STEEP  OP  NARROW  PROFILE  NEAR  THE  SURFACE, 
SUCH  AS  AN  ION  IMPLANT.  IN  THIS  CASE  THE  FIRST  SILICON  REGION,  WHO'S  GRID 
SPACING  CAN  be  SPECIFIED  by  the  USER  thru  THE  OYSI  PARAMETER.  IS  EXTENDED  To 
COVER  THE  STEEP  OR  NARROW  PART  OF  THE  DISTRIBUTION.  IN  ALL  BUT  A FEW  CASES, 
GRID  SPACINGS  OF  LESS  Than  0,01  MICRONS  IN  THIS  REGION  ARE  UNNECESSARY,  IN 
FACT,  GRID  SPACINGS  OF  0.02  AND  0,03  ARE  OFTEN  MORE  THAN  ADEQUATE. 

IF  POSSIBLE,  The  PROGRAM  wILL  USE  LESS  Than  ThE  maximum  number  OF  GRID 
POINTS  FOR  THE  DISTRIBUTION,  AND  SINCE  THE  TIME  OF  CALCULATION  AND  OUTPUT  IS 
DEPENOENT  UPON  THE  NUMBER  OF  POINTS  In  THE  DISTRIBUTION,  THE  USER  SHOULD 
CHOOSE  THE  GRID  PARAMETERS  SO  THAT  THE  MINIMUM  NUMBER  OF  POINTS  ARE  USEO. 

GRID  CAROS  MAY  APPEAR  ANYWHERE  IN  A PROCESSING  SEQUENCE • THIS  ALLOWS 
THE  USER  TO  MODIFY  THE  GRID  SPACINGS  AND  BOUNRIES  TF  NECESSARY  TO  KEEP  SUF- 
FICIENT ACCURACY  WHILE  MINIMIZING  ThE  TIME  OF  CALCULATION,  WHILE  GRIO  CAROS 
MAY  APPEAR  ANYWHERE,  AT  LEAST  ONE  MUST  APPEAR  EITHER  IMMEDIATELY  AFTER  THE 
TITLE  or  SUBSTRATE  CARDS  IN  ORDER  TO  INITIALIZE  THE  GRIO.  THE  TOTAL  DEPTH 
OF  THE  SIMULATION  SPECIFIED  BY  Ymax,  may  BE  EXTENDED  BY  LATER  GRIO  CAROS. 
SILICON  AODED  IN  THIS  WAY  HAS  A CONCENTRATION  IN  THE  ADDED  SILICON  THAT  IS 
EQUAL  TO  THE  CONCENTRATION  IN  THE  LAST  POINT  IN  THE  ORIGINAL  SILICON. 


I«-SURFACE 

I«-Intebfacf 

->I 

I 07 

i o vs  i 

I 

..  T . 

GREATER  than  OR  equal  TO  P*DYSI  } 

i OXIOE 

x 

i SI  REGION 

1 1 

ST  REGION  2 I 

...400  POINTS— 

»•••> 
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GRID  (DYSI»«N»J  (,OPTH»«N»)  t,YMAX»«N>l 


OYSI  I SPACING  BETWEEN  GRIO  POINTS  IN  MICRONS  IN  THE  REGION  FROM  THE 
SILICON  SURFACE  TO  THE  DEPTH  SPECIFIED  BY  DPTH,  OEFAULT  VALUE 
• 0.01 

l OPTH  I THE  OEPTH  IN  MICRONS  OF  THE  SILICON  REGION  THAT  BEGINS  AT  THE 

SILICON  SURFACE  AND  WHOSE  GRID  SPACING  IS  GIVEN  BY  DYSI, 
DEFAULT  VALUE  • 0,01 

YMAX  I THE  MAXIMUM  CALCULATION  THICKNESS  OF  THE  SILICON  IN  MICRONS, 
THERE  IS  NO  OEFAULT  VALUE.  THIS  PARAMETER  IS  REQUIRED  ON  THE 
FIRST  GRID  CARD  IN  THE  PROCESSING  SEQUENCE  AND  CAN  BE  EXTENDED 
SUBSEQUENTLY.  GRID  ALLOCATIONS  BETWEEN  DPTH  AND  YMAX  ARE 
EITHER  2*OYSI  OR  THE  MINIMUM  VALUE  NEEDED  TO  FIT  THE  400  POINT 
LIMIT, 


SUBSTRATE  CARO 


THE  SUBSTRATE  CARD  PROVIDES  THE  PROGRAM  WITH  INFORMATION  ABOUT  THE 
SUBSTRATE.  SUCH  AS  THE  CRYSTALLINE  ORIENTATION  AND  ANY  INITIAL  UNIFORM 
IMPURITY  CONCENTRATION.  SINCE  this  INFORMATION  is  NEEDEO  before  ANY  proces- 
SING  CAN  BEGIN,  SUBSTRATE  CARDS  MUST  APPEAR  IMMEDIATELY  AFTER  EITHER  THE 
TITLE  OR  INITIAL  GRID  CARO,  SUPREM  WILL  CURRENTLY  HANDLE  ONLY  SILICON  SUB- 
STRATES WITH  EITHER  «1U»  OR  w 1 81) * ORIENTATIONS, 

• 

IN  MANY  SIMULATIONS,  THERE  IS  AN  INITIAL  UNIFORM  SUBSTRATE  CONCENTRATION 
THAT  IS  USED  only  TO  PROVIDE  junction  depth  INFORMATION  AND  IS  not  OTHERWISE 
SIGNIFICANT.  when  THIS  IS  THE  CASE,  MUCH  COMPUTATION  TIME  ANO  PRINTED  OUT- 
PUT CAN  BE  SAVED  BY  SPECIFYING  THE  TYPE  OF  INITIAL  SUBSTRATE  IMPURITY  AS 
EITHER  P-  OR  N-TYPE,  (♦  OR  -)  AND  NOT  AS  THE  ACTUAL  ELEMENT  TYPE,  B,  P,  SB 
OR  AS.  IF  USEO  IN  THIS  WAY,  THE  IMPURITY  CONCENTRATION  SPECIFIED  IS  USED  • 
ONLY  IN  the  CALCULATION  of  OUTPUT  INFORMATION  SUCH  AS  THE  JUNCTION  DEPTHS, 
SHEET  RESISTIVITIES,  THRESHOLD  VOLTAGES  etc.  AND  IS  NOT  ACTUALLY  processed 
DURING  HIGH  TEMPERATURE  PROCESSING  STEPS,  PLEASE  NOTE  THAT  THIS  CONCENTRA- 
TION REMAINS  AT  THE  SAME  UNIFORM  LEVEL  THROUGHOUT  THE  EXISTING  SILICON,  EVEN 
INTRINSIC  EPI-LAYERS. 


SUBSTRATE  ORNT»«N>  C,ELEM«<E»  C,CONC»«N>n 

ORNT  I THE  CRYSTALINE  ORIENTATION  OF  THE  wafer.  ONLY  Fill*  or  <180* 
ORIENTATIONS  ARE  ALLOWED. 


ELEM  s THE  TYPE  OF  IMPURITY  INITIALLY  IN  THE  SUBSTRATE  WHOS  UNIFORM 
CONCENTRATION  IS  SPECIFIED  BY  CONC.  MAY  BE  ACTUAL  IMPURITY 
DESIGNATION  SUCH  AS  B,  P,  SB,  OR  AS,  OR  JUST  THE  IMPURITY  TYPE 
SUCH  AS  ♦ FO»  P-TYPE  OR  - FOR  N-TYPE. 

CONC  J The  UNIFORM  IMPURITY  CONCENTRATION  OF  THE  element  specified  by 
THE  ELEM  parameter.  DEFAULT  VALUE  « 0.0 


• - THIS  FEATURE  CAN  SUBSTANTIALLY  REDUCE  CPU  TIME, 
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8.  INPUT/OUTPUT  cards 


THE  INPUT/OUTPUT  CARDS  ARE  THE  PRINT,  PLOT,  SAVE  ANQ  LOAD  CARDS.  THE 
FIRST  TWO  OF  THESE,  THE  PRINT  AND  PLOT  CARDS,  ARE  USED  FOR  CONTROLLING  THE 
HARDCOPY  OUTPUT  OP  THE  PROGRAM,  THE  VALUES  ASSIGNED  TO  THE  PARAMETERS  OF 
THE  PRINT  ANO  PLOT  CARDS  DO  NOT  CHANGE  UNTIL  THEY  ARE  MOOIFIED  IN  SUBSEQUENT 
PRINT  OR  PLOT  CAROS.  THE  OTHER  TWO  INpijT/OUTPUT  CARDS,  THE  SAVE  ANO  LOAD, 
ARE  HSEO  TO  SAVE  THE  CURRENT  STATE  OP  THE  PROCESS  UNDER  SIMULATION  OR  TO 
LOAD  IN  A PREVIOUSLY  SAVED  PROCESS  AND  RESUME  PROCESSING  FROM  THE  POINT  AT 
WHICH  IT  WAS  SAVED. 


PRINT  CARDS 


The  print  Caros  CONTROL  SUPREME  PRINTED  NUMERICAL  OUTPUT,  WHICH  IS  OF 
TWO  BASIC  TYPES.  the  FIRST  type,  CALLED  THE  HEADING  INFORMATION,  IS  CONT- 
ROLLED by  the  print  card's  head  parameter,  the  heaoing  contains  information 

ABOUT  THE  STEP  CARO'S  PARAMETERS  ANO  OTHER  STEP  RELATED  VARIABLES  SUCH  AS 
THE  AMOUNT  OF  OXIOE  OR  SILICON  GROWN  OR  ETCHEO . THE  HEAOING  ALSO  LISTS  THE 
AMOUNT  OF  OXIDE  PRESENT,  THE  TOTAL  SURFACE  CONCENTRATION,  JUNCTION  DEPTHS, 
SHEET  RESISTIVITIES,  TOTAL  AND  INDIVIDUAL  INTEGRATED  IMPURITY  CONCENTRATIONS 
ANO,  IF  THE  GATE  MATERIAL  IS  SPECIFIED,  THE  THRESHOLO  VOLTAGE,  THE  SECOND 
TYPE  OF  PRINTED  NUMERICAL  OUTPUT  IS  A LIST  OF  THE  DEPTH  ANO  CONCENTRATION  AT 
EACH  GRID  POINT  FOR  BOTH  THE  TOTAL  AND/OR  INDIVIDUAL  IMPURITY  DISTRIBUTIONS. 
THESE  ARE  CONTROLLED  BY  THE  PRINT  CARD'S  TOTL  AND  IDIV  PARAMETERS, 

THE  PRINT  CARD  PARAMETERS  ARE  ALL  LOGICAL  FLAGS  AND  SO  MUST  BE  SET  TRUE 
OR  FALSE  BEFORE  THE  OUTPUT  CAN  BE  GENERATED.  SINCE  THE  OUTPUT  ROUTINES  ARE 
CALLED  IMMEDIATELY  after  the  COMPLETION  OF  A PROCESSING  STEP,  IF  ANY  OF  THE 
PRINT  OPTIONS  ARE  to  Be  CHANGED  FOB  THAT  STEP'S  OUTPUT,  THE  PRINT  CARO  THAT 
HAKES  THOSE  CHANGES  MUST  APPEAR  BEFORE  That  STEP  CARO.  PRINT  CARO  OPTIONS, 
flNCE  SET,  apply  FOR  all  SUBSEQUENT  STEPS  UNTIL  CHAN&EO  BY  LATER  PRINT  CARDS. 


PRINT  (HE  AO**L  *]  t,TOTL««L*J  I,IOIV««L*3 

HEAO  | CONTROLS  THE  PRINTING  OF  THE  HEADING  INFORMATION.  DEFAULT 
VALUE  ■ FALSE. 

TOTL  I CONTROLS  THE  PRINTING  OF  THE  NET  IMPURITY  CONCENTRATION  VALUES 
VS,  DEPTH  AT  EACH  GRIO  POINT.  DEFAULT  VALUE  • FALSE. 

IOIV  I CONTROLS  THE  PRINTING  OF  The  INDIVIDUAL  IMPURITY  CONCENTRATION 
VALUES  VS.  OEPTm  AT  EACH  GRID  POINT  FOR  EACH  ELEMENT  PRESENT, 

default  value  ■ false. 
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PLOT  card 


plot  CARDS  CONTROL  the  plots  OP  PIMPURTTY  CONCENTR AT  I OS  VS.  DEPTH,  THEY 
HAY  APPEAR  AT  ANY  POINT  IN  A PROCESSING  SEQUENCE  AFTER  THE  INITIAL  CHID  ANO 

substrate  cardss.  the  standard  plot  output  for  suprem  is  done  on  a line 

PRINTER.  SOME  SYSTEMS  THAT  HAVE  LINE  OR  POINT  PLOTTERS  may  OUTPUT  THE  PLOTS 
TO  THEM  INSTEAD  OF  THE  LINE  PRINTER,  CHECK  hITH  TmE  PERSON  RESPONSIBLE  FOR 
MAINTANENCE  OF  SUPREM,  ONLY  FIVE  OF  THE  PARAMETERS  LISTED  BELOw  ARE  USED 
NHEN  PLOTTING  ON  THE  LINE  PRINTER,  THE  OTHER  ThO,  AXIS  ANO  DATA,  ARE  DES- 
CRIBED IN  CASE  YOUR  SYSTEM  USES  A PLOTTER  WITH  SUPREM, 

The  LINE  PRINTER  plots  ALL  take  up  the  SAME  AMOUNT  OF  SPACE  ON  A PAGE  IN 
ORDER  THAT  DIFFERENT  PLOTS  CAN  BE  OVERLAYED.  THIS  MEANS  IS  THAT  THE  OXIDE 
IS  PLOTTED  IN  A FIXED  AMOUNT  OF  SPACE  INDEPENDENT  OF  THE  AMOUNT  OF  OXIDE. 

IF  NO  OXIDE  EXISTS  THE  IS  STILL  ADJUSTED.  THE  SILICON  IS  ALSO  PLOTTEO  IN 
ThF  SAME  amount  OF  SPACE  ON  Each  PAGE,  HOWEVER  FOUR  TIMES  more  AREA  IS  USED 
Than  FOR  The  oxiOE,  The  amount  OF  SILICON  that  Is  PLOTTEO,  fl.E.  The  PLOT 
WINDOW),  is  SPECIFIED  BY  THfc  PLOT  CARO'S  WIND  PARAMETER,  THE  MAXIMUM  X-AXlS 
value  however  may  be  LARGER  than  the  wINOOw  SINCE  THE  FOUR  MINOR  DIVISIONS 
OF  The  SILICON  AXIS  ARE  LIMITED  TO  WIDTHS  FROM  THE  FOLLOWING  SERIESl  ,01, 
.7.2,  .C25,  .04,  ,eis,  ,PB,  .1,  .2,  ETC.  ThF  Y-AXIS  REPRESENTS  THE  LOG  OF 
ThF  CONCENTRATION  AND  IT'S  MINIMUM  LOG  VALUE  IS  SPECIFIED  BY  THE  CHIN  PARA- 
METER. the  number  of  decades  plotted  is  specified  bv  thf  noec  parameter. 

THE  OTHER  four  plot  CARO  PARAMETERS  are  LOGICAL  FLAGS  WHICH  DETERMINE 

what  information  is  to  be  plotted,  the  totl  parameter  controls  the  plotting 

OF  THE  TOTAL  or  net  CONCENTRATION  AND  THE  JDIV  PARAMETER  CONTROLS  THE  PLOT- 
TING OF  The  INDIVIDUAL  IMPURITY  CONCENTRATIONS,  TH£  AXIS  PARAMETER  CONTROLS 
The  plotting  of  The  AXIS,  whICM  ALLOWS  THE  USER  TO  PLOT  SEVERAL  CURVES  ON 
The  SAME  SHEFT  WITHOUT  REPLOTTInG  THE  AXIS  AT  EACH  STEP,  (THIS  IS  NOT  POS- 
SIBLE ON  MOST  LINE  PRINTERS  SO  THE  AXIS  PARAMETER  IS  IGNORED  BV  THE  LINE 

printer  plot  routine),  tme  data  parameter  controls  the  print  of  some  step 
data  at  the  bottom  of  the  plots,  such  as  the  oxioe  thickness,  the  step  time 

ANO  TEMPERATURE,  (AGAIN  THIS  DOES  NOT  APPLY  TO  LINE  PRINTER  PLOTS), 


PLOT  (wInD»«N»J  f,CMlN»«N>)  f,NOEC«<N»j  (,TOTL««L»)  t,IDIVi«L») 

♦ ( AXIS*«L»)  t,OATA««L*l 

wino  i amount  of  silicon  plotted  in  microns,  default  value  • a. 

C m I n I LOG  OF  the  minimum  CONCENTRATION  plotted,  default  value  • 14 

NOEC  « NUMBER  OF  OECA0E  OF  CONCENTRATION  PLOTTED,  DEFAULT  value  ■ 7 

TOTL  I CONTROLS  the  PLOTING  OF  THE  NET  IMPURITY  CONCENTRATION  VS, 
DEPTH,  DEFAULT  VALUE  * FALSE, 

IOIV  j CONTROLS  The  PLOYING  OF  THE  INDIVIDUAL  IHPuRtTY  CONCENTRATIONS 
VS.  OEPTM  FOR  EACH  ELEMENT  PRESENT,  DEFAULT  VALUE  • FALSE. 

I AXIS  ! CONTROLS  THE  PLOTING  OF  AXIS.  DEFAULT  VALUE  • FALSE, 

DATA  j controls  THE  PLOTING  OF  STEP  DATA  BELOw  THE  PLOTTED  CURVES. 

default  value  « false. 


! 


SAVE  AND  LOAD  CARDS 


the  SAVE  and  LOAD  CARDS  are  used  to  STORE  AND  READ  the  INFORMATION 
NECESSARY  TO  RESTART  A PROCESSING  SEQUENCE.  THEIR  MAIN  USE  IS  FOR  SIMULA- 
TION RUN  SPLITS,  IN  WHICH  FOR  EXAMPLE,  SEVERAL  EMITTER  DIFFUSION  TIMES  MAY 
BE  SIMULATED  WITHOUT  RERUNNING  THE  ENTIRE  PROCESS  EOR  EACH  DIFFUSION  TIME. 
THE  SAVE  CARDS  MAY  APPEAR  AT  ANY  POINT  IN  THE  PROCESSING  SEQUENCE  AFTER  THE 
INITIAL  TITLE,  grid  and  SUBSTRATE  CARDS,  THE  LOAD  CARD  SHOULD  APPEAR  IMME- 
DIATELY AFTER  THE  INITIAL  TITLE,  GRIO  AND  SUBSTRATE  CARDS. 

THE  SAVE  and  LOAD  CARDS  HAVE  ONLY  TWO  PARAMETERS.  ONE  SPECIFIES  THE 
FORTRAN  LOGICAL  UNIT  NUMBER  THAT  IS  ASSIGNED  TO  THE  DATA  FILE  IN  THE  JOB 
CONTROL  LANGUAGE  OECK  WHILE  THE  OTHER  SPECIFIES  THE  TYPE  OF  OATA  TRANSFER. 


SAVE  LUNM»<N»  , TYPE««FT* 
LOAD  LUNM*<N»  , TYPE«<FT> 


LUNN  | A FORTRAN  ITO  LOGICAL  UNIT  NUMBER  WHICH  HAS  BEEN  ASSIGNED  TO 
A DATA  FILE  IN  THE  J.C.L  DECK, 


TYPE  J DESCRIBES  THE  TYPE  OF  OATA  TRANSFER  AS  EITHER  ASCII  OR  BINARY 

type,  information  transfer  in  binary  form  is  about  twenty 
FASTER  Than  ASCII  ANO  THE  DISC  STORAGE  IS  MUCH  LESS. 
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C.  PROCESS/MODEL  CAROS 


STEP  CARO 


STEP  CAROS  ARE  USED  TO  SPECIPY  THE  PROCESSING  STEP  TYPES  AND  PARAMETERS. 
THERE  ARE  CURRENTLY  SIX  DIFFERENT  types  of  PROCESSING  steps  that  can  be 
MODELED  by  SUPREMOl  ION  IMPLANTATION,  PREDEPOSITION,  OXIDATION,  EPITAXIAL 
GROWTH,  LOW  TEMPERATURE  oxide  deposition  ano  etching,  drive-ins  are  modeled 
BY  OXIDATIONS  IN  NITROGEN  OR  NEUTRAL  AMBIENTS.  EACH  STEP  IN  A PROCESSING 
SEQUENCE  USES  THE  IMPURITY  DISTRIBUTIONS  THAT  RESULTED  FROM  THE  PREVIOUS 
STEP  AS  THE  STARTING  POINT  FOR  ITS  CALCULATIONS. 

EACH  TYPE  OF  STEP  HAS  IT'S  OWN  PARAMETER  LIST,  THOUGH  THEY  MAY  USE  MANY 
OF  THE  SAME  PARAMETERS.  THE  VALUE  OF  ANY  STEP  PARAMETER  THAT  IS  NOT  SPEC- 
IFIED IS  SET  TO  ZERO,  WITH  THE  EXCEPTION  OF  TEMPERATURE,  WHICH  WILL  USE  THE 
LAST  SPECIFIED  TEMPERATURE.  FOLLOWING  IS  A DESCRIPTION  OF  EACH  STEP  TYPE 
AND  IT'S  ASSOCIATED  PARAMETER  LIST. 


ION  IMPLANTATION 


ION  IMPLANT  STEPS  may  BE  SPECIFIED  IN  TWO  WAYSI 

(1)  IF  THE  IMPLANT  ENERGY  IS  SPECIFIED  THEN  THE  PROGRAM  USES  INTERNALLY 
STORED  INFORMATION  TO  CALCULATE  THE  PROFILE,  IN  THIS  CASE  THE  PROFILE  IS  A 
TWO-SIDED  GAUSSIAN  FOR  THE  ELEMENTS,  ARSENIC,  PHOSPHORUS  AND  ANTIMONY,  AND  A 
MODIFIED  PEARSON  TYPE-IV  DISTRIBUTION  FOR  BORON. 

tZ)  IF  THE  RANGE  AND  STANDARD  OEVATION  ARE  SPECIFIED  THEN  3UPREM  USES 
THESE  VALUES  TO  CALCULATE  A SIMPLE  GAUSSIAN  DISTRIBUTION  FOR  THE  ELEMENT. 


STEP  TYPE-IMPL,  ELFM»«F>,  DOSE»«N*.  (AKEV»«N>)  OR  CRANG««N*,  STDV««N») 

♦ t,MOOL"«H*I 

ELEM  | THE  IMPURITY  ELEMENT  TO  BE  IMPLANTED  CB,  P,  SB,  AS), 

DOSE  I THE  IMPLANT  DOSE  IN  ATOMS/SQUARE  CENTIMETER. 

AKEV  « THE  IMPLANT  ENERGY  IN  KEV, 

RANG  « THE  RANGE  of  the  SIMPLE  GAUSSIAN  DISTRIBUTION  IN  MICRONS. 

STDV  I THE  STANDARD  DEVIATION  OF  THE  SIMPLE  GAUSSIAN  DISTRIBUTION  IN 
MICRONS, 

MOOL  I the  name  OF  an  ELEMENT  or  spfcial  purpose  MODEL, 
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ETCH  STEP 


I 


ETCHES  IN  SUPREM  CAN  8E  EITHER  HIGH  OR  LOW  TEMPERATURE.  A LOW  TEMPERA- 
TURE ETCH  IS  ONE  IN  WHICH  THE  TEMPERATURE  IS  BELOW  2P0  DEGREES  CENTIGRADE. 

IN  LOW  TEMPERATURE  ETCHES  THE  AMOUNT  OP  OXIDE  ETCHED  AWAT  IS  EQUAL  TO  THE 
ETCH  RATE  TIMES  THE  STEP  TIME.  ANT  LOW  TEMPERATURE  ETCH  STOPS  WMEN  THE  ALL 
OF  THE  OXIDE  HAS  BEEN  REMOVED.  IE  EITHER  THE  TIME  OR  ETCH  RATE  PARAMETER  IS 
ZERO  OR  NOT  SPECIFIED  THEN  ALL  OF  THE  OXIDE  IS  REMOVEO. 

IN  HIGH  TEMPERATURE  ETCHES.  THE  OXIDE  IS  ASSUMED  TO  BE  ETCHED  AWAY  IN  A 
NEGLIGABLE  amount  OF  time,  AND  SO  The  AMOUNT  of  SILICON  REMOVED  18  equal  TO 
THE  step  time  TIHES  The  ETCH  RATE.  IN  THE  CURRENT  VERSION  of  SUPREM  the 
IMPURITIES  in  THE  REMAINING  material  IS  not  REDISTRIBUTED,  I,E‘.  IS  ASSUMED 
TO  BE  NEGLIGABLE. 


STEP  TYPE«ETCH,  TEMPbwN*  t,TIME««N»,  ERTE»«N>)  t,MOOL»«M*J 

TEMP  ! THE  SILICON  TEMPERATURE  DURING  THE  STEP  IN  DEGREES  CENTIGRADE. 
TIME  I THE  TOTAL  STEP  TIME  IN  MINUTES. 

ERTE  t THE  ETCH  RATE  IN  MICRONS/MINUTE. 

MOOL  I THE  NAME  OF  AN  ELEMENT  OR  SPECIAL  PURPOSE  MODEL. 


LOW  TEMPERATURE  OXIDE  DEPOSITION  STEP 


IN  TMIS  STEP,  DOPED  OR  UNDOPEO  OXIDE  IS  DEPOSITED  ON  THE  SILICON  SURFACE 
OR  ON  TOP  OF  ANY  EXISTING  OXIDE,  THE  AMOUNT  OF  OXIDE  DEPOSITED  IS  EQUAL  TO 

the  STEP  TIME  times  The  growth  rate,  IF  the  deposited  oxide  is  to  be  doped 
then  THE  IMPURITY  is  specified  by  the  elem  parameter  while  conc  specifies 

IT'S  UNIFORM  IMPURITY  CONCENTRATION, 


STEP  TYPE  »DEPO  , T IME  »«N* , GRTE»«N>  l,ELEM««E*,  CONCxNk)  l,MODL»<M»J 
TIME  I THE  TOTAL  STEP  time  IN  MINUTES. 

GBTE  « THE  OXIOE  GROWTH  RATE  IN  M ICRONS /M I NUTE . 

ELEM  | THE  OXIOE  DOPING  IMPURITY  (B,  P,  SB,  AS). 

conc  i the  impurity  concentration  in  the  deposited  oxide  later  in 
ATOMS/CUBIC  CENTIMETER. 

MOOL  I THE  name  OF  AN  ELEMENT  or  special  purpose  model. 
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OXIDATION  STEP  (AND  DRIVE-IN) 


THE  OXIDATION  STEP  IS  A HIGH  TEMPERATURE  STEP  WHICH  INCLUDES  THE  REDIS- 
TRIBUTION OP  THE  IMPURITIES  PRESENT  DUE  TO  DIFFUSION  AND  EVAPORATION.  IN 
ADDITION  IF  THE  AMBIENT  INDICATED  BY  THE  OXIDATION  MODEL  SPECIFIED  BY  THE 
MODL  PARAMETER  IS  NOT  A NEUTRAL  ONE,  THEN  IMPURITY  REDISTRIBUTION  DUE  TO  THE 
GROWTH  OF  OXIDE  ALSO  OCCURS.  THE  TEMP  PARAMETER  SPECIFIES  THE  WAFER  TEMPER- 
ATURE AT  THE  BEGINNING  OF  THE  STEP.  THE  FINAL  TEMPERATURE  IS  EQUAL  TO  THE 
SPECIFIED  INITIAL  TEMPERATURE  PLUS,  THE  TEMPERATURE  RATE  OF  CHANGE  SPECIFIED 
BY  TRTE,  TIMES  THE  TOTAL  STEP  TIME.  THE  OXIDATION  AMBIENT  IS  SPECIFIED  BY 
AN  OXIDATION  MODEL  CARD  REFERENCE  (SEE  OXIDATION  MODELS).  OTHER  MODEL  CARDS 
MAY  ALSO  BE  REFERENCED  IN  THE  SAME  CARO  BUT  ONLY  ONE  OXIDATION  MODEL  REFER- 
ENCED IS  ALLOWED.  IF  NO  OX  10 AT  ION  MODEL  IS  REFERENCED,  THEN  A NITB  MODEL  IS 
ASSUMED. 


STEP  TYPE  aOXIO , T1ME»«N»  [,TEMP»«N»J  l,TRTE»«N>]  (,MODL»<M>) 

TIME  S THE  total  STEP  TIME  IN  MINUTES. 

TEMP  J THE  WAFER  TEMPERATURE  AT  THE  START  OF  THE  STEP  IN  DECREES 

CENTIGRADE.  DEFAULT  VALUE  ■ LAST  SPECIFIED  STEP  TEMPERATURE. 

TRTE  J THE  TEMPERATURE  RATE  OF  CHANGE  IN  DEGREES/MINUTE.  DEFAULT 
VALUE  ■ ». 

MODL  : THE  NAME  OF  AN  OXIDATION,  ELEMENT  OR  SPECIAL  PURPOSE  MODEL. 


~ — — n 


PREDEPOSITION  STEP  (NON-IMPLANT) 


THE  PREDEPOSITION  STEP  MODELS  THE  INRODUCTION  OF  AN  IMPURITY  ELEMENT 
FROM  * CONSTANT  SOURCE  AT  THE  SILICON  SURFACE.  SUPREM  DOES  NOT  CARE  ABOUT 
OETAILS  OF  THE  IMPURITY  SOURCE,  THE  CDNC  PARAMETER  SPECIFIES  THE  SOURCE 
CONCENTRATION  AT  THE  SURFACE  AND  THE  SOURCE  ELEMENT  TYPE  IS  INDICATED  BY 
THE  ELEM  PARAMETER,  THE  SILICON  TEMPERATURE  AT  THE  START  OF  THE  STEP  IS 
SPECIFIED  BY  THE  TEMP  PARAMETER  MHILE  THE  FINAL  TEMPERATURE  IS  EQUAL  TO  THE 

INITIAL  TEMPERATURE,  plus  the  temperature-rate-of-change,  specified  by  trte 
times  The  total  step  time. 


STEP  TYPE«POEP,  ElEM««E  * , CONC»«N*  TIMEayNa  f,TEMP»«N»J  f,TRTE««N*J 
♦ I,M0DL"«M»3 

ELEM  i THE  IMPURITY  ELEMENT  BEING  INTRODUCED  (B,  P,  SB,  AS), 

CONC  I THE  surface  gas  CONCENTRATION  OF  THE  IMPURITY  SPECIFIED  by 
ELEM  IN  ATOMS/CUBIC  CENTIMETER, 

TIME  I THE  TOTAL  STEP  TIME  IN  MINUTES. 

TEMP  | THE  WAFER  TEMPERATURE  AT  THE  BEGINNING  OF  THE  STEP  IN  DEGREES 
CENTIGRADE.  DEFAULT  VALUE  « LAST  SPECIFIED  STEP  TEMPERATURE, 

TRTE  I THE  TEMPERATURE  PATE  OF  CHANGE  IN  DEGREES  CENT IGR AOE /M INUTE , 
DEFAULT  VALUE  • B, 


WOOL  I THE  NAME  DF  AN  ELEMENT  OR  SPECIAL  PURPOSE  MOOEL 


EPITAXIAL  GROWTH  STEP 


THIS  STEP  MODELS  THE  EPITAXIAL  GROWTH  OF  A SILICON  LAYER  AT  HIGH  TEMPER- 
ATURES, ON  TOP  OF  AN  INITIAL  LATER  OF  SILICON,  THE  ADOEO  SILICON  MAT  DOPED 
OR  INTRINSIC,  IF  IT  IS  TO  BE  OOPED,  THE  IMPURITY  ELEMENT  IS  SPECIFIED  BY  THE 
ELEM  PARAMETER  AND  THE  BULK  GAS  CONCENTRATION  (NOT  THE  SURFACE)  IS  SPECIFIED 
BY  THE  CONC  PARAMETER.  THE  AMOUNT  OF  SILICON  GROWN  IS  EQUAL  TO  THE  GROWTH 
RATE,  SPECIFIED  BY  GRTE,  TIMES  THE  TOTAL  STEP  TIME.  THE  WAFER  TEMPERATURE 
AT  THE  BEGINNING  OF  THF  STEP  IS  SPECIFIED  BY  THE  TEMP  PARAMETER.  THE  FINAL 
TEMPERATURE  IS  the  INITIAL  TEMPERATURE  PLUS,  the  TEMPERATURE-RATE-OF-CHANGE 
SPECIFIED  BY  TRTE,  TIMES  THE  TOTAL  STEP  TIME.  THE  MODL  PARAMETER  MAY  BE 
USEO  TO  SPECIFY  ANY  DESIRED  element,  EPITAXY  or  special  PURPOSE  MODELS. 


STEP  TYPE«EPIT,  TIMEw<N>,  GRTEwwNP  (,TEMP»«N>)  ( , TRTEwwNwJ  (,MOOL»«M>) 

* (ELEM««E»,  CONC»«N>) 

TIME  S ThE  TOTAL  STEP  T I MF  IN  MINUTES. 

GRTE  J THE  GROWTH  RATE  OF  THE  EPI-LAYER  IN  MICRONS/MINUTE, 

TEMP  i THF  WAFER  TEMPERATURE  AT  THE  8EGINNIN  OF  THE  STEP  IN  DEGREES 
CENTIGRADE.  DEFAULT  VALUE  • LAST  SPECIFIED  STEP  TEMPERATURE. 

TRTE  : THE  TEMPERATURE  RATE  OF  CHANGE  IN  DEGREES/MINUTE.  DEFAULT 
VALUE  ■ P. 

MOOL  S THE  NAME  OF  A ELEMENT,  EPITAXY  OR  SPECIAL  PURPOSE  MODELS, 

TION  MODEL. 


ELEM  J THE  IMPURITY  ELEMENT  IN  THE  GROWN  SILICON  (B,  P,  SB,  AS). 


MODEL  CARD 


THE  MODEL  CAROS  ALLOW  THE  USER  TO  MODIFY  INTERNAL  PROGRAM  COEFFICIENTS. 
THERE  ARE  FOUR  MAIN  TYPES  OF  MODEL  CAROS  AND  TWO  OF  THESE,  ELEMENT  AND 
OXIDATION,  ARE  OIVIDEO  INTO  FOUR  SUBTYPES.  EACH  OF  THESE  TYPES  OR  SUBTYPES 
ARE  IDENTIFIED  BY  A UNIQUE  THREE  CHARACTER  NAME.  IN  AODITION  THERE  ARE 
FIVE  OF  EACH  OF  THE  MODEL  TYPES  OR  SUBTYPES,  IDENTIFIED  BY  A SINGLE  DIGIT 
NUMBER  FOLLOWING  THE  THREE  CHARACTER  IDENTIFIER.  THE  MODEL  TYPES,  SUBTYPES 
AND  THEIR  NAMES  ARE  LISTED  BELOW. 


TYPE 

SUBTYPE 

NAME 

(«M») 

ELEMENT 

MODELS 

(BORON} 

MBOtt 

* 

• 

1 

m 

S 

(PHOSPHORUS) 

MPH* 

« 

■ 

1 

m 

3 

(ANTIMONY) 

MSB# 

« 

■ 

l 

m 

3 

(ARSENIC) 

HASP 

« 

■ 

1 

m 

3 

OXIDATION  MODELS 

(STEAM) 

STM* 

« 

■ 

1 

m 

3 

(WET) 

WET# 

« 

■ 

1 

m 

5 

(DRY) 

DRY* 

« 

■ 

1 

m 

3 

(NITROGEN) 

NIT* 

N 

■ 

1 

m 

3 

EPITAXY 

MODELS 

mmm 

EPI* 

« 

■ 

I 

m 

3 

SPECIAL 

PURPOSE 

MODELS  --- 

SPM* 

« 

■ 

1 

m 

5 

A PARTICULAR  MOOFL  IS  REFERENCED  IN  A STEP  CARD  RY  ASSIGNING  THE  MOOEL 
Name  TO  A MODL  PARAMETER  (I.E,  M0DL»MA92).  THIS  RESULTS  IN  THE  PARAMETER 
VALUES  SPECIFIED  ON  THE  MODEL  CARO  OF  THAT  NAME  BEING  ASSIGNED  TO  THE  COR- 
RESPONDING PROGRAM  VARIABLES.  THESE  N£W  PARAMETER  VALUES  WILL  BE  USED  IN 
ALl  SUBSEQUENT  STEPS  UNLESS  REASSIGNED  BY  ANOTHER  REFERENCE  TO  A MODEL  CARD 
OF  THE  SAME  TYPE  D«  SUBTYPE,  IF  A STEP  REFERENCES  A MODEL  WITH  THE  NUMBER 
ZERO  AS  THE  IDENTIFYING  DIGIT,  THEN  ALL  OF  THE  PROGRAM  VARIABLES  ASSOCIATED 
WITH  MODELS  DF  THAT  TYPE  OR  SUBTYPE  ARE  RESET  TO  THEIR  OFFAULT  VALUES, 

FOLLOWING  IS  A DESCRIPTION  OF  EACH  OF  THF  MODEL  TYPES  AND  THEIR  PARA- 
METER LISTS  t 


ELEMENT  MODELS 


the  ELEMENT  MOOELS  ALLOW  the  user  to  change  those  PARAMETERS  which  are 
ELEMENT  specific,  there  are  four  ELEMENT  MODEL  subtypes  corresponding  to 

each  OF  THE  ALLOWABLE  IMPURITY  ELEMENTS!  BORON,  PHOSPHORUS,  ANTIMONY  ANO 
ARSENIC.  The  ELEMENT  THAT  THE  MODEL  CARD'S  PARAMETERS  APPLY  TO  IS  SPECIFIEO 
BY  THE  NAME  PARAMETER  ASSIGNMENT,  THE  ELEMENT  ANO  THE  CORRESPOND ING  MODEL 
NAMES  ARES 


BORON 

! 

MBOt , 

MB02, 

MB03. 

MB04, 

HBOS 

PHOSPHORUS 

1 

MPHJ, 

MPH5!, 

MPH3, 

MPH4, 

MPH5 

ANTIMONY 

1 

MSB1  , 

MSB?, 

MSB3, 

MSB4, 

M3B5 

ARSENIC 

! 

H AS  J , 

MAS1, 

MAS3, 

MAS4, 

MASS 

THE  COEFFICIENTS  THAT  CAN  BE  ACCESSED  THROUGH  THE  ELEMENT  MODEL  CAROS 
ARE  LISTED  BELOW,  ALONG  WITH  THE  WAY  THE  MODEL  CARO'S  PARAMETERS  ARE  USED 


TO  CALCULATE  THEM, 

OXIDE  DIFFUSION  COEFFICIENT  ■ DOXP 
SEGREGATION  COEFFICIENT  • 3EG0 

surface  transport  coefficient  ■ stcb 

MOVING  BOUNORY  FLUX  TERM  • HSFP 

SILICON  diffusion  coefficient  a 

fNON-QXIDIZINGl  • DSXN 

(FOR  BORON  IN  DRY  OXYGEN)  a DSXN 
(FOR  BORON  IN  WET  OXYGEN)  ■ DSXN 

(ALL  OTHERS  IN  DRY  OXYGEN)  a DSXD 
(ALL  OTHERS  IN  WET  OXYGEN)  « DSXW 


* EXP  (-EOXP/KT) 

* EXP  (-SEGE/KT) 

* EXP(-STCEZKT) 

* EXP (-HSFE/KT) 

* tXP(-ESXN/KT) 

* EXP (-ESXN/KT)  ♦ DSXD  * EXP  (-ESXD/KT) 

* EXP(-ESXN/KT)  ♦ DSXW  • EXP(-ESXW/«T) 

a EXP (-ESXD/KT) 

* EXP (-E3XW/KT) 


WHERE  KT  ■ BOLTZMAN'S  constant  a the  TEMPERATURE  IN  KELVIN  DEGREES, 


THE  PARAMETER  BETA  IS  USED  TO  MODIFY  THE  SILICON  DIFFUSION  COEFFICIENT 
OF  IMPURITIES  BY  A MULTIPLICATIVE  FACTOR  CALCULATED  BY, 

FACTOR  ■ ( 1 ♦ BETA  * (N  OR  P)/NIEF  )/(l  * BETA) 

WHERE  (N  OR  P)  IS  THE  FREE  CARRIER  CONCENTRATION  AT  A GRID 
POINT  AT  THE  PROCESS  TEMPERATURE,  AND  NIEF  IS  THE 
INTRINSIC  CARRIER  CONCENTRATION  AT  THE  PROCESS  TEMPER- 
ATURE. 

IF  THE  IMPURITY  CONCENTRATION  IS  LESS  THAN  IPX  OF  NJEF  THEN  THE  NUMBER 
CARRIERS  EQUALS  NIEF  AND  THUS  FACTOR  ■ 1. 
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THE  CLUSTERING  OF  ARSENIC  AT  HIGH  CONCENTRATION  LEVELS  IS  CONTROLLED  BV 
THE  FOLLOWING  THREE  PARAMETERS. 

(EQUILIBRIUM  CLUSTERING  COEFFICIENT)  ■ KEQ0*EXP (-KEQE/KT) 

(DECLUSTERING  COEFFICIENT)  ■ KDCBaEXP  (-KOCE/KT) 

(NUMBER  OF  ATOMS  PER  CLUSTER)  a CL3T 


MODEL  NAME»«M»  (,nOX0KN»  t,EOX0a«N»J  ( i DSXNaaN*)  (,ESXNa<N») 

« (,OSXOa«N»)  ( ( ESXD«<N»]  (,OSXWB«N»)  (,ESXMa<N») 

* ( , SEGB  *<N»]  ( , SEGEaaNa]  (,HSF0b«n>)  ( , H8FE«<N») 

* ( , STC0a<N»J  (,STCEa«N>)  (, BET  Aa«N»]  (,CLSTa«N*) 

* ( t KEQB  a<N»]  t,«EQE»«N»l  [ , K DC0 ■ «N*)  (,KOCEa«N») 

DOXP  I PRE-EXPONENTIAL  TERM  OF  THE  SILICON  DIOXIDE  DIFFUSION  COEF- 
FICIENT in  MICRONS  SQUARED/MINUTE, 

EOX0  I ACTIVATION  ENERGY  OF  TMl  SILICON  OIOXIDE  DIFFUSION  COEFFICIENT 
IN  ELECTRON  VOLTS, 

DSXN  l PRE-EXPONENTIAL  TERM  OF  THE  INTRINSIC  SILICON  DIFFUSION  COEF- 
FICIENT IN  MICRONS  SQUARED/MINUTE. 

DSXN  I ACTIVATION  ENERGY  OF  INTRINSIC  SILICON  DIFFUSION  COEFFICIENT 
IN  ELECTRON  volts. 

OSXD  I PRE-EXPONENTIAL  term  of  TMF  SILICON  DIFFUSION  DRY  OXIDATION 
ENHANCEMENT  COEFFICIENT  in  MICRONS  SQUARED/MINUTE, 

E3X0  I ACTIVATION  ENERGY  OF  THE  SILICON  DIFFUSION  DRY  OXIDATION 
ENHANCEMENT  COEFFICIENT  IN  ELECTRON  VOLTS, 

SEG0  I PRE-EXPONFNTIAL  TERM  OF  THE  SEGREGATION  COEFFICIENT.  DEFINED 
THE  CONCENTRATION  IN  THE  SILICON  OVER  THE  CONCENTRATION  IN 
THE  OXIDE, 

SEGE  I THE  ACTIVATION  ENERGY  OF  THE  SEGREGATION  COEFFICIENT  IN 
ELECTRON  volts. 

STC0  I THF.  PRE-EXPONENTIAL  CONSTANT  of  the  surface  TRANSPORT  COEF- 
FICIENT IN  MICRONS/MINUTE, 

STCE  I THE  ACTIVATION  ENERGY  OF  THE  SURFACE  TRANSPORT  COEFFICIENT 
IN  ELECTRON  VOLTS. 

HSF0  1 PRE-EXPONENTIAL  CONSTANT  of  the  MOVING  80UN0RY  FLUX  in 
MICRONS/MINUTE. 

HSFE  I PRE-EXPONENTIAL  CONSTANT  OF  THE  MOVING  BOUNDRY  FLUX  IN 

electron  volts. 

BETA  | COEFFICIENT  RELATING  THE  EFFECT  OF  UNCHARGED  VACANCIES  ON  THE 
DIFFUSION  COEFFICIENT. 

CLST  i NUMBER  OF  IMPuRITYS  ATOMS  PER  CLUSTER, 
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KEQ0  | PRE»E*P0NENTIAL  CONSTANT  or  THE  EQUILIBRIUM  CLUSTERING  COEF 
FICIENT. 

KEQE  I ACTIVATION  ENERGY  OF  THE  EQUILIBRIUM  CLUSTERING  COEFFICIENT 
KDC0  | PRE-EXPONENTIAL  CONSTANT  OF  the  DECLUSTEPING  COEFFICIENT. 
KOCE  : ACTIVATION  ENERGY  OF  THE  OECLUSTERING  coefficient. 


OXIDATION  MODELS 


THE  OXIDATION  MOOElS  HAVE  TWO  JUNCTIONS;  1)  TO  SPECIFY  THE  TYPE  OF 
AMBIENT  THAT  THE  STEP  IS  BEING  CARRIED  OUT  IN,  AND,  2)  TO  MODIFY  THE  OXIDE 
GROWTH  RA7E  PARAMETERS  FOR  THAT  AMBIENT,  THE  OXIDATION  MODEL  NAMES  ALLOW 
DISTINCTION  BETWEEN  FOUR  AMBIENTS,  THE  ambient  TYPES  and  THE  CORRESPONDING 
MODEL  NAMES  ARE  LISTED  BELOW, 

STEAM  I STM  1 , STM2,  STM3,  STMA,  STM5 

WET  02  I WET  1 , WET2,  WETS,  WET4,  WET5 

DRY  02  I DR Y 1 , DRY2,  DRY3,  DRY4,  DRY3 

NITROGEN  I DRY  1 , DRY2,  OR Y3 , DRY4,  DRY5 

THE  COEFFICIENTS  THAT  ARE  CAN  MODIFIED  8Y  THE  OXIDATION  MODEL  PARAMETERS 
ARE  THE  LINEAR  AND  PARABOLIC  OXIDE  GROWTH  RATES,  THE  EXPRESSIONS  FOR  THESE 
COEFFICIENTS  USING  THE  OXIDATION  MODEL  PARAMETERS  ARE  I 

LINEAR  GROWTH  RATE  ■ LRTE  * EXP  (-LREA/KT)  • PRES 

PARABOLIC  GROWTH  RATE  ■ PRTE  * EXP  (-PREA/KT)  * PRES 

WHERE  KT  ■ BOLTZMANS  CONSTANT  * TPHE  TEMPERATURE  IN  KELVIN  DEGREES, 

MODEL  NAME»«M»  [,LRTE»«N>]  [,LBEA«<N*J  t,PRTE«<N>]  (,PREA»<N>) 

♦ t , PRES««N>J 

LRTE  « PRE-EXPONENTIAL  TERM  FOR  CALCULATING  THE  LINEAR  GROWTH  RATE  IN 
MICRONS/MINUTE, 

LREA  t ACTIVATION  ENERGY  FOR  CALCULATING  THE  LINEAR  GROWTH  RATE  IN 
ELECTRON  VOLTS, 

PRTE  I PRE-EXPONENTIAL  TERM  FOR  CALCULATING  THE  PARABOLIC  GROWTH 
RATE  IN  MICRONS  SQUARED/MINUTE . 

PREA  I ACTIVATION  ENERGY  FOR  CALCULATING  THE  PARABOLIC  GROWTH  RATE  IN 
ELECTRON  volts. 

PRES  I THE  AMBIENT  PRESSURE  IN  ATMOSPHERES, 


■ 4.A  


EPITAXY  MODELS 


THE  EPITAXY  MODELS  ARE  USED  TO  MODIFY  THE  IMPURITY  PLUX  PARAMETERS  USED 
BY  THE  EPITAXIAL  GROWTH  ROUTINE,  THERE  ARE  FIVE  EPI  MOOELS,  NUMBERED  ONE 
THROUGH  FIVEI  EPI1,  EPI2,  EPI3,  EPI4,  AND  FPI5.  THE  PARAMETER  LIST  OF  THE 
EPI  MODELS  IS  DESCRIBED  BELOW. 

MODEL  NAME»<M>  [.FRACKNW]  ( i F AUT  «<N>) 

FRAC  | A MULTIPLIER  OF  The  EVAPORATION  COEFFICIENT  FOR  AUTOOOPING, 

FAUT  | TmE  FRACTION  OF  DOPANT  RELEASED  FROM  THE  INTERFACE  THAT  IS 
REABSORBED. 


SPECIAL  PURPOSE  MODELS 


THE  SPECIAL  PURPOSE  MOOELS  ARE  SO  CALLED  BECAUSE  THEIR  PARAMETERS  ARE 
ONES  THAT  OON'T  FIT  IN  WITH  ANY  OF  THE  OTHER  MODELS,  AND  IN  SOME  CASES  ARE 
NOT  strictly  MODEL  PARAMETERS  IN  THE  SAME  SENSE  AS  THE  OTHERS,  THERE  ARE 
FIVE  SUCH  SPECIAL  PURPOSE  MODELS,  NUMBERED  ONE  THROUGH  FIVE,  THEIR  PARAMETER 
LIST  IS  DESCRIBED  BELOW. 

SPECIAL  PURPOSE  MODELS  I SPMJ,  SPM?,  SPM3,  SPM4,  SPM5 

MODEL  NAME*<M>  t , N I EF »«N > J (,nIEA»«N>)  [,GATEa«G>]  (,QSSQ»<N>) 

* I , CBLK  »<N>) 

NIEF  | THE  PRE-EXPONENTIAL  CONSTANT  USED  TO  CALCULATE  THE  EFFECTIVE 
INTRINSIC  CARRIER  CONCENTRATION,  NI,  AS  A FUNCTION  OF  TEMPER- 
ATURE. (I.E,  NI  ■ NIEF  * EXPf-N IFA/KT)  • (T  **  1.5)  WHERE  K IS 
BOLTIMAN'S  CONSTANT  and  T IS  IN  KELVIN  DEGREES).  DEFAULT 
VALUE  ■ 2.09717E16  ATOMS/CUBTC  CENTIMETER. 


N IEA  l The  ACTIVATION  ENERGY  USED  IN  THE  ABOVE  EXPRESSION  FOR  CALCU- 
LATING THE  EFFECTIVE  INTRINSIC  CARRIER  CONCENTRATION.  DEFAULT 
VALUE  ■ 2,561839  EV. 

GATE  I THE  GATE  MATERIAL  USEO  TO  CALCULATE  THE  THRESHOLD  VOLTAGE, 

there  arf  three  types  that  may  be  specified,  aluminum,  p*,  and 

N*.  Their  SYMBOLS  («G>)  ARE  RESPECTIVELY,  AL,  ♦,  AND 

QSSQ  I THE  NUMBFR  OF  SURFACE  STATES  AT  THF  S I L I CON/OX I DE  INTERFACE  IN 
STATES/SOUARE  CENTIMETER  USED  FOR  CALCULATING  THE  THRESMOLO 

voltages,  default  value  • P , 

CBL*  I INDICATES  THE  TYPE  SUL*  CONTACT. 

C8L*  ■ a.  FOR  CONTACT  AT  THE  SUBSTRATE  (LAST  SI  POINT) 

CBLK  » 1.  FOR  CONTACT  AT  THE  SURFACE  (FIRST  SI  POINT) 

CBLK  • «N>  ARBITRARY  CONCENTRATION  AT  BULK  CONTACT. 

POSITIVE  FOR  P-TYPF , NEGATIVE  FOR  N-TVPE. 


• nnm*niiHr"“*i— 


DEFAULT  parameter  values 


BELOW  ARE  THE  DEFAULT  VALUES  OF  THE  VARIOUS  MODEL  PARAMETERS  USED  IN 

suprem,  those  parameters  that  are  orientation  dependent  and  whose  <100* 
value  is  NOT  known  USE  the  «)11>  value  as  the  default. 


NAME 

ORNT,  I 

BORON 

I 

PHOSPHORUS 

I 

ANTIMONY 

I 

ARSENIC 

I 

SEGREGATION  COEFFICIENT! 

SEGA 

<U1»  I 

1126. A 

I 

10.0 

I 

10.0 

I 

10.0 

I 

SEGA 

<1A0>  I 

2208.0 

I 

1 A , 0 

I 

10.0 

I 

10.0 

I 

SFGE 

<111*  I 

A.  91 

I 

0.0 

I 

0.0 

I 

0.0 

I 

SEGE 

<100*  I 

0.96 

I 

0.0 

I 

0.0 

I 

0.0 

I 

surface 

TRANSPORT  COEFFICIENT! 

STCA 

I 

1 .674E7 

I 

9.0E5 

I 

1.5E3 

I 

9. AES 

I 

STCE 

I 

2.4B1 

I 

1.99 

I 

1 . 0 A 

I 

1,99 

I 

MOV  INI, 

FOUNDRY 

FLUX  TERM! 

hSFA 

I 

A.l 

I 

0.1 

I 

0.1 

I 

0.1 

I 

HSFE 

I 

A.  A 

I 

0.0 

I 

0.0 

I 

0.A 

I 

SILICON 

dioxide 

DIFFUSION 

COEFFICIENT! 

DOXA 

I 

1 .896E6 

I 

4.56E7 

I 

7.66E25 

I 

1 .05E10 

I 

EOXA 

I 

3.33 

I 

3.5 

I 

8.75 

I 

4,09 

I 

SILICON 

DIFFUSION  COEFFICIENT 

f NEUTRAL  OR 

NON 

I-OXIDIZING  AMBIENT)! 

DSTN 

<IU>  I 

3.3324E9 

I 

2 , 3 1 E 1 0 

I 

7.74E1A 

I 

1 .44AE1 1 

I 

OSXN 

< 1 AA»  I 

3.332AE9 

I 

2 . 3 1 E 1 A 

I 

7.74E10 

I 

1 .440E1  1 

I 

ESXN 

< 1 1 1 » I 

3.42565 

I 

3,66 

I 

3.96 

I 

4.00 

I 

ESXN 

<1AA*  I 

3.42565 

I 

3.66 

I 

3.98 

I 

4.00 

I 

SILICON 

DIFFUSION  COEFFICIENT 

(DRY  02  OXIDIZING  AMBIENT)! 

OSXO 

<m»  i 

1 . 6 A8E 1 

I 

2 , 3 I E 1 A 

I 

7.74E10 

I 

1.440E11 

I 

OSXO 

<1AA»  1 

5.5A2E2 

I 

A.ISSEIA 

I 

7.74E10 

I 

1 .440E1  1 

I 

tsxo 

<ui*  i 

1 .453A 

I 

3,66 

I 

3.90 

I 

4,00 

I 

£5X0 

<)A0»  I 

1 , 69A65 

I 

3.66 

I 

3.98 

I 

4,00 

I 

SILICON 

DIFFUSION  COEFFICIENT 

(WET  02  OXIDIZING  AMBIENT)! 

D5XW 

<111*  I 

1 , S A8E 1 

I 

2 , 3 1 E l A 

I 

7.74E10 

I 

1 .440E1 1 

I 

0SX4 

< 1 00*  I 

5.SA2E2 

I 

4.15SE10 

I 

7.74E10 

I 

1 .440E1 1 

I 

FSXW 

<111*  I 

1 .ASA 

I 

3.66 

I 

3. 98 

I 

A, 08 

I 

ESXw 

< 1 AA»  I 

1.690 

I 

3.66 

I 

3.96 

T 

4,00 

I 

8FTAI 

RFTA 

I 

3.0 

I 

1.0 

I 

1.0 

I 

100.0 

X 

CLUSTERING  SI7EI 

CL  ST 

I 

0 t A 

I 

0.0 

I 

0.0 

I 

4.A 

I 

EQUILIBRIUM  CLUSTERING  COEFFICIENT! 

k E Q A 

I 

A,  A 

I 

0.0 

I 

0.0 

I 

2.72E-17 

I 

kEQE 

I 

A.0 

I 

A.  A 

I 

0.A 

I 

•0.305 

I 

NAME  OPNT. 

J BORON 

I PHOSPHORUS 

I 

ANTIMONY 

I 

ARSENIC 

OECUISTERING 

KOC0 

rate: 

I 0.0 

I 0.0 

! 

0.0 

I 

4.BE8 

KDCE 

I 0.0 

! 0.0 

I 

0.0 

I 

2.« 

NAME 

ORNT 

STEAM 

I 

WET 

I 

DRY 

I 

NEUTRAL 

LINEAR 

OXIOE 

growth  rate: 

lrte 

<1U> 

I 2.717E6 

I 

2.717E6 

I 

l .038E5 

I 

0.0 

LRTE 

< 1 00> 

I 1.617E6 

I 

1 .617E6 

I 

6.18114 

I 

0.0 

L RE  A 

I 2.85 

I 

2,05 

I 

2.00 

I 

0.0 

LBE  A 

<i00> 

I 2.05 

I 

2,05 

I 

2.00 

I 

0,0 

PARABOLIC  OXIDE  GROWTH  RATE: 

PRTE 

«t  i i» 

I «.43 

I 

8.43 

I 

12,87 

I 

0.0 

PPTE 

<i00> 

I 6,43 

I 

6.43 

I 

12.87 

I 

0.0 

PREA 

<U1» 

I 0 . 7 fl 

I 

0.78 

I 

1.23 

I 

0,0 

PREA 

<i«0> 

I 0 . 7 B 

I 

0.78 

I 

1.23 

I 

0.0 

r 1 


1 1 — — ■ | | — — I 


REFERENCES  FOR  DEFAULT  PARAMFTER  VALUES 


SEGREGATION  COEFFICIENTS* 

BORON* 

O.  A.  ANTONIADIS.  A.  G.  GONZALEZ  AND  R.  W.  DUTTON,  J.  ELECTROCHEM, 

SOC..  (125),  P.  813,  (1978). 

ALL  OTHER  IMPURITIES  j 

h.  F,  wolf,  SEMICONDUCTORS,  InTERSCIENCE,  1971,  P,  361. 

SURFACE  TRANSPORT  COEFFICIENTS* 

P.  H.  LANGER  ANO  J.  T.  GOLDSTEIN,  J,  ELECTROCHEM  SOC.,  APR.  1974. 

MOVING  BOUNORY  FLU*  TERM* 

SEE  PAGE  21  OF  THIS  TECHNICAL  REPORT. 

SILICON  DIOXIDE  DIFFUSION  COEFFICIENT* 

M.  GHEZZO  AND  0.  M.  BROWN,  J.  ELECTROCHEM  SOC,,  JAN.  1973, 

SILICON  DIFFUSION  COEFFICIENT  (NFUTRAL  OR  NON. OX  I DI Z I NG  AMBIENT)* 

BORON* 

0.  A.  ANTONIADIS.  A.  G.  GONZALEZ  AND  R.  W.  DUTTON,  J.  ELECTROCHEM  SOC., 
(125),  P.  813,  (1978). 

PHOSPHORUS  I 

R.  a.  FAIR  AND  J.  C.  C.  TSAI,  J.  ELECTROCHEM  SOC.,  (124),  P.  1107, 
(1977). 

ANTIMONY  * 

H.  F.  WOLF,  SEMICONDUCTORS,  InTERSCIENCE,  1971,  PR.  153  * 363, 

ARSENIC* 

T.  L.  CHIU  AND  H.  N.  GOSH,  IBM  J.  RES.  PEV , , 15,  472,  (1971). 

SILICON  DIFFUSION  COEFFICIENT  (DRY  OXIDIZING  AMBIENT)* 

BORON  I 

SAME  REFERENCE  AS  FOR  NEUTRAL  AMBIENT  ABOVE. 

PHOSPHORUS* 

R.  B.  FAIR,  PRIVATE  COMMUNICATION. 

ANTIMONY  * 

OIFFUSIVITY  SAME  AS  for  neutral  ambient. 

ARSENIC! 

DIFFUSIVITY  SAME  AS  FOR  NEUTRAL  AMBIENT, 

SILICON  DIFFUSION  COEFFICIENT  (wET  02  OXIDIZING  AMBIENT)* 

OIFFUSIVITY  SAME  AS  FQR  DRY  82  AMBIENT  FOR  ALL  ELEMENTS. 

BETA* 

SFE  PAGE  13  OF  THIS  TECHNICAL  REPORT. 


ALL  ARSENIC  CLUSTERING  INFORMATION 

F.  F,  MQRE  HE  AO , RRIVATE  COMMUNICATION. 


LINEAR  OXIDE  GROWTH  RATE  COEFFICIENT* 

B , E.  OtAL,  J.  ELtCTROCHEM  SOC.,  (125),  P.  57B,  (1976). 

PARABOLIC  OXIDE  GROWTH  RATE  COEFFICIENT* 

SAME  REFERENCE  AS  FOR  LINEAR  GROWTH  RATE  ABOVE. 
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APPENDIX  2 


Following  are  two  examples  of  processes  simulated  by  SUPREM  II. 
The  first  example  is  a CMOS  process  where  the  p-well  fabrication 
is  simulated  and  the  final  n-channel  threshold  voltage  is  calcu- 
lated, first  without  and  then  with  a threshold  voltage  tailoring 
phosphorus  implant.  The  second  example  is  a bipolar  process  including 
an  arsenic  buried  collector,  a boron  base  and  a phosphorus  emitter. 
This  process  is  a variant  of  a Stanford  University  standard  process. 

Both  of  the  following  two  examples  were  designed  so  that  the  use 
of  the  SUPREM  II  simulator  could  be  illustrated  under  a wide  variety 
of  processes.  It  is  not  implied  that  any  of  these  processes  would 
result  into  useful  electronic  devices. 
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•••  STANFORD  UNIVERSITY  PROCESS  ENGINEERING  MODELS  PROGRAM  bbb 


VERSION  0-02  ••• 


1.. ..TITL 

2. . ..GRID 

3. . ..5.BS 


CMOS  P-WE  .L  SIMULATION 
OYSIbPI.01,  OPTHB0 , 6 , YM**b2,5 
ORNT  » l 00 , ELFMb-,  CONC»lElS 


4. . ..COMM 

5. . ..5TEP 


STARTING  OXIDE  THICKNESS  OF  320  A. 
T YPE  »DEPO , TIMEbI,  GRTEbp.052 


6. ...PLOT  T0TL«Y.  CM INb l 4 , NDECbA,  WIND. 1.8 


2. . ..COMM 

8. . ..5TEP 


P-wELL  implant 

TypEbIMPl,  ELEMbP,  DOSE ■ 5E 1 2 , AKEVb200 


9. ...COMM 

12. . ..pLOT 

11. . ..PRINT 


---  STOP  PLOTTING,  start  PRINTING  INFORMATION 

TOTL«N 

ME  AD ■ V 


12. . ..COMM 

13. . ..5TEP 


DIVE-IN  IN  N2  FOR  J.5  HOURS 

TYPEbOXIO,  TEMPfliep,  TIMEbBP,  MODL»NIT0 


14. ...STEP  TYPE«ETCm,  TEmPb25 


15. .  . .COMM 

1 6 .. ..  GR 10 


---EXTEND  GRID  SPACE--- 

OYSI «2 .0 1 5 » OPTHBl.d,  v M A V b 7 ,0 


IT. ...COMM 

18. . ..MODEL 

19. . ..5TEP 


OR  I vE  - IN  for  15  HOURS  In  \9i  X DRr  o? 
NAME«ORvi,  PRESb.I 

TYPE*OXln,  TEmPb 1 1 P0 , TIMEbBPP,  M0DL«0PY1 


22  ...  . COMM 

21. ...STEP 


FIFlD  Ox  ID  GROhTh  in  hFT  0?  FOB  ShOuRS 
TYPEbOXIO,  TEMPBH25,  TIMEBJ00,  MOOLRMET0 


22. ...STEP  TYPEbETCm,  TFmPb25 


23. . ..COMM 

24. . ..5TEP 

23. . ..5TEP 

28. . ..PLOT 
2T... .STEP 


GATE  OXIDATION  AT  K»0f>  c 

TYPEbOXIO,  TFMPBl  upp,  TIMEbS,  MOOL  “DRY0 
TYPEbOXIO,  TFmPbIOPO,  T I ME  »5 , mqOlbmETP 
TOTL«Y,  w I NO • 6 , 

TYPEbOXIO,  TFMPb10PP,  TIME«5,  MOnL^DRYP 


28. ...SAVE  FlLEiPMl ,TYPE»B 


29. . ..COMM 

30. . ..MODEL 

31. . ..5.EP 


ANNEAL  CALCULATE  THRESHOLD  voltage--- 

NAME*SPm  l , GATEbAL,  0SSQ»4E 10 , CBLKsl 

TYPEbOXIO,  TEMPb10P0,  TIMEb30,  MODL»NIT0,  MQOLbSPMI 


32. .  . .COMM 

33. .  . .MODEL 

34. . ..5.EP 


THRESHOLD  TAILORING  IMPLANT 
NAMEbSPMI,  CBLKB6E15 

TYPE«IMPL,  ELEMbP,  00SEb5E11,  AKEVb90,  MODLbSPMI 


35, . , ,EnD 


1 •„  * 


CMOS  P-WElL  SIMULATION 
P-NELl  IMPLANT 
STEP  » 2 TIME  ■ 


P MINUTES 


CMOS  P-WELL  SIMULATION 
DIVE-IN  IN  N2  FOB  1,9  HOUBS 
STEP  R 3 


NEUTRAL  AMBIENT  OBIVE-IN 

TOTAL  step  TIME  • 90. 0 MINUTES 

INITIAL  TEMPERATURE  • 1100. 0P  DEGREES  C, 

OXIDE  THICKNESS  ■ 5.0000E.02  MICRONS 

1 OXIDE  1 SILICON  J l SURFACE  I 

1 DIFFUSION  1 DIFFUSION  1 SEGREGATION  1 TRANSPORT 

l COEFFICIENT  1 COEFFICIENT  i COEFFICIENT  1 COEFFICIENT  1 


BORON  1 2.09716E-07  1 8.90324E.04  j ,60145  1 1.3U32E-02  1 

SURFACE  CONCENTRATION  ■ 4,3 1 0697E* 1 6 AT0MS/CMT3 

Junction  oepth  i sheet  resistance 

......................  i. ........ 

1,60157  MICRONS  1 3722,72  OHMS/SOUARE 

1 66619,9  UHMS/SOUARE 


NET  ACTIVE  CONCENTRATION 


OXIOE 

CHARGE 

• 

3.114949E«10 

IS 

.636 

1 

OF 

total 

SILICON 

CHARGE 

• 

4.8S4020EM2 

IS 

99.4 

X 

OF 

total 

total 

CHARGE 

■ 

4.8SS170EM2 

IS 

90.3 

s 

OF 

INITIAL 

INITIAL 

Charge 

8 

5.071 P22E*12 

CHEMICAL  CONCENTRATION  OF  BORON 


OXIDE 

CHARGE 

B 

3 . 1 1 4050EMP 

TS 

,6?3 

X 

of 

total 

SILICON 

CHARGE 

a 

4 , 966288C* 1 2 

IS 

99,4 

X 

OF 

TOTAL 

total 

CHARGE 

a 

4,997437E*12 

IS 

1 00. 

X 

OF 

INITIAL 

INITIAL 

charge 

a 

4 , W97  420E*  1 2 

i 

I 
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CMOS  P-WELL  SIMULATION 
OIVE-TN  IN  n2  POP  t.S  HOURS 
STEP  « 4 


ETCH  STEP 

ETCh  TEMPERATURE  ■ 25.0  DECREES  C. 

ETCH  TIME  » .0  MINUTES 

ETCH  RATE  ■ .0000  MICBONS/MINUTE 

OXIDE  THICKNESS  ■ 0.0000E*00  MICRONS 

SURFACE  CONCENTRATION  ■ 4.3106O7E*18  AT0MS/CMT3 

JUNCTION  DEPTH  J SHEET  RESISTANCE 



1.49137  MICRONS  1 3722.72  OHMS/SQUARE 

1 66919,0  QHMS/SQUARC 


NET  ACTIVE  CONCENTRATION 

OXIOE  CHARGE  ■ 0.000fl00E*00  IS  0.000E*00  X OP  TOTAL 

SILICON  CHARGE  ■ 4 . 834P20E ♦ 1 ? IS  100.0  X OP  TOTAL 

TOTAL  CHARGE  ■ 4.8S4020E*12  IS  99.4  X OP  INITIAL 

INITIAL  CHARGE  ■ 4,9951 799*12 


CHEMICAL  CONCENTRATION  OF  BORON 

OXIOE  CHARGE  ■ 0.0000006*00  IS  0,0P0E*P0  X OP  TOTAL 

SILICON  CHARGE  • 4 ,966J80E*12  IS  1*0.0  X OP  TOTAL 

TOTAL  CHARGE  ■ 4 . 0A6008E* 1 2 IS  99,4  X OP  INITIAL 

INITIAL  CHARGE  ■ 4.997437E*12 


CMOS  P-wELL  SIMULATION 

DRIvE-IN  FOR  15  HOURS  IN  10X  DRY  02 

STEP  « 5 


OXIDATION  IN  DRY  OXYGEN 

TOTAL  STEP  TIME  • 902.0  MINUTES 

INITIAL  TEMPERATURE  • 1100.00  DEGREES  C. 

OXIDE  THICKNESS  ■ .1404  MICROS 

LINEAR  OXIOE  GROWTH  rate  ■ 3.0074l0E-»4  M I CRONS / M I NUTE 

parabolic  oxide  growth  rate  ■ s.qsbspbe-bs  micronst2/minute 

OXIDE  GRQMTH  PRESSURE  ■ 9.999999E-02  ATMOSPHERES 


1 OXIOE  l 

1 DIFFUSION  1 
1 COEFFICIENT  1 


SILICON 

DIFFUSION 

COEFFICIENT 


SEGREGATION 

COEFFICIENT 


l SURFACE  l 

1 TRANSPORT 
1 COEFFICIENT  1 


BORON  1 2.P9716E-07  1 1,235226-03  1 .66145 

SURFACE  CONCENTRATION  ■ 1.9l527lE*16  ATOMS/CMtS 

junction  depth  i sheet  resistance 




1.31132E-B2  1 


4.P3260  MICRONS  1 4067,52 

1 19710.6 


NET  ACTIVE  CONCFNTRATION 


OMMS/SOUARfc 

OHMS/SOUARE 


OXIOE 

CHARGE 

■ 

5.k)P3394t*l  1 

IS 

10.9 

X 

OF 

total 

SILICON 

CHARGE 

■ 

4, 081520 t*12 

IS 

89.1 

X 

OF 

total 

TOTAL 

CHARGE 

■ 

4.501 B59E*  1 2 

IS 

94.4 

X 

OF 

INITIAL 

INITIAL 

charge 

■ 

4. 85412 4E *12 

CHEMICAL 

CONCENTRATION  of  boron 

OXIDE 

CHARGE 

S 

5.003394E*1 1 

IS 

10.4 

X 

OF 

total 

SILICON 

CHARGE 

a 

4.31 1637E*12 

is 

49.6 

X 

OF 

total 

TOTAL 

ChaRGF 

a 

4.8119766*12 

IS 

96.9 

X 

OF 

INITIAL 

INITIAL  CHARGE  ■ 4.9663916*1? 
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CMOS  R-WELL  SIMULATION 

field  n*IO  GROWTH  IN  wET  02  FOR  5H0URS 

STEP  * 6 


OXIDATION  IN  WET  OXYGEN 

TOTAL  STEP  TIME  ■ 300,0  MINUTES 

INITIAL  TEMPERATURE  ■ 1025. 02  DEGREES  C. 

OXIDE  THICKNESS  ■ 1.203  MICRONS 

LINEAR  OXIDE  GROWTH  RATE  ■ 2.11S309E-02  MICRONS/MINUTE 

PARABOLIC  OXIDE  GROWTH  RATE  ■ 5.71853OE-03  MICR0NST2/MINUTE 

OX IOE  GROWTH  PRESSURE  » . 8 A2 1 05  ATMOSPHERES 


1 OXIDE 
1 DIFFUSION 
1 COEFFICIENT 


1 SILICON 
1 OIFFUSION 
l COEFFICIENT 


1 t 

1 segregation  1 

1 COEFFICIENT  1 


BORON  1 3.74197E-08  1 S.18287E-04  i .41393  1 

SURFACE  CONCENTRATION  ■ 3.349655F*15  ATOMS/CHtS 


JUNCTION  DEPTH  1 

_______  i 

3.6*438  MICRONS  l 

1 


SHEET  RESISTANCE 


5788.51  OHMS/SOUARE 

21049.1  OHMS/SOUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE 

CHARGE 

■ 

1 .7?5380E*12 

IS 

37,8 

I 

OF 

TOTAL 

SILICON 

CHARGE 

■ 

2.859700EM2 

IS 

62.4 

X 

OF 

total 

total 

CHARGE 

■ 

4, 58507 9E *12 

IS 

100. 

X 

OF 

INITIAL 

INITIAL 

charge 

• 

4.5B1S59E*!? 

CHEMICAL 

CONCENTRATION  OF  80R0N 

OXIDE 

CHARGE 

a 

1 ,7253B0E*I2 

IS 

36.0 

X 

OF 

total 

SILICON 

CHARGF 

a 

3.071S98EM2 

IS 

64,0 

X 

OF 

total 

total 

CHARGE 

a 

4 , 798978E*  1 2 

IS 

99.7 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

a 

4,81  1 976E*  1 2 
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SURFACE  1 

TRANSPORT 
COEFFICIENT  1 


3.90494E-03  1 


CMOS  P-wELL  SIMULATION 

FIELD  OX  ID  GROWTH  IN  WET  02  FOR  SHOURS 
STEP  * 7 


t 


! 


I 

I 

! 


: 


ETCH  STEP 

ETCH  TEMPERATURE  ■ 
ETCH  TIME  a 
ETCH  RATE  ■ 
OXIOE  THICKNESS 


25. P DEGREES  C. 

.0  MINUTES 

.8002  MICRONS/MINUTE 
a 0 . 0000E*80  MICRONS 


SURFACE  CONCENTRATION  a 3.349635E*15  AT0MS/CMt3 


JUNCTION  DEPTH  1 

_ _____________  i 

3.68437  MICRONS  I 

1 


SHEET  RESISTANCE 


5788,55  OHMS/SQUARE 

21840.1  OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIOE 

CHARGE 

• 

0,e00000E»00 

IS 

0.800E»00 

X 

OF 

total 

SILICON 

CHARGE 

■ 

2 . 859700E*  1 2 

IS 

180.0 

X 

OF 

TOTAL 

TOTAL 

CHARGE 

■ 

2 ,859700E*12 

IS 

62.4 

X 

OF 

INITIAL 

I ITIAL 

CHARGE 

B 

4 , 585079E*  1 2 

C EMICAL  CONCENTRATION  OF  BORON 


1 

OXIOE 

CHARGE 

• 

0 , 000000E *00 

IS 

0 ,000E*00  X 

OF 

total 

* 

SILICON 

CHARGE 

8 

3.07l598E*12 

IS 

100,0  % 

OF 

total 

TOTAL 

charge 

8 

3.87lS98E»12 

IS 

64.0  X 

OF 

INITIAL 

♦ 

INITIAL 

charge 

a 

4.7*6978E«1» 

I 


I , 
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CMOS  P-WELL  SIMULATION 
GATE  OXIDATION  AT  1800  C 
STEP  N 8 


OXIDATION  IN  DP V OXYGEN 

TOTAL  STEP  TIME  » 5.0  MINUTES 

INITIAL  TEMPERATURE  ■ 1000,000  DEGREES  C. 

OXIDE  THICKNESS  ■ 2.0675E-02  MICRONS 

LINEAR  OXIDE  GROWTH  RATE  ■ 8.019000E-03  MI CRONS /MINUTE 
PARABOLIC  OXIOE  GROWTH  RATE  ■ 1.705008E-04  M ICRONS t2/M I NUTE 
OXIDE  GROWTH  PRESSURE  ■ t. 00000  ATMOSPHERES 

I OXIOE  t SILICON  I 1 SURFACE  J 

1 DIFFUSION  I DIFFUSION  J SEGREGATION  J TRANSPORT 

1 COEFFICIENT  { COEFFICIENT  J COEFFICIENT  1 COEFFICIENT  l 


BORON  1 2 , 0 1 368E-00  i 2.03949E-04  I .34974  1 2.52S26E-03  t 

SURFACE  CONCENTRATION  » 1 , 7 24042E* 1 5 ATOMS/CMt3 

JUNCTION  DEPTH  1 SHEET  RESISTANCE 

— J 

3 , 67869  MICRONS  1 5B21.S1  OHMS/SQUARE 

I 2 1 0 B 1 . 9 OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIOE 

charge 

■ 

1 ,1  1 3799E* 1 0 

IS 

.390 

X 

OF 

TOTAL 

SILICON 

charge 

• 

2.8441B4EM2 

IS 

99.0 

X 

OF 

total 

total 

CHARGE 

■ 

2.855322E*12 

IS 

99.8 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

• 

2 , 899700E*  1 2 

CHEMICAL  CONCENTRATION  OF  BORON 


OXIOE 

charge 

• 

1 .1  1 3799E*  1 0 

IS 

,383 

X 

OF 

TOTAL 

SILICON 

charge 

• 

3.055481EM2 

IS 

09.0 

X 

OF 

total 

TOTAL 

Charge 

• 

3.0«6Sl9E*i? 

IS 

99.8 

X 

OF 

INITIAL 

INITIAL 

charge 

■ 

3.07lS98E*12 

BORON 


2.013681-08  1 2.03949E-04  | 


.34974 


2.52620E-B3  1 


SURFACE  CONCENTRATION  • 6.7580766*14  ATOMS/CMfS 


JUNCTION  DEPTH 
3.B5247  MIC 


SHEET  RESISTANCE 


MICRONS  l 3872. 50 
1 21077.7 


OHMS/SQUARE 

OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIOE 

CHARGE 

■ 

3.4229246*10 

IS 

1.20 

X 

OF 

total 

SILICON 

CHARGE 

• 

2.8213651*12 

IS 

98.8 

X 

OF 

total 

total 

CHARGE 

s 

2 . 835594E*  1 2 

IS 

100, 

X 

OF 

INITIAL 

initial 

CHARGE 

s 

2.8533226*12 

chemical 

, CONCENTRATION  OF  BORON 

OXIOE 

CHARGE 

■ 

3.422924E*10 

IS 

1.12 

X 

OF 

TOTAL 

SILICON 

charge 

■ 

3.0304036*12 

IS 

98.9 

X 

OF 

total 

total 

CHARGE 

■ 

3.0648321*12 

IS 

99.9 

X 

OF 

INITIAL 

initial 

charge 

• 

3.0668196*12 

CMOS  P-WELL  SIMULATION 

gate  oxidation  at  tone  c 

STEP  II  10 


OXIDATION  IN  DRV  OXVGEN 

total  STEP  time  • 5,0  MINUTES 

INITIAL  TEMPERATURE  ■ 1000,0100  DEGREES  C, 

OXIDE  THICKNESS  ■ S.0046E-02  MICRONS 


LINEAR  OXlOE  GROWTH  RATE  ■ 
PARABOLIC  OXIDE  GROWTH  RATE  ■ 
OXIDP  GROWTH  PRESSURE  ■ 


8.019060E-04  MICRONS/MINUTE 
1.765125E.04  MICRONStS/MINUTE 
1,00000  ATMOSPHERES 


1 OXIDE 
1 OIFPUSION 
1 COEFFICIENT 


1 SILICON 
1 DIFFUSION 
1 COEFFICIENT 


1 1 
l SEGREGATION  1 
1 COEFFICIENT  1 


BORON  1 2.0136RE-0R  2 2.03949E-04  J ,34974  1 

SURFACE  CONCENTRATION  « 3,119767E*15  AT0MS/CMf3 

Junction  depth  i sheet  resistance 

... i 

3,65270  MICRONS  1 5882,82  OHMS/SOUARE 

1 21067,8  OHMS/SOUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE 

CHARGE 

■ 

3.929135E*10 

IS 

1.36 

X 

OF 

total 

SILICON 

CHARGE 

S 

2.815937E*12 

IS 

98.6 

X 

OF 

total 

total 

CHARGE 

9 

2 , 8 55228E ♦ 1 2 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

9 

2 .855594E*!  2 

CHEMICAL  CONCENTRATION  OF  80R0N 


OXlOE 

charge 

■ 

3.92913SE*10 

IS 

1,26 

X 

OF 

TOTAL 

SILICON 

CHARGE 

■ 

3 .0251 67E*1 2 

IS 

98.7 

X 

OF 

total 

TOTAL 

CHARGE 

t 

3.064A59EM2 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

3.064632E*12 

SURFACE  1 

TRANSPORT 
COEFFICIENT  1 


2.52626E-03  I 


I 


CMOS  P-wELl  SIMULATION 
GATE  OXIDATION  AT  1*00  C 
STEP  • 10  TIME  » 5.0  MINUTES. 


i 


L 1 


DEPTH  1 
(UM)  1 


• .08 


2.00 


A. 00 


8.00 


0,00 


CONCENTRATION  (LOG  ATOMS/CC) 


4 

15 

16 

17 

16 
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CMOS  P-WELL  SIMULATION 

ANNEAL  — CALCULATE  THRESHOLD  VOLTAGE  — « 
STEP  P U 


NEUTRAL  AMBIENT  DPIVE-IN 

total  STEP  time  ■ 3(1 ,0  MINUTES 

INITIAL  TEMPERATURE  « 1000,000  DEGREES  C. 

OXIDE  THICKNESS  • S.0046E-02  MICRONS 


1 OXIDE 
l DIFFUSION 
l COEFFICIENT 


1 SILICON 
1 DIFFUSION 
1 COEFFICIENT 


t 

1 SEGREGATION 
l COEFFICIENT 


SURFACE  I 

TRANSPORT 
COEFFICIENT  1 


BORON  1 2 , 0 1 368E.0B  I 9.16197E.05  1 ,34074 


2.32626E-03  1 


SURFACE  CONCENTRATION  ■ 3,636077E*13  ATOMS/CMt3 


GATE  MATERIAL  ■ 
OXIDE  THICKNESS  • 
THRESHOLD  VOLTAGE  • 


ALUMINUM  SILICON  UNDER  GATE  ■ P - TVPE 
809,5  ANG.  CAPACITANCE/AREA  ■ 4.26E-04  PF/UM2 
.58  VOLTS  AT  SURFACE  STATES  ■ 4.00EM0 


JUNCTION  OEPTH  1 SHEET  RESISTANCE 

....j.._.  

3,65608  MICRONS  l 3887,67  OHMS/SQUARE 

1 21109,7  OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE 

CHARGE 

a 

4 .229726EM0 

IS 

1.46 

X 

OF 

total 

SILICON 

CHARGE 

a 

2.812097E*12 

IS 

08.5 

X 

OF 

total 

TOTAL 

CHARGE 

a 

2.854394E*12 

IS 

100,0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

a 

2 . 85522BE* 1 2 

CHEMICAL  CONCENTRATION  DF  BORON 


OXIDE 

CHARGE 

8 

4,229726E*10 

IS 

1.38 

X 

OF 

total 

SILICON 

CHARGE 

8 

3.022125EM2 

IS 

98.6 

X 

OF 

total 

total 

CHARGE 

8 

3 , 06 4 4 22E+  J 2 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

8 

3 . 064459E* 1 2 
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CMOS  P-WELL  SIMULATION 


ANNEAL 
1 STEP 
l 

OEPTH  1 

CUM ) 1 

14 

...  CALCULATE 
• 11  TIME 

15 

THRESHOLO  VOLTAGE  — 

■ 30.0  MINUTES. 

CONCENTRATION  (LOG  ATOMS/CC) 

16 

IT 

18 
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STANFORD  UNIV  CALIF  STANFORD  ELECTRONICS  LABS  F/G  9/2 

SUPREM  II  --  A PROGRAM  FOR  IC  PROCESS  MODELING  AND  SIMULATION. (U) 
JUN  78  D A ANTONIADIS.  S E HANSEN.  R W DUTTON  DAAG-29-77-C-0006 
SU-SEL-78-020  ARO-14589.2-EL  NL 


CMOS  P.MELL  SIMULATION 
THRESHOLD  TAILORING  IMPLANT 
STEP  m 12 


TON  IMPLANT  (GAUSSIAN  APPROXIMATION) 
IMPLANTEO  IMPURITY  ■ PHOSPHORUS 

implanted  dose  ■ s.ppR0Pee*u 

IMPLANT  ENERGY  ■ 90,0000 

BANGF  ■ 9 , 664024E-02 

STANDARD  OEVIATION  ■ 4.174892E-02 

PEAK  CONCENTRATION  ■ 3.023804F*16 


SURFACE  CONCENTRATION  ■ .4 , P60335E* 1 6 ATOMS/CMtS 


GATE  MATERIAL  ■ 
OXIDE  THICKNESS  ■ 
THRESHOLD  VOLTAGE  ■ 

JUNCTION  DEPTH 


ALUMINUM  SILICON  under  GATE  ■ N - TYPE 
809.5  ANG.  CAPACITANCE/AREA  ■ 4.26E-04  PF/UM2 
.15  VOLTS  AT  SURFACE  STATES  • 4,00E*10 

1 SHEET  RESISTANCE 


.101351 

3.65608 


MICRONS  1 
MICRONS  1 
1 


24786.9 

6092.49 

21109.4 


OHMS/SQUARE 

OHMS/SQUARE 

OHMS7SOUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE  CHARGE  4 

l .l*1700E*tl 

IS 

3.79 

X 

OF 

TOTAL 

1 

SILICON  CHARGE  ■ 

3 . 002370E*  1 2 

IS 

96.2 

X 

OF 

TOTAL 

TOTAL  CHARGE  » 

3.1 20349E*  1 2 

IS 

109. 

X 

OF 

INITIAL 

INITIAL  CHARGE  ■ 

2,854394E*12 

CHEMICAL  CONCENTRATION  OF  RORON 

0 

OXIDE  CHARGE  ■ 

4 . 2297  26E*  1 0 

IS 

1.38 

X 

OF 

TOTAL 

SILICON  CHARGE  * 

3.022125EM2 

IS 

98.6 

X 

OF 

total 

TOTAL  CHARGE  » 

3 , 064  422E ♦ 1 2 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL  CHARGE  ■ 

3.084422E*12 

CHEMICAL  CONCENTRATION  OF  PHOSPHORUS 

1 t | 

OXIOE  CHARGE  ■ 

1 ,413504E*U 

IS 

28.3 

X 

OF 

TOTAL 

SILICON  CHARGE  ■ 

3 . 58307  3E*  1 1 

IS 

FI. 7 

X 

OF 

TOTAL 

I 

TOTAL  CHARGE  ■ 

4.996577E*!! 

IS 

0 . 000E  *00 

X 

OF 

INITIAL 

initial  Charge  ■ 0.000000E*00 
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***  STANFORD  UNIVERSITY  PROCESS  ENGINEERING  MODELS  FROGRAM  •». 

version  0-02 


1.. . .TITLE  STANFORO  KITCHIP  PROCESS  8/5/78 

2. . ..GRID  YMAXbS,  OPTHbj,  DYSIb.P1 

3. . ..5.BS  ELEmbB  ,C0NCbiE15  ,ORnTb100 

A. ...PRINT  HEAD* Y 

5. . ..COMM  BURIED  LAYER 

6. . ..5.EP  TYPEbIMPL,  DDSE ■ 1 . 5E 1 5 . AKEVb30,  ELEMbAS 

7. . ..plot  totlby,  windba 

8. . ..5.EP  TYPEbDYIO,  TEMPbISSJ),  TIMEbISB,  mool-oryp 

9. . ..PLOT  TOTL'N 

IP. ...STEP  TYPE«ETCM,  TEMPb25 

11.. ..  COMM  EPITAXY 

12. . ..PLOT  TOTL«Y,  bINObB 

13. . ..5.EP  TYPEbEPIT,  TEMP. 1000,  TIMEbji,  GRTEb.5,  ElFMbaS,  CONCb1E15 

14. . ..PLOT  TOTLbN 

15. . ..COMM  ISOLATION 

18..  ..5.EP  TYPE»OXlO,  TEMPbUPP,  TIME.60,  MODLbDRYP 

17. . ..5.EP  TYPE'OYIO,  TEMPb1200,  TIMEbGP,  MODL"ORY0 

18. . ..plot  totl»t 

19. . ..5.EP  TYPE*ETCM,  TEMP.25 

20. . ..COMM  BASE  DIFFUSION 


21. . ..PLOT 

22. . ..5TEP 

23. . ..PLOT 

24. . ..5TEP 

25 .. . . 5TEP 

26. . ..5TEP 

27. . ..5TEP 
2«. ...PLOT 


M I NO ■ 4 

TYPEbPOEP  ,TEMPb95P  ,TIMEb30  .ELEMbB  ,CONCb1.2E20 
totl»n 

TYPEbOKID  ,TFMPb1000,TImEb45  ,MODLbORY0 
TYPE«0YI0  ,TEMPb10P0,TIMEb80  ,MOOL»WET0 
TYPEbOxIO  , TCMRb 1000, TIME b5  ,MOOL«DRY0 


TYPEbOXIO 
TOTl»Y 


,TEmPb10P0,TIMEb10  .MOOLPNITR 


29, ...STEP 

TYPEbETCM  ,TFMPb25 

‘ 

3P. ...SAVE 

FILEbEBNCF  .TYPEbB 

B 

1 | 

31. ...COMM 

EMITTER  DIFFUSION 

1 i 

32. ...MODEL 

NAMEbMPhP  , STC0B0 

t j 

33. ...STEP 

TYPEbPOEP  ,TFMPb1025 

, T I ME  *33 

, ElEMbP  , COnC ■ 1 , 1 SE2 1 

34. ...PLOT 

TOTlbn 

35. ...STEP 

type«oxio  ,tempbB00 

, T IM£b 10 

.MQDLbNITB  ,MODL«MPH2 

3 STEP 

TYPEbOXIO  ,TEMPb900 

.TIMEbip 

.MODLBWET0 

1 . 

37. ...SAVE 

FILE«EBnCL,  TYPEbB 

' 

38. . . .PLOT 

totlby 

f 

1 

1 

39. ...STEP 

TYPEBOXIO  ,TEMPb900 

, T I ME ■ 1 0 

.MOOLbNITB 

V * 

40. ...END 
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STANFORD  KITCHIP  PROCESS 
BURIED  LAYER 
STEP  • 1 


6/3/70 


ION  IMPLANT  (GAUSSIAN  APPROXIMATION) 


: . 


IMPLANTED  IMPURITY 
IMPLANTEO  DOSE 
IMPLANT  ENERGY 
RANGE 

standard  deviation 

PEAK  CONCENTRATION 


ARSENIC 
1 .50P000E413 
9(1 ,000(9 
2.6SB498E-02 
7.8821B5E-03 
7 , 789628E  *20 


SURFACE  CONCENTRATION  • 9.97082BE*14  ATQMS/CMtS 


NET  ACTIVE  CONCENTRATION 


OXIOE 

charge  ■ 

0.000000E*00 

IS 

0 .000E+00 

X 

OF 

TOTAL 

SILICON 

CHARGE  ■ 

1 .449207EM5 

IS 

100.0 

X 

OF 

TOTAL 

TOTAL 

CHARGE  • 

t ,449?07E*15 

IS 

2.89BE*05 

X 

OF 

INITIAL 

INITIAL  CHARGE  ■ S.000000EM1 


CHEMICAL  CONCENTRATION  OF  80R0N 


Jr 

OXIOE 

CHARGE 

■ 

0 ,000000E*00 

IS 

0.000E*00 

X 

OF 

total 

r- 

SILICON 

CHARGE 

« 

5.000000EM1 

IS 

100.0 

X 

OF 

TOTAL 

TOTAL 

CHARGE 

■ 

3.000000EM1 

IS 

100.0 

X 

OF 

INITIAL 

t 

INITIAL 

CHARGE 

a 

5.000000E*11 

i 

CHEMICAL 

, CONCENTRATION  OF  ARSENIC 

i 

OXIOE 

CHARGE 

a 

0.00000PE+00 

IS 

0 ,000E*00 

X 

OF 

TOTAL 

SILICON 

CHARGE 

a 

1.440723EM3 

IS 

100. 

X 

OF 

TOTAL 

• 

total 

CHARGE 

a 

1 .448723EM5 

IS 

0.000E*60 

X 

OF 

INITIAL 

initial 

CHARGE 

a 

0.000000E*00 

-92- 


STANFORD  KJTCMIP  PROCESS 
BURIED  LAVER 
STEP  « 2 


6/5/78 


OXIDATION  IN  DRV  OXVGEN 

TOTAL  STEB  TIME  a 180,0  MINUTES 

INITIAL  TEMPERATURE  ■ 1 250 a 00  DEGREES  C, 

OXIDE  THICKNESS  ■ .4177  MICRONS 

linear  OXIDE  GRONTH  RATE  ■ 1.577485E-02  MICRONS/MINUTE 

RaBaROLZC  OXIDE  GROhTM  RATE  ■ 1.113008E-03  mICR0nSt2/mINUTE 

OXIDE  GRQmTh  PRESSURE  ■ 1.00000  ATMOSPHERES 


l OXIDE 
1 oiefusion 

1 COEEEICIENT 


1 SILICON  1 1 SURFACE  1 

l OIEFUSION  1 SEGREGATION  1 TRANSPORT 

{ COEFFICIENT  j COEFFICIENT  1 COEFFICIENT  1 


BORON  1 3.95869E-06  1 1.68161E-02  1 1,4700 

■■•■•■•■■■•■•■••■■■■■■■■■•■•■■•■■••■•••••••■••••••■••■■I 

ARSENIC  l 6 , 93 1 86E-07  J 4.55201E-P3  ! 10.0000 


1 .10339  I 

immmmmmmmmmmmmmmmrnt 

1 .23421  i 


SURFACE  CONCENTRATION  a -7 . 765 1 7 1 E * 1 8 AT0MS/CMt3 
JUNCTION  DEPTH  1 SHEET  RESISTANCE 


28,4196 


OHMS/SOUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE 

charge 

■ 

8.700552EM13 

IS 

5.92 

X 

OF 

total 

SILICON 

charge 

■ 

1 .382655EM5 

IS 

94.1 

X 

OF 

TOTAL 

total 

Charge 

a 

1 ,469560E*15 

IS 

101  . 

X 

OF 

INITIAL 

initial 

charge 

a 

1 ,449207E*15 

Ch£miC*l 

, CONCENTRATION  OF  BORON 

oxide 

Charge 

a 

2 ,784715E*10 

IS 

5.60 

X 

OF 

total 

silicon 

Charge 

a 

4 , 6 98  3 7 2E • 1 1 

IS 

94.4 

X 

OF 

total 

total 

Charge 

a 

4 . 976843E*  1 1 

IS 

99,5 

X 

OF 

INITIAL 

Initial 

charge 

a 

5.0P0000EM1 

chemical 

, CONCENTRATION  of  ARSENIC 

OXIDE 

Charge 

a 

S.703345E*13 

IS 

5,92 

X 

OF 

total 

SILICON 

charge 

a 

t .383032E*15 

IS 

94.1 

X 

OF 

total 

total 

CHARGE 

a 

1 .470P66EMS 

IS 

101  . 

X 

OF 

INITIAL 

initial 

CHARGE 

a 

1 .448723E*15 

6/3/78 


OEPTM 

(UN) 

•.42 


14 


.08 


1.80 


2.00 


3.00 


STANFORD  KITCHIP  PROCESS 
BURIED  LAVER 

STEP  • 2 TIME  • 180.0  MINUTES. 


4.00 


13 


CONCENTRATION  (LOG  ATOMS/CC) 
16  17  18  19 


20 


21 
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9TANF0R0  KITCHIP  PROCESS  8/3/70 

BURIED  LAVER 
STEP  • 3 


ETCH  STEP 
ETCH  TEMPERATURE 
ETCH  TIME 
ETCH  RATE 
OXIDE  THICKNESS 


29.0  DEGREES  C. 

.0  MINUTES 

.0000  MICRONS/KINUTE 
• 0 , P000E*00  MICRONS 


surface  concentration  ■ -7,765i7ie*i8  atoms/cmts 

junction  depth  i sheet  resistance 

i 

1 28.4198  OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE 

CHARGE 

■ 

0 ,00000PE*00 

IS 

0,000E*00 

1 

OF 

TOTAL 

SILICON 

CHARGE 

■ 

1 .382535E*15 

IS 

100.0 

X 

OF 

total 

total 

CHARGE 

• 

1 .38253SEM3 

IS 

9A.1 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

• 

1 .A6938PEM5 

chemical  concentration  of  boron 


OXIOE 

CHARGE 

■ 

0,000000E*00 

IS 

0.080E»00 

X 

OF 

TOTAL 

SILICON 

charge 

■ 

4.698372EM1 

IS 

100.0 

X 

OF 

total 

TOTAL 

charge 

s 

4 . 698372E*  1 1 

IS 

04.4 

X 

OF 

INITIAL 

initial 

CHARGE 

■ 

4.97  60  4 3EM1 

CHEMICAL  CONCENTRATION  OF  ARSENIC 


OXIOE 

CHARGE 

■ 

0 ,000000E*00 

IS 

0 , 000E*00  X 

OF 

TOTAL 

SILICON 

CHARGE 

■ 

1 .383P30EM3 

IS 

100.0  X 

OF 

TOTAL 

TOTAL 

CHARGE 

■ 

1 .383P32EM5 

IS 

94.1  X 

OF 

INITIAL 

INITIAL 

charge 

■ 

1 .470066EM3 

STANFORD  KITCHIP  PROCESS 

EPITAXY 

STEP  « 4 


8/3/78 


i 

i 

. 

l 

: 

L i 

if. 

u 

11 


EFt-LAYER  IMPURITY  ■ ARSENIC 
TEMPERATURE  ■ 10*0.0  DEGREES  C 
TIME  ■ U.0  MINUTES 

GROWTH  RATE  ■ .300  MICRONS/MINUTE 

SILICON  AOOEO  « 3.5000  MICRONS 
GAS-PHASE  CONCENTRATION  ■ l.000000E*13 


FAUT  • 


BORON 


,000000E*00  FRAC  ■ ,000000E»0P  FINT 


, 0000P0E*00 


l OXIDE  1 SILICON  1 

1 DIFFUSION  1 DIFFUSION  J 

1 COEFFICIENT  l COEFFICIENT  1 


1 SURFACE  1 

SEGREGATION  1 TRANSFORT 
COEFFICIENT  | COEFFICIENT  l 


1 3,958896-06  1 9.16197E-05  1 1.4706  i 2.52828E-03  I 


ARSENIC  I 8.93188E-07  1 1.F1692E-03  1 10.0000  t 1.19298E-02  1 


SURFACE  CONCENTRATION  ■ - 1 . 000000E* 1 3 AT0MS/CMt3 

JUNCTION  OEPTH  l SHEET  RESISTANCE 

I - — 

) 28.2615  OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIOE 

CHARGE 

• 

0.000P00E*00 

IS 

0.000E+00 

X 

OF 

total 

SILICON 

CHARGE 

• 

1 .379372E*15 

IS 

100.M 

X 

OF 

TOTAL 

TOTAL 

CHARGE 

• 

1 .379372EM5 

IS 

99.8 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

« 

1 .382555EM5 

CHEMICAL 

, CONCENTRATION  OF  SORON 

OXIOE 

CHARGE 

• 

0 .0000006*00 

IS 

0.000E«00 

X 

OF 

total 

SILICON 

CHARGE 

• 

4.702829E*11 

IS 

100.0 

X 

OF 

total 

TOTAL 

CHARGE 

■ 

4 , 702629E* 1 1 

IS 

100. 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

4.698372EM1 

CHEMICAL 

, CONCENTRATION  OF  ARSENIC 

OXIDE 

CHARGE 

■ 

0.0000P0E*00 

IS 

0.000ET00 

X 

OF 

total 

SILICON 

CHARGE 

■ 

1 .379850EM3 

IS 

100.0 

X 

OF 

TOTAL 

total 

CHARGE 

• 

1.379B50EMS 

IS 

99.8 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

1 .383032E*  13 
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STANFORD  KITCHIP  PROCESS  6/3/78 

ISOLATION 

STEP  R 3 


OXIDATION  IN  DRY  OXYGEN 

TOTAL  STEP  TINS  ■ 60.0  MINUTES 

INITIAL  temperature  • UOO.eB  decrees  c. 

OXIDE  THICKNESS  ■ .1104  MICRONS 

LINEAR  OXIDE  GROWTH  RATE  ■ 3.007418E-03  MICRONS/MINUTE 

PARABOLIC  OXIDE  GROWTH  RATE  « 3.986828E-04  M I CRONS T2 /M I NUTE 

OXIDE  GROWTH  PRESSURE  ■ 1,00000  ATMOSPHERES 

1 OXIDE  1 SILICON  1 1 SURFACE  1 

1 DIFFUSION  1 DIFFUSION  J SEGREGATION  1 TRANSPORT 

1 COEFFICIENT  1 COEFFICIENT  l COEFFICIENT  \ COEFFICIENT  1 

BORON  1 2.0971SE-07  1 1.23522E-03  1 ,66149  | 1.3U32E-02  1 


ARSENIC  I 1.18402E-06  1 1.52575E-04  j 10.O00O  1 4.47011E-02  t 

SURFACE  CONCENTRATION  ■ - 1 , 383323E* 1 9 ATOMS/CMtS 
JUNCTION  DEPTH  1 SHEET  RESISTANCE 

1 27,4748  OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE 

charge 

■ 

1 .3604396*09 

IS 

9.863E-05 

X 

OF 

total 

SILICON 

CHARGE 

■ 

1 .3793286*13 

IS 

lOO.O 

X 

OF 

total 

total 

CHARGE 

■ 

1 . 3793306*13 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

1 ,379372E*15 

CHEMICAL  CONCENTRATION  OF  BORON 


OXIDE 

CHARGE 

■ 

4.494211E-10 

IS 

9. 4726-20 

X 

OF 

total 

SILICON 

CHARGE 

■ 

4.7026366*11 

IS 

10O.0 

X 

OF 

TOTAL 

total 

CHARGE 

■ 

4.702«36E*11 

IS 

100. 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

4 , 702629E* 1 1 

CHEMICAL 

CONCENTRATION  OF  ARSENIC 

0XI0E 

CHARGE 

■ 

1.3604396*09 

IS 

9.860E-05 

X 

OF 

TOTAL 

SILICON 

CHARGE 

■ 

1 .3798136*13 

IS 

100.0 

X 

OF 

total 

total 

CHARGE 

■ 

1 ,379817E*15 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

1 ,379830E*18 
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Stanford  kitchir  process  6/3/78 

ISOLATION 
STEP  * 6 


OXIDATION  IN  DRY  OXYGEN 

TOTAL  STEP  TINE  ■ 90.0  minutes 

INITIAL  TEMPERATURE  • 1200.00  DEGREES  C. 

OXIDE  THICKNESS  ■ .2674  MICRONS 


LINEAR 

DXIOE 

GROWTH  RATE 

• 

9.425977E-03  MICRONS/MINUTE 

PARABOLIC 

OXIDE 

GROWTH  RATE 

■ 

8.088042E-04  HICRONST2/MINUTE 

OXIDE  GROWTH  PRESSURE 

■ 

1.00000  ATMOSPHERES 

1 

OXIDE 

1 

SILICON  1 

t 

SURFACE 

I 

1 

OIEEUSJON 

1 

DIFFUSION  J 

SEGREGATION  j 

TRANSPORT 

I 

COEFFICIENT 

1 

COEFFICIENT  l 

COEFFICIENT  l 

COEFFICIENT 

1 

BORON 

1 

1 .38B96E-06 

1 

7.263l«E-03  l 

1.1473  1 

3.44300E-02 

1 

ARSENIC 

1 

1 .95738E-07 

1 

1.38489E-03  1 

10.0000  1 

.14000 

1 

SUREACE  CONCENTRATION  ■ -1.1997826*13  AT0MS/CMt3 
JUNCTION  DEPTH  J SHEET  RESISTANCE 


1 


l 24.7982 

OHMS/SQUARE 

NET  ACTIVE  CONCENTRATION 

OXIDE 

CHARGE  * 

3 , 329523E*09 

IS 

2.4I4E-04 

X 

OF 

total 

SILICON 

CHARGE  * 

1 ,379357E*15 

IS 

100.0 

X 

OF 

TOTAL 

TOTAL 

CHARGE  ■ 

1 ,3793606*15 

IS 

100. 

X 

OF 

INITIAL 

INITIAL 

CHARGE  * 

1 , 379330E* 1 5 

CHEMICAL 

CONCENTRATION  OF  BORON 

OXIDE 

CHARGE  ■ 

19582.8 

IS 

4.163E-08 

X 

OF 

total 

SILICON 

CHARGE  ■ 

4.7021716*11 

IS 

100.0 

X 

OF 

total 

total 

CHARGE  ■ 

4.7021726*11 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE  ■ 

4, 7826366*11 

CHEMICAL 

CONCENTRATION  OF  ARSENIC 

OXIDE 

CHARGE  ■ 

3.3295436*09 

IS 

2.413E«04 

X 

OF 

total 

SILICON 

CHARGE  ■ 

1 .379831E*15 

IS 

100.0 

X 

OF 

total 

TOTAL 

CHARGE  ■ 

1.3790346*15 

IS 

100. 

X 

OF 

INITIAL 

INITIAL 

CHARGE  ■ 

1 .3798176*15 

STANFORD  KITCHIP  PROCESS  8/3/78 

ISOLATION 

STEP  * 7 


ETCH  STEP 

ETCH  TEMPERATURE  » 23.0  DEGREES  C. 

ETCH  TIME  ■ .0  MINUTES 

ETCH  RATE  » ,00010  MICR0N3/MINUTS 

OXIDE  THICKNESS  « 0.0000E+00  MICRONS 

SURFACE  CONCENTRATION  ■ -1 ,199782E*15  ATOMS/CMT3 

JUNCTION  DEPTH  l SHEET  RESISTANCE 

1 

I 24,7981  OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE 

CHARGE 

■ 

0 ,000000E*00 

IS 

0,00PE*00 

X 

OF 

total 

SILICON 

Charge 

■ 

1 ,379353E*1S 

IS 

100,0 

X 

OF 

total 

TOTAL 

CHARGE 

■ 

1 ,379353E*15 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

a 

1 .379360E+15 

CHEMICAL  CONCENTR 

AT  ION  OF  BORON 

oxide  CHARGE  ■ 

0 ,000000E«00 

IS 

0.000E*00  X 

OF 

TOTAL 

SILICON  CHARGE  ■ 

4.702186EM  1 

IS 

100.0  X 

OF 

TOTAL 

TOTAL  CHARGE  ■ 

4,7021 88E*  1 1 

IS 

100.  X 

OF 

INITIAL 

INITIAL  CHARGE  * 

4.702172EM1 

CHEMICAL 

CONCEN 

TRAT 

ION  OF  arsenic 

OXIDE 

CHARGE 

a 

0 

, 000000E*00 

IS 

0 , 000E*00 

X 

OF 

total 

SILICON 

CHARGE 

a 

i 

,379827E*15 

IS 

100.0 

X 

OF 

total 

TOTAL 

CHARGE 

a 

i 

.379827EM3 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

a 

i 

,379834E*15 

STANFORD  kJTChxp  PROCESS  8/3/78 

BASE  DIFFUSION 
STEP  R 8 


GASEOUS  PPEDEPPSIUOn 

total  STEP  time  ■ JP.P  minutes 

INITIAL  temperature  ■ 950.088  DEGREES  C. 

OXIDE  Thickness  ■ 0 ,9000E*0fl  microns 

PR£OEpOSI T ION  impurity  ■ BO«ON 

GAS  CONC.  AT  INTERFACE  • 1.2000PPF«20  ATO*S/CC 

1 OXIDE  l SILICON  l l Surface  l 

1 DIFFUSION  l OIFFUSION  1 SEGREGATION  l TRANSPORT 

: coefficient  j coefficient  j coefficient  j coefficient  j 


BORON  I 5.40470E-P9  I 2.SS8S8E-05  J .244S7  1 1.0000 


ARSENIC  : 7.43S46E-11  1 2.22365E-06  J 10.000P  1 5.88350E-03  i 

SURFACE  CONCENTRATION  ■ 1 . 19871 9E *20  ATQMS/CmtS 


Junction  depth  ; 

SHEET 

RESISTANCE 

.388777  -ICRONS  j 

59.3145 

24.8011 

OHMS/SOUARE 

Oh-S/SOI/ARE 

NET  active  concentration 


OXIDE 

C-ARGF 

■ 

3 . 0000P:»E*0P 

IS 

0,000Et00  X 

OF 

TOTAL 

SILICON 

Charge 

t 

3.263K40EMS 

IS 

100.0  X 

OF 

total 

total 

Charge 

a 

3 . 28354  0E ♦ 1 5 

IS 

237.  X 

OF 

initial 

initial 

charge 

■ 

1 ,379353E»15 

Chemical 

concentration  of  eCRON 

CXIOE 

C-ARGE 

• 

0 ,?P00PPE*0P 

IS 

!>,000E*00 

X 

OF 

total 

SILICON 

Charge 

s 

1 ,884S53E*15 

IS 

100.0 

X 

OF 

total 

total 

Charge 

■ 

1 .884K53EM5 

IS 

4.006F*05 

X 

OF 

initial 

Initial 

C-aRGE 

■ 

4 , 702 1 88E • 1 1 

che-ical 

CCNCEntbation  OF  arsenic 

OXIDE 

charge 

■ 

0,00PP0PE«00 

IS 

0,000E*00 

X 

OF 

total 

SILICON 

C-aRGE 

a 

1 .379S2SE.15 

IS 

100.0 

X 

OF 

total 

total 

C-ARGE 

a 

1 .379828E*1S 

IS 

103.0 

X 

OF 

Initial 

Initial 

Charge 

8 

1 ,379827E*15 

i 


I 


STANFORD  KITCHIP  PROCESS 
BASE  DIFFUSION 
STEP  * 9 


S/3/78 


OXIDATION  IN  DRV  OXYGEN 

TOTAL  STEP  TIME  * 43. a MINUTES 

INITIAL  TEMPERATURE  ■ 10F0.000  DECREES  C. 

OXIDE  THICKNESS  • 4.7P21E-02  MICRONS 

LINEAR  OXIDE  GROWTH  RATE  ■ 7.994831E-04  MICRONS/MINUTE 
PARABOLIC  OXIDE  GROWTH  RATE  ■ 1.007331E-04  MICRONSt2/MINUTE 
OXIDE  GROWTH  PRESSURE  ■ 1.P0000  ATMOSPHERES 

1 OXIOE  I SILICON  I l SURFACE  1 

1 DIFFUSION  1 DIFFUSION  J SEGREGATION  l TRANSPORT 

1 COEFFICIENT  l COEFFICIENT  1 COEFFICIENT  1 COEFFICIENT  l 


2.01360E-06  l 2.P 3949E-04  | 


.34974 


ARSENIC 


4.60947E-10  1 1.01692E-03  1 10,0000 


2.32620E-0S  1 
1.19298E-02  1 


SURFACE  CONCENTRATION  • 2.542623E*19  ATOMS/CMt3 


JUNCTION  DEPTH 
1,03444  MIC 


sheet  resistance 


MICRONS  I 73.6273 
1 24.0069 


OHMS/SQUARE 

OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIOE 

charge 

■ 

2.616383EM4 

IS 

9.47 

X 

OF 

total 

SILICON 

CHARGE 

• 

2 ,6936286*1 3 

IS 

90.5 

X 

OF 

total 

TOTAL 

CHARGE 

• 

2 , 9734B6E* 1 5 

IS 

91.2 

X 

OF 

initial 

INITIAL 

CHARGE 

■ 

3.263640E*15 

CHEMICAL 

CONCENTRATION  of  BORON 

OXIOE 

CHARGE 

• 

2 , 81 0390E* 1 4 

IS 

17.6 

X 

OF 

total 

SILICON 

CHARGE 

■ 

1 .31S187EM3 

IS 

82.4 

X 

OF 

TOTAL 

TOTAL 

CHARGE 

■ 

1 .396826EM3 

IS 

04.7 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

I .6046536*15 

Chemical 

CONCENTRATION  OF  ARSENIC 

OXIOE 

CHARGE 

S 

7.2611 106*08 

IS 

3.277E-05 

X 

OF 

total 

SILICON 

CHARGE 

■ 

1 . 3790206*13 

IS 

100.0 

X 

OF 

total 

total 

CHARGE 

• 

1.3790296*13 

IS 

100. 

X 

OF 

initial 

initial 

CHARGE 

• 

1 ,379626E*1 3 

*v  .» 

i . , . 


1 


I 

p 

I 


i 

I 

*1 

: 


I 


i 


STANFORD  K ITCH  IP  PROCESS  6/5/70 

BASE  DIFFUSION 
STEP  • 10 


OXIDATION  IN  WET  OX VSEN 

TOTAL  STEP  TIHE  « 60.0  MINUTES 

INITIAL  TEMPERATURE  ■ 1000,000  DEGREES  C. 

OXIDE  THICKNESS  • , 4235  MICRONS 


LINEAR 

OXIDE 

GROWTH  RATE 

■ 

1.472099E.02  MICRONS/MINUTE 

PARABOLIC 

OXIOE 

GROWTH  RATE 

■ 

5.149767E-03  MICR0NSt2/MlNuTE 

OXIDE  GROWTH  PRESSURE 

■ 

,842109  ATMOSPHERES 

1 

OXIDE 

1 

SILICON  1 

1 

surface 

1 

1 

DIFFUSION 

1 

DIFFUSION  1 

SEGREGATION  { 

TRANSPORT 

» 

COEFFICIENT 

1 

COEFFICIENT  1 

COEFFICIENT  1 

coefficient 

1 

BORON 

l 

2.01306E.0S 

1 

2.03949E-04  l 

.34974  1 

2.92626E*03 

1 

ARSENIC 

l 

4 . 60947E* 1 0 

1 

1.01692E-05  1 

10.0000  1 

1 .1929BE-02 

1 

SURFACE  CONCENTRATION  ■ 3,871102E«18  ATOMS/CMtS 

JUNCTION  DEPTH  1 SHEET  RESISTANCE 



1.20266  MICRONS  1 127.610  OHMS/SQUARE 

1 24,6019  OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


oxide 

CHARGE 

9 

9.046204E+14 

IS 

30.6 

X 

OF 

total 

SILICON 

CHARGE 

9 

2 .028090E+ 1 9 

IS 

69.2 

X 

OF 

total 

total 

CHARGE 

9 

2.932671EM5 

IS 

96.6 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

9 

2.975466EM5 

CHEMICAL  CONCENTRATION  OF  BORON 


OXIOE 

CHARGE 

a 

9 , 0463  B8E  + 1 4 

IS 

96.2 

X 

OF 

total 

SILICON 

CHARGE 

a 

6 . 49469  IE*  1 4 

IS 

41.6 

X 

OF 

total 

TOTAL 

CHARGE 

a 

1 .994104E*19 

IS 

97.3 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

a 

1 .59662SEM5 

CHEMICAL  CONCENTRATION  OF  ARSENIC 


OXIOE 

CHARGE 

a 

1 . 84S629E* 1 0 

IS 

1 .336E-03 

X 

OF 

total 

SILICON 

CHARGE 

a 

1 ,379007E*1  9 

IS 

100.0 

X 

OF 

total 

TOTAL 

CHARGE 

a 

1 ,379625E*19 

IS 

100.0 

X 

OF 

INITIAL 

initial 

CHARGE 

a 

1.379029EM9 
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STANFORD  KITCMIP  PROCESS 
BASE  OIFFUSION 
STEP  a 11 


6/9/7# 


OXIDATION  IN  DRV  OXYGEN 

TOTAL  STEP  TINE  ■ 9.0  minutes 

INITIAL  TEMPERATURE  ■ 100P.000  DEGREES  C. 

OXIDE  THICKNESS  ■ ,4243  MICRONS 

LINEAR  OXIDE  GROWTH  RATE  ■ 7.889412E.04  MICRONS/MINUTE 

PARABOLIC  OXIDE  GROWTH  RATE  ■ 1.795914E.04  MICRONST2/MINUTC 

OXIDE  GROWTH  PRESSURE  ■ 1.00P00  ATMOSPHERES 

1 OXIDE  1 SILICON  1 1 SURFACE  1 

1 DIFFUSION  1 DIFFUSION  J SEGREGATION  J TRANSPORT 

l COEFFICIENT  1 COEFFICIENT  1 COEFFICIENT  1 COEFFICIENT  J 


BORON 

ARSENIC 


t 2.01368E-08  1 2.03949E-04  1 


.34974 


4.60947E-10  1 1.01692E-09  1 10.0000 


SURFACE  CONCENTRATION  ■ 9,871 862E+1 8 AT0MS/CMt3 


I 2.52626E-03  I 


1.1929BE-02  i 


JUNCTION  DEPTH 
1.22169  MIC 


SHEET  RESISTANCE 


MICRONS  1 128.319 

1 24.7969 


.OHMS/SQUARE 

OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIOE 

CHARGE  ■ 

9.127785EM4 

IS 

31.1 

X 

OF 

TOTAL 

SILICON 

CHARGE  ■ 

2. 01 97  51 E*  19 

IS 

68.9 

X 

OF 

TOTAL 

TOTAL 

CHARGE  ■ 

2.932929E+19 

IS 

100.0 

X 

OF 

INITIAL 

initial 

CHARGE  a 

2.932671EMS 

CHEMICAL 

, CONCENTRATION  OF  BORON 

OXIOE 

CHARGE  a 

9.127969EM4 

IS 

98.7 

X 

OF 

total 

SILICON 

CHARGE  ■ 

6.41  1452EM4 

IS 

41.3 

X 

OF 

TOTAL 

TOTAL 

CHARGE  ■ 

1 .553942EM5 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE  a 

1 .994104EM9 

1 

CHEMICAL 

CONCENTRATION  OF  ARSENIC 

OXIDE 

CHARGE  ■ 

1 .837133EM0 

IS 

1.331E-03 

X 

OF 

total 

SILICON 

CHARGE  • 

1 .379832EM9 

IS 

100.0 

X 

OF 

total 

TOTAL 

CHARGE  ■ 

1 , 379890E* 1 5 

IS 

100. 

X 

OF 

initial 

initial 

CHARGE  ■ 

1 ,379829E*19 

STANFORD  KITCHIP  PROCESS 
BASE  DIFFUSION 
STEP  « 12 


6/5/78 


neutral  ambient  drive-in 

total  STEP  TIME  ■ 10.0  MINUTES 


INITIAL 

TEMPERATURE  ■ 1 000.000  DEGREES 

c. 

OXIDE  THICKNESS 

■ .4243  MICRONS 

i 

OXIDE  1 

SILICON 

i 

1 

SURFACE 

1 

1 

OIFFUS ION  t 

DIFFUSION 

l 

SEGREGATION 

1 

TRANSPORT 

1 

COEFFICIENT  1 

COEFFICIENT 

1 

COEFFICIENT 

1 

COEFFICIENT 

» 

BORON 

1 

2,01 3S8E-08  1 

9.16197E-09 

» 

.34974 

1 

2.32626E-03 

1 

ARSENIC 

1 

4.S0947E-10  1 

1 . 0 1 892E-05 

1 

10.0000 

1 

1.1 929BE-02 

1 

SURFACE  CONCENTRATION  ■ 6,B91920E*18  ATOMS/CMt3 


Junction  oepth  i sheet  resistance 

- ------ 

1.24112  MICRONS  1 128, SBS  OHMS/SQUARE 

1 24,7988  OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE 

CHARGE 

■ 

9 . 1020386*1  4 

IS 

31.2 

X 

OF 

total 

SILICON 

CHARGE 

■ 

2,0102006*15 

IS 

08.8 

X 

OF 

total 

total 

CHARGE 

■ 

2,9324046*15 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

2.9325296*15 

CHEMICAL  CONCENTRATION  OF  BORON 


OXIDE 

CHARGE 

■ 

9 , 1 6282 1 E* 1 4 

IS 

59.0 

X 

OF 

total 

SILICON 

CHARGE 

■ 

0.3759806*14 

IS 

41.0 

X 

OF 

TOTAL 

TOTAL 

CHARGE 

■ 

1 ,553881E*15 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

charge 

■ 

1.5539426*15 

chemical  concentration  of  arsenic 


OXIOE 

CHARGE 

« 

1 ,B29183E*10 

IS 

1 , 320E-03 

X 

OF 

total 

SILICON 

charge 

■ 

1 .379B32E*! 5 

IS 

100.0 

X 

OF 

total 

total 

CHARGE 

■ 

1 ,379850E*15 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

charge 

■ 

1 ,379850E*15 
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STANFORD  KITCHIP  PROCESS 


S/S/78 


ETCH  STEP 

ETCH  TEMPERATURE  • 23.8  DEGREES  C. 

ETCH  TIME  » .8  MINUTES 

ETCH  RATE  • ,8888  MICRONS/MINUTE 

OXIDE  THICKNESS  ■ 8 ,0000E*08  MICRONS 

SURFACE  CONCENTRATION  ■ 8,8919288*18  ATOMS/CMtS 

JUNCTION  DEPTH  1 SHEET  RESISTANCE 



1.23930  MICRONS  | 128.389  OHMS/SQUARE 

I 24.7996  OHMS/SOUARE 


NET  ACTIVE  CONCENTRATION 


OXIOE 

CHARGE 

■ 

8 ,888888E*08 

IS 

8.888E*88 

X 

OF 

TOTAL 

SILICON 

CHARGE 

• 

2.0!6]S1E*19 

IS 

188.8 

X 

OF 

total 

total 

CHARGE 

■ 

2.016181E*19 

IS 

88.8 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

2.932464E*15 

chemical 

, CONCENTRATION  OF  BORON 

OXIOE 

CHARGE 

a 

0 ,8088886*08 

IS 

8.0008*88 

X 

OF 

TOTAL 

SILICON 

CHARGE 

a 

6.3759746*14 

IS 

188.8 

X 

OF 

total 

total 

CHARGE 

a 

6.375974E*14 

IS 

A 1 . 8 

X 

OF 

INITIAL 

initial 

CHARGE 

a 

1 ,553881E*15 

chemical 

CONCENTRATION  of  arsenic 

OXIDE 

CHARGE 

a 

8,8888B8E*08 

IS 

0 . 000E*80 

X 

OF 

total 

SILICON 

CHARGE 

a 

1 .3798148*15 

IS 

100.0 

X 

OF 

total 

total 

CHARGE 

a 

1 ,379814E*15 

IS 

100.0 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

• 

1 . 379858E*  1 5 
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STANFORD  RITCMIR  PROCESS  6/3/7» 

BASE  DIFFUSION 

J STEP  • 13  TINE  » .0  MINUTES. 

DEPTH  { CONCENTRATION  (U06  AT0MS/CC1 


V"\ 


STANPORO  KITCHIP  PROCESS 
EMITTER  DIFFUSION 
STEP  • 14 


6/9/76 


GASEOUS  PREDEPOSITION 

TOTAL  STEP  TIME  ■ 33.6  MINUTES 

INITIAL  TEMPERATURE  « 1029,00  DEGREES  C, 

OXIDE  THICKNESS  ■ 0.0000E400  MICRONS 

PREOEPOS JTION  IMPURITY  ■ PHOSPHORUS 

GAS  CONC.  AT  INTERFACE  ■ t,190000E*2l  ATOMS/CC 


l OXIOE 
1 DIFFUSION 
1 COEFFICIENT 


1 SILICON  I 

1 DIFFUSION  1 
1 COEFFICIENT  l 


1 surface  1 

SEGREGATION  1 TRANSPORT 
COEFFICIENT  l COEFFICIENT  l 


BORON 

• ••1 

l 

3.74197E-06  1 

1.67164E-04  1 

.41393 

1*91 

1 

3.90494E-03 

1 

ARSENIC 

l 

1.08752E-09  | 

2.08126E-05  1 

10,0000 

1 

1,691 77E.02 

1 

PHOSPHORUS 

1 

1.17679E-06  1 

2.967091.04  1 

10.0000 

l 

1.0000 

1 

SURFACE  CONCENTRATION  ■ .3.390586E*20  AT0MS/CMt3 


JUNCTION  DEPTH 


SHEET  RESISTANCE 


1,06376  MICRONS  i 6.73261 

1,59913  MICRONS  1 26066.9 

i 24.7797 


NET  ACTIVE  CONCENTRATION 


OHMS/SQUARE 

OMMS/SQUARE 

OHMS/SQUARE 


OXIOE  CHARGE  ■ 
SILICON  CHARGE  » 
TOTAL  CHARGE  ■ 
INITIAL  CHARGE  ■ 


0 .0000P0E*00  IS 
1.370081E*16  IS 
1.370001EM6  IS 
2.016181E«!9 


0.000E*00  X OF  TOTAL 

100.0  X OF  TOTAL 

660.  X Of  INITIAL 


CHEMICAL  CONCENTRATION  OF  BORON 


OXIOE  charge 

SILICON  CHARGE 
TOTAL  CHARGE 
INITIAL  CHARGE 


0.000P00E+00  IS 
6.1251 95E* 1 4 IS 
6 , 1291 95E*1 4 IS 
6.375974E414 


0 . 000E*00  X OF  TOTAL 

100.0  X OF  TOTAL 

96.1  X OF  INITIAL 


CHEMICAL  CONCENTRATION  OF  ARSENIC 


OXIOE  CHARGE 
SILICON  CHARGF 
TOTAL  CHARGE 
INITIAL  CHARGE 


0.000000 E*00  IS 
1.379605E*1S  IS 
1.379009EMS  IS 
1.379P14EM9 


0,000E*00  X OF  TOTAL 

100.0  X OF  TOTAL 

100.0  X OF  initial 


CHEMICAL  CONCENTRATION  OF  PHOSPHORUS 


OXIOE  CHARGE 
SILICON  CHARGE 


0.000P00E*0P 
2 ,973232(416 


TOTAL  CHARGE  • 
INITIAL  CHARGE  ■ 


2 . 97  3232E*  1 6 
0.000000E400 


IS 

IS 


0.000E*00  X OF  TOTAL 
100.0  X OF  TOTAL 


IS  0.000t>00  X OF  INITIAL 
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STANFORD  KITCHIP  PROCESS  8/3/78 

emitter  diffusion 


1 STEP  ■ 

14 

TIME  • 33.0  MINUTES. 

DEPTH 

1 

l 

CONCENTRATION 

(LOG  ATOMS/CC) 

CUM) 

A A 

l 

14 

15 

IS  17 

IB  19 

90 

21 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

l 

1 

l 

l 

1 

1 

1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

i 1 

1 

1 

1 

l 

1 

1 

l 

l 

l 

l 

1 

1 

1 

1 

l 

1 

1 

1 

1 

1 

1 1 

1 1 

J * 

l 

1 

1 

1 1 

1 1 

1 * 

1 

1 

1 

1 1 

1 1 

| * 

1 

1 

1 

1 1 

1 1 

J * 

1 

1 

1 

1 1 

1 1 

• 1 

1 

1 

1 

1 1 

1 1 

* 1 

1 

1 

l 

1 1 

S l 

* 1 

1 

1 

1 

1 1 

l * 

1 

t 

1 

1 

1 1 

1 * 1 

1 

l 

l 

l 

J • 1 

1 1 

l 

1 

1 

1 

J * 1 

1 1 

1 

1 

1 

1 

i * i 

1 1 

1 

1 

1 

1 

A 1 1 

1 1 

i 

1 

1 

i * 

1 1 

1 1 

1 

l 

| * 

l 

1 

• 1 

1 1 

1 1 

1 1 

1 1 

1 

1 

1 

1 

1 

• 

1 1 

1 1 

1 

1 

1 

• 

1 1 

1 1 

1 

1 

1 

• 

1 1 

1 1 

1 

1 

1 

• 

1 1 

1 1 

1 

1 

1 

• 

1 1 

1 1 

l 

1 

1 

• 

1 1 

1 1 

1 

1 

1 

• 

1 1 

1 1 

l 

1 

1 

• 

1 I 

1 1 

1 

1 

1 

• 

l 1 

1 1 

1 

1 

1 

* 

1 1 

1 1 

1 

1 

1 

1* 

1 1 

1 1 

l 

1 

1 

1 

1 * 

1 

I i 

# 1 i 

1 1 

1 1 

» 

l 

1 

1 

1 

1 

* 1 1 

1 l 

1 

1 

1 

1 

• 1 1 

1 1 

1 

1 

1 

1 

1 • 1 

1 1 

1 

1 

1 

1 

1 • 1 

1 1 

1 

1 

1 

1 

1 * 1 

1 1 

1 

l 

l 

1 

1 • 1 

1 1 

1 

1 

1 

1 

l • 

1 1 

1 

1 

4,08  * 
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STANFORD  KJTCHTP  PROCESS 
EMITTER  DIFFUSION 
STEP  • 15 


NEUTRAL  AMBIENT 

DRIVE-IN 

TOTAL  STEP 

TIME 

■ 10.0  MINUTES 

INITIAL  TEMPERATURE  » 900.000  DEGREES  C 

OXIDE  THICKNESS 

■ 0 ,0000E*00  MICRONS 

1 

OXIDE  1 

SILICON  1 

1 

surface 

1 

i 

1 

DIFFUSION  1 

DIFFUSION  1 

SEGREGATION 

1 

TRANSPORT 

l 

COEFFICIENT  i 

COEFFICIENT  1 

COEFFICIENT 

1 

COEFFICIENT 

1 

BORON 

1 

i 

1.29677E-09  1 

6.39B43E-06  1 

.16589 

1 

3.67548E-04 

1 

ARSENIC 

l 

1.03189E-11  1 

8 , 27 1 44C-0  7 | 

10.0000 

1 

2.54185E-03 

1 

PHOSPHORUS 

1 

4,1 9S32E-0S  1 

4.36609E-06  t 

10.0000 

1 

0.00000E400 

1 

SURFACE  CONCENTRATION  ■ -3.322567E*?0  ATDMS/CMT3 


JUNCTION  DEPTH  1 SHEET  RESISTANCE 


1 .08766  MICRONS  1 6,71579  OHMS/SOUARE 

1.56485  MICRONS  J 26019.8  0**575 OUARE 

t 24,7739  OHMS/SQUARE 


NET  ACTIVE  CONCENTRATION 


OXIDE 

Charge 

■ 

0 .00P000fe*00 

IS 

0 , 00PE*0B  X 

OF 

total 

SILICON 

CHARGE 

■ 

1 .3729B4E*16 

IS 

100.0  X 

OF 

total 

TOTAL 

Charge 

■ 

1 ,3729B4E*16 

IS 

100.  X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

1 ,370081E*16 

Chemical 

CONCENTRATION  OF  boron 

oxinE 

CHARGE 

■ 

0 ,000000E*00 

TS 

0 ,000E»00 

X 

OF 

total 

SILICON 

CHARGE 

• 

0.12! 1 20E* 1 4 

IS 

100.0 

X 

OF 

total 

total 

CHARGE 

■ 

6. 121 120E*1 4 

IS 

99.9 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

■ 

6.125195EM4 

CHEMICAL 

CONCENTRATION  OF  ARSENIC 

OXIDE 

CHARGE 

■ 

0 ,000000E*00 

IS 

0 . 000E*00 

X 

OF 

total 

SILICON 

CHARGE 

■ 

1 ,379B04E*15 

IS 

100.0 

X 

OF 

total 

total 

CHARGE 

■ 

1 ,37980 4E*15 

IS 

100.0 

X 

OF 

initial 

initial 

charge 

9 

1 .379B09EM5 

chemical 

CONCENTRATION  OF  PHOSPHORUS 

oxide 

CHARGE 

■ 

0 ,000000E»00 

IS 

0,000E«00 

X 

OF 

total 

SILICON 

Charge 

t 

2.97323JEM6 

IS 

100.0 

X 

OF 

total 

total 

CHARGE 

■ 

2.973231EM6 

19 

100.0 

X 

OF 

INITIAL 

initial 

CHARGE 

• 

2.973232EM6 
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STANFORD  KITCHJP  PROCESS 
EMITTER  DIFFUSION 
STEP  M 18 


6/3/78 


f i 

OXIDATION  IN  NET  OXVGEN 


1 • 

TOTAL  STEP  TIME  » 

IP. P 

MINUTES 

1 

INITIAL  TEMPERATURE  ■ 

9PP 

.P0P 

DEGREES 

c. 

OXIDE  thickness  ■ 

• 

1777 

MICRONS 

l'  ■*. 

LINEAR  OXIDE  GROWTH 

RATE 

• 

3.628432E.P2 

MICRONS/MINUTE 

PARABOLIC  OXIOE  GROWTH 

RATE 

■ 

4,381 483E-03 

MICR0NST2/MINUTE 

OXIDE  GROWTH  PRESSURE 

■ 

,842105 

ATMOSPHERES 

1 OXIDE 

i 

SILICON 

1 1 

SURFACE  1 

1 DIFFUSION 

i 

DIFFUSION 

1 SEGREGATION  1 

TRANSPORT 

1 COEFFICIENT 

1 

coefficient 

1 COEFFICIENT  1 

COEFFICIENT  1 

BORON 

• m m a 

1 

1.29677E-09  1 

3,661436-05  1 

,18589 

1 

3.67548E-04 

1 

ARSENIC 

l 

1.P3185E-11  1 

4.27144E-B7  1 

10.0000 

l 

2.34185E-03 

l 

PHOSPHORUS 

1 

4.19632E-08  1 

7.86003E-06  1 

10.0000 

l 

0.00000E*00 

1 

SURFACE  CONCENTRATION  ■ -5 . 1 97843E*20  ATOMS/CMtS 


JUNCTION  OEPTH  1 SHEET  RESISTANCE 


1.00799  MICRONS  l 8.25491  OHMS/SOUARE 

1 .49109  MICRONS  1 24321.2  OHMS/SQUARE 

1 24,7715  OHMS/SOUARE 


NET  ACTIVE  CONCENTRATION 


: 


OXIDE  CHARGE  ■ 

3 , 463763E*  1 5 

IS 

23.4 

X 

OF 

total 

SILICON  CHARGE  « 

1 . 135358E*1« 

IS 

76.6 

X 

OF 

total 

^ total  CHARGE  ■ 

1 ,481735E*16 

IS 

108. 

X 

OF 

INITIAL 

i? 

INITIAL  CHARGE  ■ 

1 ,372984E*16 

r- 

CHEMICAL  CONCENTRATION  OF  BORON 

t 

OXIDE  CHARGE  ■ 

4 . 604792E* 1 3 

IS 

7.53 

X 

OF 

total 

SILICON  CHARGE  « 

5, 63921 3E *14 

IS 

92.5 

X 

OF 

total 

\\ 

TOTAL  CHARGE  ■ 

6.099«92E*t4 

IS 

99.6 

X 

OF 

INITIAL 

INITIAL  CHARGE  • 

6,1211 20E*  1 4 

R . 1 

» 

CHEMICAL  CONCENTRATION  OF  ARSENIC 

OXIDE  CHARGE  » 

3.248391E*09 

IS 

2.354E-04 

X 

OF 

total 

SILICON  CHARGE  ■ 

1 ,379806E*13 

IS 

100.0 

X 

OF 

total 

TOTAL  CHARGE  ■ 

1 ,379809E*15 

IS 

100. 

X 

OF 

INITIAL 

INITIAL  CHARGE  ■ 

1 . 379804E*15 

CHEMICAL  CONCENTRATION  OF  PHOSPHORUS 


U 


OXIDE 

CHARGE 

■ 

6.676233E*13 

IS 

22.3 

X 

OF 

total 

SILICON 

CHARGE 

■ 

2.322669E*16 

IS 

77.7 

X 

OF 

total 

TOTAL 

CHARGE 

a 

2 , 990994E*  1 8 

IS 

101. 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

a 

2.973231E*10 
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STANFORD  KITCHJP  PROCESS 
EMITTER  DIFFUSION 
STEP  P 17 


6/3/70 


NEUTRAL  AMBIENT  DRIVE-IN 

TOTAL  STEP  TIME  ■ 10.0  MINUTES 

INITIAL  TEMPERATURE  ■ 900.000  DECREES  C. 

OXIDE  THICKNESS  ■ .1777  MICRONS 


l 

OXIDE  l 

SILICON  l 

l 

SURFACE 

1 

1 

DIFFUSION  1 

DIFFUSION  1 

segregation 

1 

TRANSPORT 

1 

COEFFICIENT  1 

COEFFICIENT  1 

COEFFICIENT 

1 

COEFFICIENT 

1 

BORON 

1 

1.29677E-09  l 

S.39843E-06  | 

.16589 

1 

3.67546E-04 

1 

ARSENIC 

1 

1.03185E-11  1 

A.27144E-07  1 

10.0000 

l 

2.54185E-03 

1 

PHOSPHORUS 

l 

4 , 1 9632E.0S  1 

4 , 36089E.06  1 

10.0000 

1 

0,00000E*00 

l 

SURFACE  CONCENTRATION  ■ -4,416993E*20  AT0MS/CMT3 


JUNCTION  DEPTH  1 SHEET  RESISTANCE 

....... ------------- 1----- ........ ......... 


1.00908  MICRONS  1 8.20940 

OHMS/SQUARE 

1.49193  MICRONS  1 24343.0 

OHMS/SOUARE 

1 24.7713 

OHMS/SQUARE 

NET  ACTIVE  CONCENTRATION 

OXIDE  CHARGE  ■ 3 . 40324SE*  1 5 

IS 

23.0 

X 

OF 

TOTAL 

SILICON  CHARGE  ■ 1 . 1 41 ?75E* 1 6 

IS 

77.0 

X 

OF 

total 

TOTAL  CHARGE  ■ 1.401600EM6 

INITIAL  CHARGE  ■ 1 . 48 1 735E*  1 6 

IS 

100.0 

X 

OF 

INITIAL 

CHEMICAL 

CONCENTRATION  OF  80R0N 

OXIDE 

CHARGE 

■ 

4 . 670280E* 1 3 

IS 

7,68 

X 

OF 

total 

SILICON 

CHARGE 

■ 

5 , 62945BE*  1 4 

IS 

92.3 

X 

OF 

total 

total 

CHARGE 

■ 

6.096486EM4 

IS 

99.9 

X 

OF 

INITIAL 

INITIAL 

CHARGE 

a 

6 ,099692E*14 

CHEMICAL 

CONCENTRATION  OF  ARSENIC 

OXIDE 

CHARGE 

a 

3,1 66797E*09 

IS 

2.295E-04 

X 

OF 

TOTAL 

SILICON 

CHARGE 

■ 

1 .379806EM5 

IS 

100.0 

X 

OF 

TOTAL 

total 

CHARGE 

■ 

1 .379809EM5 

IS 

100.0 

X 

OF 

INITIAL 

initial 

CHARGE 

a 

1 .3790P9EM5 

CHEMICAL 

CONCENTRATION  OF  PHOSPHORUS 

OXIDE 

CHARGE 

a 

6 . 490390E* 1 5 

IS 

21.7 

X 

OF 

TOTAL 

SILICON 

CHARGE 

a 

2 , 34 1 234E*  1 6 

IS 

78.3 

X 

OF 

total 

TOTAL 

CHARGE 

a 

2 , 990293E*  1 6 

IS 

100.0 

X 

OF 

INI  T I 

INITIAL  CHARGE  ■ 2 , 990294E*  1 6 

-114- 


1 STEP  • 17  TIME  ■ 10,0  MINUTES, 


OEPTM  1 CONCENTRATION  (LOG  ATOMS/CC) 
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